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Der Einflu6 einer durch Elektronen bewirkten Gitterdeformation 
auf optische und elektrische Eigenschaften von Kristallen’) 


Von S. 1. PEKAR 


1. Einleitung 


Dieser Artikel ist der Theorie der nichtmetallischen Kristalle, d. h. der Iso- 
latoren und Halbleiter gewidmet. Die moderne Theorie der Halbleiter laiBt 
sich in zwei Teile gliedern: die phainomenologische Theorie und die Mikro- 
theorie der Elektronenzustiinde im Kristall. Die erstere benutzt beispielsweise 
das OHMsche Gesetz oder die allgemeinere Gleichung fiir elektrische Leitung 
und Diffusion der ,,Ladungstrager“ 

dV dN 


= b D 
4) ebN a e aa (1) 


dabei ist N die Konzentration der Ladungstriger, e ihre Ladung, 6 ihre 
_ Beweglichkeit, D der Diffusionskoeffizient, V das Potential des elektrischen 
Feldes und J die Stromdichte. Ferner ist fiir die phinomenologische Theorie 
die Potssonsche Gleichung kennzeichnend: 


AV =—4zo(N), (2) 


hier ist o9(N) die Raumladungsdichte der Ladungstriger. Zur phainomeno- 
logischen Theorie muB man auch die Temperaturabhiingigkeit der Kon- 
zentration der Ladungstriger rechnen; sie hat die Form 
fy le 

Wiz ha PRE» (3) 
wobei W die thermische Dissoziationsenergie der Ladungstrager ist (Fall des 
thermischen Gleichgewichts), und wobei der Ausdruck fiir die Anzahl der 
Rekombinationen von Ladungstriigern in der Volumeneinheit und der Sekunde 


pNN, (4) 


lautet; hierbei ist 6 der Rekombinationskoeffizient und NV, die Konzentration 
der Rekombinationszentren. 

Die oben angegebenen Gleichungen kann man insofern als phinomenologisch 
bezeichnen, als bei ihrer Begriindung der Mechanismus der Mikroprozesse 
vernachliassigt wird. Dies ergibt sich schon daraus, dai diese Gleichungen 


1) Ungekiirzte Ubersetzung des in Uspechi Fiz. Nauk 50, 197, 1953 erschienenen Artikels. 
25 Zeitschrift ,,Fortschritte der Physik‘‘ 
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gebunden sind. Die meisten Erscheinungen werden in dieser nullten Naherung 
behandelt. Nur in wenigen Fiillen wird als kleine Stérung auch die Wechsel- 
wirkung zwischen Elektronen und Atomschwingungen eingefiihrt. Dies wird 
beispielsweise erforderlich beim Problem der Beweglichkeit und der freien 
Weglinge der Leitungselektronen, denn diese GréBen werden unendlich, 
wenn man sich auf die oben erwihnte nullte Niherung beschrankt. 
Hinsichtlich der ersten der erwihnten vereinfachenden Annahmen mu man 
betonen, da8 in Metallen die Ersetzung des Mehrelektronenproblems durch 
ein Einelektronenproblem nicht hinreichend gerechtfertigt werden konnte, 
Die gewohnlich in den Lehrbiichern angefiihrten Begriindungen fiir diese 
Ersetzung kénnen nicht als korrekt angesehen werden [PEKAR (1948a, 
1952a)]. Auch bei Behandlung der Kigenelektronen’) eines Isolators 1aBt 
sich diese Ersetzung nicht rechtfertigen. Handelt es sich jedoch um die 
Leitungselektronen eines Halbleiters oder Isolators oder auch um die Elek- 
tronen eines lokalen ,,Zentrums‘‘, die mit den Kernen weitaus schwéacher 
gekoppelt sind als die Kigenelektronen im Isolator oder Halbleiter, so ist 
die bekannte adiabatische Naherung brauchbar: man nimmt an, daf die 
schwach gebundenen Leitungselektronen (oder die Elektronen des lokalen 
Zentrums) sich hinreichend langsam bewegen, so daf der Zustand der stark 
gebundenen Higenelektronen des Kristalls dieser Bewegung der Leitungs- 
elektronen adiabatisch folgen kann. In diesem Fall braucht man gemaB der 
adiabatischen Naherung nur die Bewegung der schwach gebundenen lang- 
samen Elektronen zu behandeln, und das Vorhandensein der Eigenelektronen 
ist nur von HinfluB auf die potentielle Energie fiir diese langsamen Elektronen 
[PEKAR (1948a, 1951, 1952a, 1952b)]. 

In nichtmetallischen Kristallen 1aéBt sich also in bezug auf die schwach ge- 
bundenen Elektronen die erste Voraussetzung zum Teil rechtfertigen: die 
CouLomBsche Wechselwirkung der schwach gebundenen Elektronen mit den 
EKigenelektronen des Kristalls 1aBt sich durch Einfiihrung eines zweckmiabig 
gewihlten auBeren Potentials beriicksichtigen. Im Fall der Leitungselektronen 
im idealen Isolator ist dieses Feld periodisch [PEKAR (1951)]|. Weiter unten 
wird eben von diesen schwach gebundenen Elektronen die Rede sein, auf die 
die erste vereinfachende Voraussetzung anwendbar ist. Die meisten der unten 
angegebenen Ergebnisse kann man auch gewinnen, indem man an Stelle der 
ersten Voraussetzung die Niiherung von HEITLER-LONDON-HEISENBERG 
(HLH)?) benutzt. Sie besteht darin, da’ man die Mehrelektronen-Wellen- 
funktion des Kristalls als antisymmetrisches Produkt der Wellenfunktionen 
der einzelnen Atome oder Molekiile konstruiert. Hierbei werden die mole- 
kularen Wellenfunktionen als bekannt, die Wechselwirkung zwischen den 
Molekiilen als schwach angesehen. 

Wenn wir uns nun der zweiten vereinfachenden Voraussetzung zuwenden, 
haben wir zu beachten, daf die Aufteilung des Kristalls in zwei kaum mit- 
einander wechselwirkende Untersysteme, niimlich die Elektronen und die 


1) Anm. d. dtsch. Red.: ,,Higenelektronen‘* sind die im Idealkristall vorhandenen 
Elektronen, ,,Leitungselektronen‘* werden z. B. beim Einbau von weiteren Atomen in 
den Kristall gebracht, falls diese dissoziieren (Donatoren), 

*) Siehe SOMMERFELD und BEeTHE (1933), Teil IIT und PrKar (1948). 
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Atomschwingungen, die Méglichkeit einer auch nur qualitativen theoretischen 
Deutung vieler Erscheinungen von vornherein ausschlieBt. Man wei bei- 
spielsweise, dal die optischen Ubergiinge der Elektronen im Kristall im all- 
gemeinen mit einer Anregung von Gitterschwingungen, d. h. mit einer Warme- 
erzeugung verbunden sind. Die hierbei freiwerdende Wirmemenge liBt sich 
als Differenz zwischen der experimentell beobachtbaren Energie iw des 
optischen Uberganges und der Energie W des entsprechenden thermischen 
Elektroneniiberganges berechnen. hq@ ist oftmals anderthalb mal so grob 
wie |W (siehe unten). Es verwandelt sich also ein bedeutender Teil der Energie 
des absorbierten Lichtquants, d. h. eine Energie von der GréBenordnung 1 eV 
in Warme. Diese Tatsache zeugt von der starken Wechselwirkung zwischen 
den Elektronen und den Atomschwingungen und liBt sich in einer Theorie, 
die diese Wechselwirkung als schwache Stérung auffaBt, nicht erkliren. Eine 
solche Theorie ist auch nicht imstande, die hiufig beobachtete groBe STOKES- 
sche Verschiebung in Kristallen (den groBen Abstand zwischen den Absorp- 
tions- und den Lumineszenzbanden, die den Elektroneniibergingen in den 
verschiedenen Richtungen entsprechen) zu erkliiren sowie die groBe Halb- 
wertsbreite der Absorptions- und Lumineszenzbanden, die manchmal die 
GréBenordnung | eV erreicht, wenn das Elektron von einem diskreten auf 
einen diskreten Term tibergeht, ferner Form und Temperaturabhingigkeit 
dieser Banden u. dgl. 

Ferner ist eine Theorie, die die Wechselwirkung zwischen Elektronen und 
Atomschwingungen als kleine Stérung ansieht, nicht in der Lage, die strahlungs- 
losen thermischen Elektroneniibergange im Kristall zu deuten; je mehr 
Schwingungsquanten namlich bei einem Ubergang erzeugt (absorbiert) werden, 
in umso héherer Ordnung verschwindet in dieser Theorie die Ubergangs- 
wahrscheinlichkeit. Uberginge, bei denen einige Dutzend Schwingungs- 
quanten erzeugt werden, besitzen nach dieser Theorie nur eine verschwindend 
kleine, zur Erklarung der beobachteten Tatsachen vdéllig unzureichende 
Wahrscheinlichkeit. Diese Beispiele erschépfen natiirlich bei weitem nicht 
die Schwierigkeiten jeder Theorie, die sich auf die zweite vereinfachende 
Annahme stiitzt. 

Im Laufe der letzten zehn Jahre wurde in der UdSSR eine neue Richtung 
der Theorie entwickelt und zu beachtlichen Erfolgen gefiihrt, die nichts tiber 
die Kleinheit der Wechselwirkung zwischen Elektronen und Atomschwin- 
gungen annimmt. Diese Wechselwirkung, oder manchmal ihr Hauptanteil, 
wird in den HAMILTONoperator des Systems schon in der nullten Naiherung 
einbezogen. Hierbei ergeben sich nicht nur quantitativ, sondern auch quali- 
tativ neue Ergebnisse. Die Form der Absorptions- und Lumineszenzbanden 
der verschiedensten lokalen Zentren im Kristal] findet eine zwanglose Er- 
klarung. Halbwertsbreite und Temperaturabhingigkeit dieser Banden, die 
Stokessche Verschiebung, das Verhiiltnis zwischen der Energie oe des 
optischen Ubergangs und der Energie W des entsprechenden thermischen 
Ubergangs lassen sich berechnen. In vielen Fiillen steht die Theorie zum 
erstenmal in quantitativer Ubereinstimmung mit dem Experiment. Quali- 
tativ neue Ergebnisse erhielt man bei der Behandlung der Quantenzustiinde 
der Leitungselektronen. 

Wir geben nun einen kurzen Abrif dieser neuen Richtung der Theorie. 
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2. Die stationaéren Zustande der lokalisierten 
schwach gebundenen Elektronen 


[PEKAR (1952b)] 


Es handelt sich um Elektronen, die in der Nahe einer Stelle des Kristall- 
gitters lokalisiert sind, an der die Periodizitit gestért ist. Man nimmt an, dafs 
diese Elektronen bedeutend schwicher gebunden sind und sich weitaus lang- 
samer bewegen als die Higenelektronen des Isolators. Hier kommen z. B. die 
an Gitterfehlstellen lokalisierten Leitungselektronen, die Valenzelektronen 
von Fremdatomen, die Eigenelektronen des Isolators, welche durch Anregung 
in Zustiinde mit groBen Wirkungsradien iibergegangen sind und infolgedessen 
schwach gebunden sind, in Frage. Ein Kriterium fiir die schwache Bindung 
eines Elektrons ist die Bedingung, daf die Frequenz des intensiv von ihm 
absorbierten Lichtes wesentlich kleiner ist als die Frequenz des von den Eigen- 
elektronen des Isolators absorbierten Lichtes. 

In diesem Fall kann man die bekannte adiabatische Naherung benutzen 
[PAULI (1933)], bei der man annimmt, daB der Zustand der Eigenelektronen 
der Bewegung der Atomkerne und der schwach gebundenen Elektronen 
adiabatisch folgen kann. Hierbei braucht man nur die konservative Bewegung 
der langsamen Teilchen, der Atomkerne und der schwach gebundenen Elek- 
tronen, zu betrachten; das Vorhandensein der Eigenelektronen des Kristalls 
kommt lediglich in der potentiellen Energie der langsamen Teilehen zum 
Ausdruck. 

Die HAmiILTonfunktion des Systems der langsamen Teilchen hat die Form 


mn 42 ON 1 ; a2 
H = — et V (r,q) +5 Yiu (a—35)- (6) 
= = x x 
Hierbei stellt das erste Glied den Operator der kinetischen Energie der schwach 
gebundenen Elektronen des lokalen Zentrums, das zweite Glied die potentielle 
Energie der Wechselwirkung dieser Elektronen unter sich und mit den Atomen 
(fonen) des Kristalls dar; r ist die Zusammenfassung der Koordinaten siimt- 
licher schwach gebundenen Elektronen des Zentrums, q die Zusammen- 
fassung simtlicher Normalkoordinaten q,, die die Atomschwingungen _be- 
schreiben, w, sind die Frequenzen der EKigenschwingungen der Atome. Das 
letzte Glied in (6) ist der Operator der Normalschwingungen der Kristall- 
atome, wobei die schwach gebundenen Elektronen wegzudenken sind, wihrend 
im tibrigen das lokale Zentrum ebenfalls vorhanden ist. 
In der Folge wird nur von den schwach gebundenen Elektronen die Rede 
sein, die wir einfach als Elektronen bezeichnen werden. Der explizite Aus- 
druck fiir V (r, q) ist unbekannt; trotzdem erhilt man konkrete Ergebnisse. 
Unter der Annahme, da die Schwingungsamplituden der Atome klein sind, 
kann man die potentielle Energie in eine Reihe nach Potenzen der qx. ent- 


wickeln : x A eae 
V (”, q) So Vo (7) oft 2 Ve (7) Gases (7) 


Zur Bestimmung der stationiiren Zustiinde und der Energieterme des Systems 
muB man die Kigenwerte und Eigenfunktionen des Operators (6) aufsuchen. 


~~ 
ses 
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Da die Atome mit ihren verhiltnismaiSig groBen Massen sich bedeutend lang- 
samer bewegen als selbst die schwach gebundenen Elektronen, kann man das 
System wiederum in zwei Untersysteme aufteilen: ein langsames, naémlich die 
Atomschwingungen, und ein schnelles, die schwach gebundenen Elektronen. 
Ferner kann man wieder annehmen, da8 der Zustand der Elektronen der 
Atombewegung adiabatisch (trigheitslos) folgt. Fiir jede augenblickliche 
Atomfiguration q, stellt sich sofort der entsprechende stationiire Elektronen- 
zustand ys (7, g) ein. Er ergibt sich aus der Wellengleichung 


h2 N 
se 2 + Vir, a| Ws (1, 7) = Es (q) ws(", Q); (8) 
4M i—1 4 


wobei s die Gesamtheit der Quantenzahlen der Elektronen und q ein viel- 
dimensionaler Parameter ist. Bei der Lésung dieser Gleichung nehmen wir 
an, da die Atome vorwiegend in der Umgebung derjenigen Lagen schwingen, 
in denen ihre potentielle Energie ein Minimum ist; in der Folge wird sich 
diese Annahme bestatigen, nachdem die Wellenfunktion des Systems be- 
stimmt ist. g,¢ seien die Koordinaten dieser Gleichgewichtslagen. Bei den 
erwihnten Atomschwingungen schwankt V (r,q) leicht um die Funktion 
V (r, qs). In nullter Naherung kann man JV (r, g) in (8) durch V (r, qs) ersetzen. 
Man erhalt dann 


h2 N ‘ 
\—sa 24 ee, a Ys (1, 7s) = Hs (qs) Ps (7 Gs) - (9) 


Die Korrektur zum Potential 


V(r, g) — V (9, 9s) = D> Vu (7) (Ge — Ges) + °° (10) 


x 


wird dann als kleine Stérung eingefiihrt, und in erster Niherung ergibt sich 


E,(q) = Es (4s) + > Vie gh@e Gua) 12° *; (11) 
Views 
Ws (7, 7) = Ys (1, Ws) oa Ey (qs) — Es (Qs) Ws (1, Ys) (Ge — Yas) (12) 
mit ; 
Vins= | Ae (r) | ws (7, qs)? ar, (13) 
Vara = [Vu (1) welts 4s) Yo (Ps %) ar. (14) 


Wir gehen nun zur Behandlung der Atomschwingungen tiber. GemiB der 
adiabatischen Niherung [PAULI (1933)] hat der Energieoperator der Atome 
die Form 


A 1 _ 
H — E, (q) + -S how, (2 3a) (15) 


bol 


Die Rolle der potentiellen Energie der Atome spielt die Funktion 


| 


F (Gd) = £5(¢) + = D hough. (16) 


3 
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Mit Einfiihrung der Elektronen des ,,Zentrums‘ erscheint also in der poten. 


tiellen Energie der Atome das Zusatzglied (11), d. h. die Gesamtenergie dieser 
Elektronen. 

Bricht man die Entwicklung auf der rechten Seite von (11) nach dem ersten, 
von q, unabhangigen Glied ab, so ergibt sich, da die Kinftthrung der Elek- 
tronen tiberhaupt ohne Kinflu8 auf die Atombewegung ist. Dies entspricht 
eben der Vernachlassigung der Wechselwirkung zwischen Elektronen und 
Atomschwingungen, d.h. der zweiten vereinfachenden Annahme der im 
vorigen Paragraphen erwihnten alten Theorie. 

Bricht man die Entwicklung auf der rechten Seite von (11) erst nach denw 
zweiten, von q, linear abhingigen Glied ab, so beginnen nach Formel (16) 
mit Einfihrung der Elektronen des Zentrums in den Kristall zusatzliche 
Krafte auf die Atome zu wirken, wobei entlang der Koordinaten qg, konstante 
verallgemeinerte Kriifte V,, in Erscheinung treten, die zu einer Verschiebung 
der Gleichgewichtslagen der Atome AnlaB geben, ihre Schwingungsfrequenzen 
aber nicht beeinflussen. Weiter unten beschriinken wir uns auf diese Naiherung, 
beriicksichtigen also in der Entwicklung (11) nur die linearen Glieder. Es 
zeigt sich, daB diese Niherung schon die Erklarung der uns interessierenden 
Effekte liefert: die Warmeerzeugung bei den optischen Ubergingen, die 
STOKESsche Verschiebung, den groBen Unterschied zwischen der Energie 
des optischen Uberganges und der thermischen Aktivierungsenergie des ent- 
sprechenden Uberganges, die groBe Halbwertsbreite der Absorptions- und 
Lumineszenzbanden usw. Beriicksichtigt man auch die quadratischen Glieder 
der Entwicklung, so kann man an den gewonnenen Ergebnissen Korrekturen 
héherer Ordnung anbringen. 

Setzt man (11) in (16) ein, so kann man leicht die Koordinaten der Atome 
berechnen, bei denen die potentielle Energie /’, (q) ein Minimum ist. Man erhalt 


Ogee ees (17) 
Das Minimum der potentiellen Energie betrigt 


7 i %, 2 
Js a F's (Js) sa a Pa h@y Gas (18) 
Die potentielle Energie der Atome lift sich schreiben 
7 | 7 
F.(q) = Js = 9 >: ho, (Vx sneha 5 (19) 


sod 


und ihr Energieoperator (15) liBt sich foleendermaBen umformen: 


. mais 0 
Jake —— Af ain 9 = hwy a. — Gas)? ial (20) 
Die Higenwerte dieses Operators sind 
‘ ; ] 
ieee: no =JIs+ D hor, (1.4 5): (21) 
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wobei die », ganze Zahlen sind; seine Eigenfunktionen lauten 
eS Myre (q) =] ] Da (Qu — Tus) s 


1 
oe (%% — % s)* 


Pn, (Ie — Yes) = Ane ‘Hn, (Ge — Ges) + 


(22) 


Hierbei sind die H,, die HERMITEschen Polynome vom Grade n,: Ay, sind 
Normierungskonstante. 

Aus der Theorie der adiabatischen Naherung geht hervor, daB die Eigenwerte 
des Operators (2) oder (15) gleichzeitig die Energieterme des Gesamtsystems 
bilden. Die Wellenfunktionen des Systems haben die Form 


POE ie ee a, (1,9) Di.. wa (9) (23) 


Neben der oben betrachteten adiabatischen Naherung benutzt man zur 
Bestimmung der Eigenwerte des Operators (6) auch die direkte Variations- 
methode. Bekanntlich sind die Eigenwerte des Operators (6) identisch mit 
den Extremwerten des Funktionals 


HY (r,9)) =| P* HW dr dq, 

st (24) 

dr = II dz, dy, dz;, dq =I dq, 
t=1 % 

unter der Zusatzbedingung der Normierung von Y. Die Funktionen Y (r, q), 
die das Funktional H [¥Y) zum Extremum machen, sind die Eigenfunktionen 
des Operators (6). 
Zur Bestimmung der angeniherten Extremwerte von H [¥| approximieren 
wir Y durch den Produktansatz 


PG, qd) = YTV OD, t. woe  (Q). (25) 


Bilden wir das Extremum des Funktionals H [ys:®s...n,...] beztiglich 
®,...n,... bei festem beliebigen ys, so erhalten wir 


es, Vues [We] BY a 
Deng (O) = TE Png (tu + = ) (26) 


“% 


wobei ®, durch Formel (22) und Vz. [ys] durch die Beziehung (13) gegeben 
werden, Setzen wir (26) in H [ps:®;...n,...] ein, so erhalten wir ein Funktional, 
das nur noch yon ws abhangt: 


] 
H [y.(r)] = J [ye(r)] + Y ho, (», + =) (27) 
Sassi, free Wk Liev hie 
J [Ws (r)| a Ws ‘ak Im As oh Vi | sd 2 pe a (28) 


Dieses Funktional mu durch Variation von y, zum Extremum gemacht 
werden. Dabei ergibt sich fiir y, folgende EULERsche Gleichung: 


aD) 


2 N V 
i 2 Ai 2 Vo ir) a. % ie Vu (| Ws = Eys. (29) 
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Die dargelegte direkte Variationsmethode liefert die giinstigste Losung inner- 
halb der Klasse der durch den Produktansatz (25) darstellbaren Funktionen, — 
Die Gleichung (29) ist identisch mit Gleichung (9); hiervon iiberzeugt man 
sich leicht, wenn man Formel (7) beriicksichtigt und ferner beachtet, dab 
Gleichung (9) nicht linear ist, sondern ein ,,self consistent*’ Problem dar- 
stellt: die in dieser Gleichung vorkommenden Parameter q,s, sind, wie aus 
den Formeln (17) und (13) hervorgeht, selbst Funktionale von ys. Aus der 
Identitit der Gleichungen (29) und (9) folgt, daB die in der direkten Variations- 
methode auftretende Funktion y,(7) mit der Funktion y, (7, gs) tiberein- 
stimmt, die in der vorher erérterten adiabatischen Naherung vorkommt. 
Ebenfalls ergibt sich die Ubereinstimmung der Groen V,, und der Funk- 
tionen ®,....... (g), die in beiden Methoden eine Rolle spielen. Durch Ein. | 


#% 


setzen von ws(7) = Ws(7, 7s) in das Funktional (28) erhalt man ferner | 


1 : od | 
J [Ys] = Be (Qs) + hOn Gis =I. (30) 


Die Energie des Systems wird durch die Extremwerte des Funktionals (27) 
gegeben; das stimmt iiberein mit dem Ausdruck (21) der adiabatischen | 
Naherung. 

Die Ergebnisse der direkten Variationsmethode sind also identisch mit denen | 
der adiabatischen Niherung, falls man sich in dieser bei Anwendung der oben | 
dargelegten Storungsrechnung auf die erste Naherung in der Energie und die | 
nullte Niherung in der Wellenfunktion des Systems beschrankt. 

Eine exakte Lésung der nichtlinearen Integro-Differentialgleichung (29) ist | 
mathematisch auBerst schwierig. Ihre Niherungsl6sung und eine explizite : 
Bestimmung von ws(r) wird unten fiir mehrere Spezialfalle durchgefiihrt, 
fiir die das Modell des lokalen Zentrums wohlbekannt und einfach ist. In 
anderen Fallen kann man ohne Benutzung eines Modells des Zentrums und 
eines expliziten Ausdrucks fiir die Funktion ws(”) eine Reihe interessanter! 
Beziehungen ableiten. 


3. Optische Elektroneniiberginge 


[PEKAR (1952b)] 


In diesem Paragraphen werden die optischen Ubergiinge der schwach ge-» 
bundenen Elektronen eines lokalen Zentrums von einem diskreten auf einenr 
diskreten Knergieterm behandelt. Wir betrachten einen Ubergang vom: 
Zustand s;...m,... in einen Zustand s,...n,... Bei diesem Ubergangy 
andern sich nicht nur die Quantenzahlen s der Elektronen, sondern auch die 
Schwingungsquantenzahlen »,, d.h. es wird Wiirme erzeugt. Die Frequenz: 
des bei diesem Ubergang absorbierten Lichts betriigt 


<b In yp) , , 
(i "i (eles ata hy EK, ee tyes | = oe S- DY Oy (Ne — Nx) , 
x 
] 
h 
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wobei w° die bei einem ,,rein elektronischen‘‘ Ubergang absorbierte Frequenz 
ist; dieser rein elektronische Ubergang ist charakterisiert durch die Be- 
dingung n;, = nx. 
Die Wahrscheinlichkeit fiir den optischen Ubergang liBt sich auf Grund der 
tiblichen quantenmechanischen Betrachtung berechnen [PEKAR (1952b)]. Die 
bei dem erwihnten Elektronenschwingungsiibergang in der Sekunde und. der 
Volumeneinheit des Isolators absorbierte Lichtenergie kann man folgender- 
mafBen schreiben: 

tS (w), (32) 


hierbei ist S (@) die spektrale Dichte der Intensitit des einfallenden Lichtes, 
und der Koeffizient + betrigt 
Anew N (s,...Nx 


ar 
) 7\2 Hs 
: 3Acn(@) [Dtesn, a.07 (2 (33) 


Dabei ist c die Lichtgeschwindigkeit im Vakuum, MW (s,...n,...) die Kon- 
zentration der Zentren, die sich im Ausgangszustand s,...n,... befinden; 
n (q@) ist der Brechungsindex des Lichts im Isolator: 


WMts.n, CN ge Wee, : M,, re (34) 
lao —s po; ere tre te. a eo th) Oly at (35) 
Man’ = 11 | Dn, (Te — Uxs,) Dui (Te — Yxs,) Ux (36) 


dabei sind r; die Koordinaten der optischen Elektronen des Zentrums, wiihrend 
der vektorielle Koeffizient f(r;) so definiert ist, daB das Produkt 


e (f (r.), €) (37) 


die Stérenergie des i-ten optischen Elektrons durch die Lichtwelle © und das 
auBere Feld der Polarisation der stark gebundenen Elektronen, das ebenfalls 
proportional € ist, darstellt. In Formel (33) ist tiber alle méglichen Richtungen 
des auf das Zentrum fallenden Lichtes gemittelt. Diese Formel unterscheidet 
sich von der bekannten Formel fiir die Absorption von Licht durch Atome 
im Vakuum durch den komplizierteren Ausdruck fiir das Matrixelement 
Mts,n, sn’, Sowie auch durch den Faktor n (@) im Nenner. 
Unter Benutzung der bekannten EINSTEINschen Uberlegungen iiber den 
Zusammenhang zwischen Absorptions- und Emissionswahrscheinlichkeit und 
in Verallgemeinerung dieser Betrachtungen auf ein Medium mit einer optischen 
Dichte »(q@) erhalten wir fiir die Energie, die in der Sekunde und in der 
Volumeneinheit beim Ubergang s8,...mj,...—> 5,...n,... ausgesandt wird, 
den Ausdruck 

4@*n(w) WN (83... - 


Y 
S — 
j 363 


2) | Mean, sen” |? (38) 


/ 


in diesem bedeutet WN (s,...n/,...) die Konzentration der Zentren im Zustand 


| Sy...m,... Die in §2 benutzte Niherungsmethode zur Behandlung des Zu- 


standes des Systems, die zu einem Produktansatz (25) fiir die Wellenfunktion 


te oem 
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des Systems und zu dem additiven Ausdruck (21) fiir die Energie fiihrt, gibt 
die Méglichkeit, die Zustinde der Elektronen, gekennzeichnet durch die 
Quantenzahl s, und die Zustiinde der harmonischen Schwingungen, deren 
jede durch eine besondere Quantenzahl n, gekennzeichnet wird, in gewissem 
Sinne getrennt zu behandeln. Hierbei braucht nicht vorausgesetzt zu werden, 
daB die Wechselwirkung zwischen Elektronen und Schwingungen klein ist, 
Weiter unten wird stets nur von zwei Elektronenzustinden die Rede sein: 
von dem Grundzustand s, und einem angeregten Zustand s,. In jedem dieser 
Zustiinde kénnen die n, beliebige Werte annehmen. ; 
Wir betrachten den Ubergang s,...n,...—> 82...N,-.., bei dem die Frequenz 


w +a, Wo w= > @, (m,—N,), (38a) 


absorbiert wird, und den Ubergang $5...nz...—> 5 ...N,-.., bei dem die 
Frequenz 0 6) (38b) 
emittiert wird. 
Die Frequenzen w° + @, und w?— a, liegen symmetrisch beiderseits der 
Frequenz w® des rein elektronischen Ubergangs. Auf Grund der Formeln (34), 
(35), (36) und (22) kann man zeigen [PEKAR (1952b)], dab 

| Ween, syn’ ae Msn, syn M4 (39) 


ist. Die Wahrscheinlichkeit der oben erwihnten absorbierenden und emittieren- 

den Ubergiinge sind also verkniipft. Nimmt man an, daB die Schwingungen 

vor jedem optischen Ubergang im thermischen Gleichgewicht stehen, so. 
ergibt sich [PEKAR (1952b)]: 


om ate NiP_S, (40) 
WM |w 4 (Oy) h NN; a NM | wo — or 


Hierbei ist V, die Konzentration der Zentren, in denen sich die Elektronen 
im Zustand s, befinden, unabhingig vom Schwingungszustand; JN, ist die: 
Konzentration der Zentren, deren Elektronen sich im Zustand s, befinden. 

Die Formel (40) zeigt, dai die Kurven fiir tn/@ und S,/o*n als Funktionen | 
der Frequenz symmetrisch zur Linie @ = @® gelegen sind. Die zu @, sym-: 
metrische Lage des Absorptions- und des Lumineszenzspektrums wurde: 
experimentell zuerst von LJOWSCHIN im Jahre 1931 beobachtet [LJOWSCHIN | 
(1951)], welcher in seinen darauffolgenden Arbeiten zeigen konnte, dal: 
diese Symmetrie eine Kigenschaft der Spektren sehr vieler Substanzen (Farb- 
stoffe) in Lésungen ist. Kine solehe Symmetrie beobachtet man auch in den 
Spektren vielatomiger Molekiile. Die oben angegebene theoretische Ableitung: 
bezieht sich zwar auf den Fall eines Kristalls, die regulire und periodische: 
Anordnung der Atome wird jedoch nirgends vorausgesetzt. Es wird nur an-. 
genommen, dafs die Atome harmonische Schwingungen um irgendwelche: 
konstanten Gleichgewichtslagen ausfiithren. Soleche Schwingungen finden aber’ 
auch in Molekiilen und sogar in Fliissigkeiten innerhalb ziemlich groBer Zeit- 
riume statt, bis infolge einer Fluktuation der Schwingungen eine Durch-, 
mischung der Fliissigkeit eintritt und sich neue Gleichgewichtslagen ein- 
stellen. Gewohnlich ist die Dauer einer solchen Durchmischung linger als: 
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die des optischen Uberganges. Die oben gewonnenen theoretischen Ergebnisse 
(sowie auch ein Teil der untenstehenden) gelten also nicht nur fiir den festen 
Korper, sondern auch fiir vielatomige Molekiile und z. T. auch fiir Fliissig- 
keiten. 

u beachten ist, daB die Formeln (39) und (40) Naherungscharakter haben 
d sich nicht hitten ableiten lassen, wenn wir uns in der Entwicklung (11) 
icht auf das lineare Glied beschriinkt, sondern auch das quadratische beriick- 
ichtigt hitten. Dies ist gleichbedeutend mit einer Beriicksichtigung der Ab- 
angigkeit der Higenfrequenzen ow, der Schwingungen vom Elektronen- 
ustand 8 


4. Form und Temperaturabhingigkeit der Absorptions- 
und Lumineszenzbanden 


ir betrachten eine Absorptionsbande, die aus einer Linienserie entsteht, 
eren Linie einen bestimmten Elektroneniibergang s,—> s,, aber den ver- 
hiedenen médglichen Schwingungsiibergingen ...n,—->...m,... ent- 
spricht. Die Frequenzen dieser Serie werden durch Formel (31) gegeben, 
enn man in dieser die Quantenzahlen s, und s, festhilt, den Zahlen n, und nj, 
aber alle méglichen ganzen positiven Werte beilegt. Die Intensitat der Linien 
aBt sich nach (33) bis (36) folgendermaBen schreiben: 


t= Gi.3, ba (41) 
abei ist: 
TEIN Spe «Myce ) 
G,.. = oe eke: 42 
tie 3Bhon (a) | Mse| Cs 


Verindert sich das Verhaltnis w/n(m) innerhalb der Bande verhaltnismaBig 
wenig, so kann man G,.., als konstanten Koeffizienten ansehen. Form und 
Halbwertsbreite der Bande werden also durch den Faktor |M nn’ |? bestimmt. 
Aus der Formel (31) geht hervor, daB eine gegebene Frequenz m von unend- 
lich vielen elementaren Quanteniibergiingen absorbiert wird. 

Es handelt sich um alle die Ubergiinge, fiir die die Summe 


D On (14, = ng) (43) 


x 


‘einen gegebenen Wert hat, den wir in der Folge mit a, bezeichnen werden. 
Um die gesamte absorbierte Intensitait der Frequenz 


o= w+ a, (44) 


‘zu erhalten, hat man die GréBe (41) iiber siimtliche Elementariibergiinge zu 
summieren, bei denen die gegebene Frequenz absorbiert wird. Diese Sum- 
mierung fiihrt auf schwierige Rechnungen. Der Kiirze halber geben wir nur 
die Endergebnisse an. 

\Die exakte Durchfiihrung der Summation gelang in dem Grenzfall geringer 
Dispersion der Normalschwingungen der Atome (Ionen) im Kristall, d. h. 
in dem Fall, da8B simtliche Higenfrequenzen wm, fast gar nicht von x abhingen 


ij 


] 
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und gleich der Grenzfrequenz Wo sind ([PEKAR (1950a)}. In diesem Fall kann | 
man die Formel (31) folgendermafen umformen : ; 
i” wo = w + PW, «| | 
dabei ist pease a 


eine ganze Zahl, die sowohl positive als auch negative Werte annehmen kann, 
Das Absorptionsspektrum besteht also aus einer Reihe diskreter Linien, die 
voneinander den gleichen Abstand wp) haben. Die gesamte in der Sekunde 
und der Volumeneinheit des Kristalls von der Linie w° -+ pm, absorbierte 
Lichtenergie betragt 
Tp S(w° + pap) 


dabei ist!) 


=e (7 ny 7) a é = pT ae é | 
Ep ola e = (1 + =) Ty(aV Np (M + 1)), (47) 
0 
ity ist der mittlere PLANCKsche Wert fiir die Oszillatorenquantenzahl n,: 
ja 
1 
Nie ae PS Ree (48) 
ekT__] 
I, (z) ist die BEssELfunktion reellen Arguments: 
1 : 
Tp (2) = Io (i?) (49) | 
a eine Konstante vom Wert 
a= 3! (Gus, — Ys,)?- (50) | 


Wenn eine Dispersion der Normalschwingungen besteht, w, aber nahezu 
gleich @, ist, so erhilt man an Stelle der oben erwihnten diskreten Linien : 
schmale Banden von endlicher Breite, und die Formel (47) gibt die integrale ° 
Absorptionsintensitit jeder dieser Banden an. Mit weiterem Ansteigen der | 
Dispersion wiichst die Breite dieser Banden, sie beginnen sich zu iiberdecken, . 
und das Absorptionsspektrum wird kontinuierlich, hat aber immer noch eine : 
Struktur: es bildet eine oszillierende Kurve mit abwechselnden Maxima und | 
Minima, wobei der Abstand zwischen den Maxima etwa wp, betrigt. Diese: 
Kurve kann man als breite Glockenkurve darstellen, iiber die eine Sinuskurve 
mit der Periode @, tiberlagert wird. Die exakte Gleichung dieser Kurve findet t 
der Leser bei PEKAR und KRIWOGLAS (1953). Mit weiter ansteigender Dis-- 
persion geht die Amplitude der Sinuskurve gegen Null, und das Spektrum: 
wird zu einer glatten Glockenkurve. Deren Form erhalt man, wenn man’ 
beriicksichtigt, daB in dem Spektralbereich A @ gerade A m/m) Banden liegen, 
t) Anm. d. dtsch. Red.: Unabhingig von PEKAR haben HuancG und Ruys (Proc. Roy.’ 
Soc. (A) 204, 406, 1951) eine Theorie angegeben, die zu einer Formel gleicher Bauart} 
fiihrt. Vgl. auch O’RourKE, Phys. Rev. 91, 265, 1953. } 

1 

| 
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von denen jede in der Sekunde eine Energie t)S absorbiert. In dem Spektral- 


TpS(@)A@ 


bereich A @ wird also die Energie absorbiert. Diese Energie schreibt 


0 
man am besten in der Form t(@)S(@)4@, wobei t(@) der gewéhnliche Ab- 
sorptionskoeffizient fiir ein kontinuierliches Spektrum ist. Es ist also 


t(w) =—. (51) 


In dieser Formel ist p nicht mehr als ganzzahliger Index aufzufassen, sondern 
als kontinuierlich verinderlicher Parameter, der mit der Frequenz @ vermége 


Die Formeln (47) und (51) bestimmen die Form des Absorptionsspektrums 
nd seine Temperaturabhingigkeit, da ja m) von der Temperatur abhiingt. 
ei der Herleitung dieser Formeln wurde vorausgesetzt, daB nur ein Dis- 
persionszweig mit der Grenzfrequenz @) vorhanden ist, und daB w, von wp 
mur schwach abweicht. Die oben erwihnte Uberdeckung der Banden tritt 
haufig schon ein, wenn die Dispersion der w, 10% von wy betragt. 

Auch der allgemeinere Fall eines Kristalls mit beliebiger Struktur der Ele- 
mentarzelle, beliebiger Anzah] von Dispersionszweigen und beliebig hoher 
Dispersion, die zu einer Uberdeckung der Banden fiihrt, 1a8t sich behandeln 
[PEKAR und KRIWoOGLAS (1953)}. Im Gegensatz zu der vorstehenden Ab- 
leitung trigt diese Behandlung jedoch Naherungscharakter und gilt nur fiir 
wei Grenzfalle exakt, namlich fiir den Fall groBer und geringer Warme- 
erzeugung beim optischen Ubergang, genauer dann, wenn die mittlere Anzahl 
der beim Ubergang erzeugten Schwingungsquanten ho, bedeutend gréBer 
bzw. bedeutend kleiner als Eins ist. 

Im Fall groBer Warmeerzeugung wird die Form des Absorptionsspektrums 
durch die Formel 


a eg arn ae TCT re 
T= Ty exp | og = li 6 J’? a it (52) 


gegeben: hierbei ist mm, die Frequenz, bei der der Absorptionskoeffizient t 
maximal ist; sie betragt 


: 
Ong = OL D> (Gs, —Gne,)? On, * (53) 
wu 4 al hes 1 5 
(Gj) T= Pe (xs, — Tus. e Wy | Mx + > |? (54) 
vr 1 27,8 55 
Rice > (us, — Tes.) Oe; 8) 
1 RK 
ey. 56 
Ny hota > (5 ) 
Pet es | 
pee ata (57) 
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Wie aus Formel (52) hervorgeht, konvergiert im Gebiet des Absorptions- 
maximums wegen der Kleinheit von @ — @m die Potenzreihe im Exponenten 
der Exponentialfunktion schnell, so da man sich auf das erste, quadratische 
Glied beschrinken kann. Mit wachsender Entfernung vom Maximum beginnt_ 
allmaihlich auch das kubische Glied eine gewisse Rolle zu spielen, so da8 di 
Absorption nach der roten Seite hin etwas schneller abnimmt als nach der | 
violetten. Eine niihere Untersuchung zeigt, daB im Fall hoher Warme- 
erzeugung das kubische Glied weitaus kleiner als Kins ist und bis zu deh | 
Punkten fortgelassen werden kann, in denen t(@) gleich der Hilfte seines ; 
Maximalwertes ist. Innerhalb dieses Gebietes wird das Absorptionsspektrum | 
durch eine Gausskurve mit der Halbwertsbreite 
pea | 
dm =2Y2In2 Vg". (58) ) 
dargestellt. Dieses Ergebnis sowie auch die Formeln (53), (54) und (57) wurden 1 
vom Verfasser zunichst fiir den Fall geringer Dispersion abgeleitet [PEKAR} 
(1950a)] und dann auf den Fall beliebiger Dispersion verallgemeinert [PEKAR } 
und KRIWOGLAS (1953)]. 
Aus den Formeln (57) und (58) geht hervor, dai das Produkt 


In2 
Tm ° 30 = Gos, 2 | 5 (59) 


konstant ist. Dies bedeutet, daB die Flache unter der Glockenkurve t(q) } 
temperaturunabhangig ist. 

In zwei Grenzfillen kann man aus den Formeln (53) und (54) die unbequemen, . 
schwer zu berechnenden unendlichen Summen eliminieren: 


A. Fall ohne Dispersion bei beliebiger Temperatur. In diesem Fall kann man 
w, durch wp, ersetzen und wm, und %, vor das Summenzeichen ziehen. Dann~ 
ergibt sich aus (53) und (54) 


g” = 205 (@m — w°) (7 aia 3). (60) | 


dw = 2Y2In2 2 Wo (Om — @®) (*. f 5): (61) 


B. Fall beliebiger Dispersion und hoher Temperatur, bei der 

: 1 

Ny, = ES ee a (62) 
ist. In diesem Fall folgt aus (53) und (54) 


2kT 
| ho — a (Wm — @p), (63); 


as I 
0M = 2y2n2 VF (Wm — @)) kT. (64) 
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Samtliche Darlegungen dieses Paragraphen bezogen sich auf eine Absorptions- 
bande. Das Symmetriegesetz (40) gibt die Méglichkeit, simtliche oben ge- 
wonnenen Ergebnisse auch fiir die Lumineszenzbande hinzuzuschreiben, die auf 
dem umgekehrten Elektroneniibergang s,—> s, beruht. Nach Formel (40) werden 
die GréBen tr/w und S,/w*n als Funktionen von @ durch spiegelbildliche 
Kurven dargestellt. Symmetriezentrum ist der Punkt w°®. Bei nicht allzu 
groBer Breite der Banden stimmen der Verlauf dieser GréBen, ihrer Halb- 
wertsbreite, die Lage der Maxima usw. ungefiihr mit denen fiir t(@) und 
S,(@) iiberein. Oben wurde gezeigt, daB die GréBe tn/w bei der Frequenz am, 
die durch Formel (53) definiert wird, ein Maximum hat. Infolgedessen hat 
die GréBe S,/e@*n ein Maximum an dem spiegelbildlichen Punkt 


1 
Om = 0 — > J (Gus, — Guay)” On: (65) 


Nimmt die erste dieser GréBen zur roten Seite hin schneller ab als zur violetten 
(infolge des kubischen Gliedes in (52)), so nimmt die zweite umgekehrt zur 
violetten Seite hin schneller ab. Die Kurven fiir die GréBen tn/w und S,/a*n 
besitzen gleiche Halbwertsbreiten. Die Halbwertsbreiten der Absorptions- 
und der Lumineszenzbande sind also ungefihr gleich. 

Die StToKeEssche Verschiebung, d.h. der Abstand zwischen dem Maximum 
yon t(@) und dem von S,(q@) betragt ungefahr 


A@ = On — On = 2(@m — w°) = > (Gus, — Uns.) Ox- (66) 

Im Fall A) haben wir also 

—— i) 
bo =2V2In2 |) Ao oy (my +5) (67) 
im Fall B) 
——~ 4/Aw 
Ow => 2yama y 42 kT. (68) 
U 


Die Formeln (67) und (68) sind dadurch bequem, daf in ihnen keine un- 
bekannten Parameter vorkommen, so daB man sie leicht mit dem Experiment 
vergleichen kann. A 
Bei dem absorbierenden und dem emittierenden optischen Ubergang, die 
zueinander spiegelbildlich sind, tritt dieselbe Warmeerzeugung, nimlich 
D> h@y (Ne — Nx). (69) 

auf. Im Gebiet des Absorptionsmaximums betraigt diese Wiairmeerzeugung 
‘nach den Formeln (31) und (53) 
hAw Uh gens 

Ge ety Shy (Ges, — Us,)” (70) 


a “ 


h(@m — w°) = 


. Dieses Ergebnis laBt eine einfache klassische Deutung zu: bedenkt man, da 
mit unseren Bezeichnungen die potentielle Energie der Atomschwingungen 
die Form 1 


— Sho, (Gx — Ins)? , (71) 


2 
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> 


hat, so ist die rechte Seite von (70) die Energie der Schwingungen, die die 
Atome infolge der momentanen Anderung der Krifte und der Verschiebung 
der Gleichgewichtslagen q,,, und q,., ausfiihren, j 
Unter der Annahme, da sich der Kristall im thermischen Gleichgewicht 
befindet, und unter Benutzung der iiblichen kanonischen Verteilung und der 
Formel (21) kann man die Wahrscheinlichkeit w (s) dafiir berechnen, daB sich 
die Elektronen des Zentrums im Zustand s befinden, wihrend die n, beliebige 
Werte haben kénnen. Diese Rechnung liefert 


So (72) 


w (8;) 


A Aa 7 | 
Wis = Js, ~Ja= be = b (om —") = 4 (ale +S). (73) | 


» 
\ a 


W,. bezeichnet man als thermische Anregungsenergie der Elektronen fiir den | 
Ubergang aus dem Grundzustand s, in den Zustand s,. Formel (73) zeigt, , 
daB W,, nicht mit den Energien der entsprechenden optischen Ubergange » 
h@®m und hw, iibereinstimmt. 

Wir miissen betonen: die oben behandelten Effekte, nimlich die Warme- . 
erzeugung bei den optischen Ubergiingen, die STOKESsche Verschiebung, , 
der Unterschied zwischen der thermischen Anregungsenergie der Elektronen | 
und der Energie des entsprechenden optischen Uberganges, die starke Ver-- 
breiterung der Absorptions- und Lumineszenzbanden, die den Ubergang von 1 
Elektronen von einem diskreten auf einen diskreten Energieterm entsprechen, 
die Temperaturabhingigkeit der Banden — alle diese Effekte lassen sich | 
auf Grund der Tatsache, daB qs, == dus, ist, verstehen. Diese Effekte beruhen | 
auf der Anderung der Gleichgewichtslagen der Atome infolge des optischen » 
Uberganges der Elektronen; sie sind mit einer Deformation des Gitters dureh ' 
die optischen Elektronen verkniipft. 

AuBerst wichtig ist, daB die oben angegebenen SchluBfolgerungen ohne; 
konkrete Annahmen tiber das Modell des absorbierenden Zentrums, tiber den: 
Typ des Kristalls (oder eines amorphen Isolators) und ohne Benutzung eines: 
expliziten Ausdrucks fiir den Elektronenanteil ys,(r) der Wellenfunktion ab- 
geleitet wurden. Diese Ergebnisse besitzen also groBe Allgemeinheit und: 
deuten auf die Existenz einer Verwandtschaft in den Absorptionsspektren: 
der verschiedenartigen lokalen Zentren in den verschiedenen Kristallen hin. 
So muB z. B. im Fall groBer Warmeerzeugung beim optischen Ubergang, 
im Fall nicht allzu kleiner Dispersion und bei Erfiillung der iibrigen Kriterient 
fiir die Anwendbarkeit unserer Theorie (Adiabasiebedingung, Harmonizitit 
der Normaischwingungen der Atome, die Méglichkeit, die quadratischen 
Glieder in (11) zu vernachliissigen) die Form der Absorptions- und Lumineszenz~ 
banden stets einer GAUsSkurve iihneln, die durch einen einzigen Parameter g’” 
oder eine Halbwertsbreite (58) gekennzeichnet wird. Die letztere hingt 
mit der StoxkEsschen Verschiebung vermége der Formeln (67) und (68° 
zusammen. 

Oben behandelten wir den Fall groBer Wiirmeerzeugung beim optischen Uber! 
gang. Auch der entgegengesetzte Grenzfall kleiner Wirmeerzeugung, d. h| 
der Fall, in dem die mittlere Anzahl der beim optischen Ubergang erzeugter) 
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Schwingungsquanten A@ bedeutend kleiner als Eins ist, ist der Behandlung 
zugiinglich. In diesem Fall zeigt es sich, daB die Form des Absorptionsspek- 
trums wesentlich von den Einzelheiten des Spektrums der Eigenfrequenzen «,, 
abhiingt. Es besteht deshalb keine Ahnlichkeit zwischen den Absorptions- 
spektren im Fall verschiedener Kristalle und verschiedener Typen von Zen- 
tren; t(@) lie sich nicht durch so einfache und allgemeine Formeln aus- 
driicken, wie oben im Fall hoher Wiirmeerzeugung. Trotzdem lieBen sich 
allgemeine Formeln fiir t(@) gewinnen. Um diesen Artikel nicht zu iiberladen, 
fiihren wir die Formeln nicht an, sondern nur den einfacheren Ausdruck fiir 
die integrale Intensitit der Bande, d.h. fiir 


T = = [t(o) w)do 


wobei das Integral iiber die ganze Bande erstreckt wird. 
Der Index y numeriere diejenigen Dispersionszweige der Normalschwingungen, 
in denen die Anzahl] der Schwingungsquanten infolge des optischen Ubergangs 
um eine Zahl qg, > 0 wichst: 

Gy De oe aa, Tile) . 


x 


» bezeichnet die Summierung innerhalb eines Dispersionszweiges. Der 


ia . 
Index / numeriere analog diejenigen Zweige der Normalschwingungen der 
Atome, in denen die Anzahl der Schwingungsquanten infolge des optischen 
Uberganges um 7, > 0 abnimmt: 

1, = SI (Mua — Mra) - 
Wir betrachten eine Absorptionsbande, die durch optische Ubergiinge mit 
vorgegebenen Zahlen g, und r; erzeugt wird. Diese Bande erstreckt sich von 
der Frequenz 


. 2 
Onin = O° + ZS Goode — yr wi? (73a) 
2 
bis zu der Frequenz 
(2 1 ; 
Omax = 0° +2 ay wy = th Wa > (73) 
2 


wobei ws, Ou und aw) F ool ’ die kleinste und die erdBte Frequenz in den ent- 


sprechenden Dispersionszweigen sind. Die integrale Intensitit dieser Bande 


betragt 


dp fr 
1 Oy G2 4 me 
TY ide ioe pelts a= Gy Pl ome Tos (73¢) 
qd: 2 LAL 
mit 
On = > (Quxs, — Yuxs,)* (Mx 1); 
on = 2 (Yuxs, — uxs,)” Nyx; (73d) 
Gg Gut On 


y g! 
a G exp a D) ? 
u a J 
der Index yw numeriert siimtliche Dispersionszweige ohne Ausnahme. 
26* 


386 S. I. PEKAR 


In dem Spezialfall, daB nur ein einziger Zweig der Normalschwingungen 
u = 1 mit den optischen Elektronen mit merklicher Intensitat in Wechsel- 
wirkung steht, erhilt man auf der violetten Seite der Linie des rein elektro- 
nischen Uberganges eine Folge nahezu aquidistanter Banden mit den inte- 
egralen Intensitaiten ' 

on 


a (73e) 
2 he 


DG OO 5 x OO we) 2 Gs 


und an der roten Seite eine entsprechende Folge von Banden mit den inte- 
gralen Intensititen 
5 oi ry | 


(73f) | 


T (00. <5: 00 va) Gy 
x 


Die ganzen Zahlen q, und 7, stellen die Ordnungszahlen der Banden in diesen 
Folgen dar, und die Formeln (73e) und (73f) geben die Abhangigkeit der 
Intensitat der Bande von ihrer Ordnungszahl sowie auch von ihrer Tem- | 
peratur an. 

Bei hohen Temperaturen ist of ~ o,, und infolgedessen sind die integralen 
Intensitaiten der Absorptionsbanden an der violetten und an der roten Seite 
der Linie w® des rein elektronischen Uberganges ungefiahr spiegelbildlich 
zueinander. 

Bei tiefen Temperaturen geht die Intensitit der Banden auf der roten Seite 
exponentiell gegen Null, und zwar mit zunehmender Abkiihlung des Kristalls 
immer stiirker (o, > 0); auf der violetten Seite dagegen geht die Intensitit 
der Banden gegen einen endlichen Grenzwert. 
Die Formeln (73¢), (73e) und (73f) wurden in einer Arbeit des Verfassers 
gemeinsam mit M. A. KRIWOGLAS abgeleitet [PEKAR und KRIWOGLAS 
(1953)]. 

Es verdient hervorgehoben.zu werden, daB man im Falle geringer Warme- 
erzeugung aus dem Absorptionsspektrum bedeutend weitgehendere Schliisse : 
iiber das Spektrum der @, ziehen kann als im Fall hoher Wiirmeerzeugung. 


5. Vergleich der berechneten Form und der Temperaturabhangigkeit 
der Banden mit dem Experiment 


Die oben abgeleiteten theoretischen Ergebnisse hinsichtlich der Form und | 
der Symmetrie der Absorptions- und Lumineszenzbanden gelten nur fiir’ 
einfache Banden, d.h. fiir solche, die einem einzigen Elektroneniibergang : 
8,58, entsprechen, wobei alle méglichen Schwingungsiibergiinge n,—> nj! 
vorkommen kénnen, Stellt das experimentell beobachtbare Spektrum eine » 
Uberlagerung mehrerer einfacher Banden dar, so hiingt seine Form wesent- - 
lich von dem Abstand zwischen den einfachen Banden, ihrer Halbwertsbreite ° 
sowie ihrem Intensititsverhaltnis ab. Zum Vergleich eines solehen zusammen- - 
gesetzten Spektrums mit der Theorie hat man es zunichst in einfache Banden ! 
zu zerlegen, was oft nur angeniihert méglich ist. Beim Vergleich der Theorie} 
mit dem Experiment hat man also solchen Fallen den Vorzug zu geben, bei) 


denen einfache Banden oder hédchstens leicht trennbare Doppelbanden zu 
beobachten sind. 
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Die Ahnlichkeit der einfachen Absorptions- und Lumineszenzbanden mit 
einer GAusskurve wurde schon lange vor der Entwicklung der Theorie von 
zahlreichen Experimentatoren bemerkt. Ohne uns allzusehr in die Mannig- 
faltigkeit des experimentellen Materials zu verlieren, geben wir hier nur 
einige typische Beispiele. 

In Abb. 1 ist das Lumineszenzspektrum von CaWO, mit einer PbhWO,- 
Beimischung dargestellt [MoskWIN (1945, 1948); BRINKMANN und VLAM 
(1949)|. Wie aus der Figur ersichtlich, liegen die eal gut auf einer 
Gausskurve. 


199 20 27 
Wellenzah/(cm-’) 


Abb.1: Lumineszenzspektrum von CaWO, mit einer Beimischung von PbWO, bei verschiedenen Tempera- 
turen. Wellenlange des Anregungslichtes 2 = 2537 A. 0: MeBwerte; ausgezogene Linien: GAUSS- 
kurve. 


In Abb. 2 sind die entsprechenden Werte fiir einen MgWQ,-Kristall dar- 
gestellt [MOSKWIN (1945, 1948), BRINKMANN und VLAM (1949), VLAM 
(1949)]. Abb. 3 zeigt das Lumineszenzspektrum eines Zn§-Cu-Kristalls 
{MoskKwIN (1945, 1948), BRINKMANN und. VLAmM (1949)]. Es besteht aus zwei 


_ einfachen Banden, einer griinen, der Hauptbande, die man gewohnlich mit der 


Kupferbeimischung verkniipft, und einer blauen, die auf dem tiberschiissigen 
Zink beruht. 

Unsere Abbildungen zeigen, daf in den entsprechenden vier Fiillen das 
Lumineszenzspektrum gut durch Formel (52) ohne das kubische Glied im 
Nenner dargestellt wird. Dies bedeutet, daB das Kriterium der hohen Wirme- 
erzeugung gut erfiillt ist; weiter unten wird dies nochmals unabhingig be- 
statigt werden. 

Jeder der erwihnten Lumineszenzbanden muf eine Bande der Fremd- 
absorption von gleicher Form und Halbwertsbreite entsprechen. Diese Ab- 
sorptionsbanden fallen jedoch schon ins ultraviolette Gebiet und werden 
durch die Absorption der Grundsubstanz verdeckt. Bis heute gelang es nicht, 
mit hinreichender Genauigkeit die oben erwihnte Fremdabsorption vom 
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gesamten Absorptionsspektrum abzutrennen; es ist also bisher nich 
méglich gewesen, das Spektrum der Fremdabsorption mit der Theorie z 
vergleichen. 
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Abb. 2: Lumineszenzspektrum von MgWQ, bei verschiedenen Temperaturen. Wellenlinge des Anregungs- 
lichts 2 = 2537 A. 0: MeBwerte; ausgezogene Linien: GAUSSkurve. 
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Abb. 3: Lumineszenzspektrum von ZnS+Cu. Wellenlinge des Anregungslichts 2 = 3650 A. o: MeBwerte; 
punktiert: GAUSSkurven; ausgezogene Linien: Summe der GAUSSkurven. 
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Zur experimentellen Bestimmung der Form der Fremdabsorptionsbanden 
mu man diejenigen Fiille heranziehen, in denen diese Bande in das Durch- 
sichtigkeitsgebiet der Grundsubstanz fillt. Als klassisches Beispiel erwihnen 
wir die allgemein bekannten F-Absorptionsbanden der Ionenkristalle. Eine 
typische F-Absorptionsbande ist in Abb. 4 dargestellt [POHL (1937)]. 
Die Form des Spektrums wird durch 

die theoretischen Formeln (47), (51) Wellenlange (mu } 
gut wiedergegeben. Im Grenzfall hoher 300 one tee oe 
Warmeerzeugung geht diese Formel 
in Formel (52) iiber. Die geringe Asym- 
metrie der Kurve (die langsamere Ab- 
nahme des Absorptionskoeffizienten an 
der violetten Seite) stimmt vollauf 
mit der Theorie iiberein. 

Nach der theoretischen Formel (59) 
darf die von der Kurve t(w) um- 
schlossene Flache nicht von der Tempe- 
ratur abhangen. Dies ist tatsichlich 
der Fall und wurde von den Experi- 
mentatoren mehrfach festgestellt. Die 
Unabhiangigkeit der Fliche von der 
Temperatur wird in Abb. 4 deutlich. 
Wir gehen nun zur Behandlung der Tem- 
peraturabhangigkeit der Halbwerts- 
breite der Bande itiber und schreiben 
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Abb. 4: F-Absorptionsbande von KBr bei ver- 


Forme] (67) folgendermaBen um: schiedenen Temperaturen. 
Xa) ] | 
2,354 - VamA o et) Z 
ek? —]} 


Diese Forme] ist anwendbar, wenn die Dispersion der Eigenfrequenzen «, 
klein ist und wenn sie alle durch die mittlere effektive Frequenz wp, ersetzt 
werden kénnen. Genauer gesagt: es ist hinreichend, da die Dispersion in 
dem Teil der Normalschwingungen klein ist, die am stirksten mit den opti- 
sehen Elektronen wechselwirken. Es handelt sich um diejenigen Normal- 
schwingungen, bei denen infolge eines optischen Uberganges die gréBte Ver- 
schiebung qs, — qvs, eintritt. 

In Abb. 5 stellt die ausgezogene Kurve die GréBe (67a) dar. 

Aus dieser Abbildung ist ersichtlich, dai bei den tiefsten Temperaturen da 
und infolgedessen auch die gesamte Form der Bande temperaturunabhingig 
ist; dies steht im Einklang mit zahlreichen experimentellen Daten. Im Gebiet 


mittlerer Temperaturen, fiir die 0.3 << ——— < 1,8 ist, hingt dm nahezu 
Wo 


linear von der Temperatur ab. Mit wachsender Temperatur nihert sich die 


; , Oh oe. 
ausgezogene Kurve asymptotisch der punktierten Kurve | -——. Die letztere 


ae 
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Abb. 5: Die ausgezogene Kurve stellt die Funktion V. 
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gibt die Formel (68) wieder, die im Fall beliebiger Dispersion und hoher 


Temperatur anwendbar ist. 


Im Fall der Ionenkristalle mit einem Gitter vom Steinsalztyp stehen die 
longitudinalen Polarisationsschwingungen in stirkster Wechselwirkung mit 
den optischen Elektronen. Die entsprechenden Frequenzen @, fiir die Alkali- 
halogenidkristalle wurden von K. B. ToLpyGo (1950) aus Daten itiber Ab- 
sorption und Dispersion des Ultrarotlichts in diesen Kristallen bestimmt. 


Fiir KBr ergab sich hay = 0,020 eV. 
Die Anwendung der Formel (67a) auf 
die F-Absorptionsbande des KBr und 
den Vergleich der Ergebnisse mit dem 
Experiment ist in Tabelle I darge- 
stellt!). 

Die StToKEssche Verschiebung Aw 
wurde hier durch Vergleich des berech- 
neten und des beobachteten dw bei 
Zimmertemperatur bestimmt. (Unten 
im §7 wird Aq theoretisch berechnet.) 


Die Ubereinstimmung der berechneten 


Tabelle I 
2,354 & Vand w = 0,323 eV 


hd@ in eV 


0,25 
293 0,37 
473 0,48 
0,63 


berechnet gemessen 


0,22 
0,37 
0,49 
0,66 


1) Anm. d. dtsch. Red.: Zur Theorie der F-Zentren-Absorption vgl. auch DEXTER, 


Phys. Rey. 83, 435, 1951. 
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d der gemessenen d@ bei drei anderen Temperaturen dient zur Bestiitigung 
er Theorie. Entsprechende Ergebnisse erhilt man auch fiir die F-Banden 
anderer Alkalihalogenidkr istalle. 

In Abb. 6 ist die beobachtete [VLAM (1949)] Temperaturabhingigkeit der 
Halbwertsbreite der Lumineszenzbanden fiir die oben erwihnten Kristalle 
CaWO, und MgWO, dargestellt. Der Verlauf der experimentellen Kurven 
stimmt qualitativ mit der in Abb. 5 dargestellten theoretischen Kurve iiber- 
ein. Zum Zwecke eines genaueren quantitativen Vergleichs der Kurven muB 
man die Eigenfrequenzen der longitudinalen Polarisationsschwingungen der 
lonen in diesen Kristallen kennen. Der Vergleich der Abb. 6 mit Abb. 5 
zeigt, daB fiir die beiden erwaihnten Wolframate Temperaturen iiber 400° K 
schon zum Gebiet der hohen Temperaturen gehéren, fiir die man an Stelle 
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Abb. 6: Abhangigkeit der Halbwertsbreite der Lumineszenz von der Temperatur bei CaWO, und MgWO,. 


der Formel (67) die Formel (68) verwenden kann. Die letztere gibt zwar 
einen etwas zu kleinen Wert fiir dw (wie durch Vergleich der punktierten 
und der ausgezogenen Kurve in Abb. 5 ersichtlich ist); sie hat jedoch den 
Vorzug, daB in ihr die Frequenzen der Atomschwingungen tiberhaupt 
nicht vorkommen; man braucht nur die STOKESsche Verschiebung d@ zu 
kennen. 

Zur Bestimmung von Aw fiir MgWO, verwenden wir die Experimente, bei 
denen in einer und derselben Kristallprobe sowohl] das Lumineszenzspektrum 
als auch das entsprechende Anregungsspektrum gemessen wurden [FRONDA 
(1944)]. Nach diesen Messungen liegt das Anregungsmaximum in der Nahe 
von 4,8 eV und das Lumineszenzmaximum bei 2,3 eV. Die StTOKESsche Ver- 
schiebung betragt also 44m = 2,5 eV. Berechnen wir den Absolutwert von 
da bei der Temperatur 451° K nach Formel (68), so erhalten wir hO@per = 
= 0,74 eV. Der entsprechende gemessene Wert ist [VLAM (1949)]: 


hOWeem = 0,84 eV. 


Die Ubereinstimmung zwischen Theorie und Experiment muB als gut an- 
gesehen werden, besonders wenn man bedenkt, daB die Formel (68) einen zu 
kleinen Wert fiir 6 liefert. 
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Die mittlere Wirmeerzeugung beim optischen Ubergang betriigt nach Formel 


(70) ie Fir die F-Absorptionsbanden des KBr-Kristalls betragt sie 0,5 eV 


2 

also 25 Quanten der Polarisationsschwingungen der Ionen’). Wir haben es als 
tatsichlich mit einem Fall hoher Warmeerzeugung zu tun. Von der gleiche 
GréRenordnung ist die Warmeerzeugung bei der #'-Absorption auch in den 
anderen Alkalihalogenidkristallen. Bei MgWO, betrug die mittlere Warme- 
erzeugung fiir die oben betrachteten optischen Uberginge 1,2 eV; das bes 
deutet ebenfalls eine hohe Wiirmeerzeugung (bei ziemlich wahrscheinlichen 
Annahmen iiber die GréBe von ha,,). Im letzten Fall ist die Warmeerzeugung 
hdher als die halbe Energie des emittierten Quants. 
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Abb. 7: Symmetrie der Absorptions- und Lumineszenzspektren einer Lésung yon Rhodamin 6 G extra | 
in Athylalkohol bei verschiedenen Temperaturen (nach LJOWSCHIN). 


In Abb. 7 sind die Absorptions- und Emissionsspektren von Molekiilen des: 
Rhodamins 6 G extra, gelést in Athylalkohol, dargestellt. Diese Abbildung: 
demonstriert ausgezeichnet die LJOWSCHINsche Symmetrie [LJOWSCHIN: 
(1951, Seite 99)]. Aus der Symmetrie folgt, daB die Absorptions- und Lumines- | 
zenzbanden einfach sind. Bei Zimmertemperatur betriigt die Halbwerts-. 
breite beider Banden 0,22 eV. Die StoKEssche Verschiebung betriigt 0,135 eV. 
Berechnen wir nach Formel (67) die mittlere effektive Frequenz der Atom- 
schwingungen, die am stiirksten mit dem optischen Elektron wechselwirken.. 
so erhalten wir 
AO, = 013 OV. 


Die mittlere Wiarmeerzeugung betriigt etwa 0,07 eV, also weniger als ha@p, 
In unserem Beispiel ist also das Kriterium der hohen Warmeerzeugung ver-+ 
letzt. Hierdurch erklirt sich auch die starke Abweichung der Form der Bandd 
von der GAusskurve?), | 

; 


stark abweichende Krgebnis von Muro, Progr. Theor. Phys. 4, 181, 1949. 
*) Anm, d.dtsch. Red.: Zu diesem Problemkreis vgl. auch das Buch von FE. I. AprRof 
witscH: Kinige Fragen zur Theorie der Lumineszenz der Kristalle, Akademie-Verlag: 
Berlin 19538. 
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Der Fall der lonenkristalle und grofer Radien der Elektronenzustiinde 


e oben dargelegten Ergebnisse wurden ohne konkrete Annahmen iiber das 
odell des lichtabsorbierenden Zentrums und ohne eine explizite Darstellung 
er Wellenfunktion yp, () der Elektronen gewonnen. Zur theoretischen Lisung 
ieler wichtiger Probleme mu man jedoch ein Modell des Zentrums und. der 
istallstruktur vorgeben und einen expliziten Ausdruck fiir y,(r) finden. 
u diesen Problemen gehéren: die Berechnung der Frequenz des rein elek- 
ronischen Uberganges «), die theoretische Berechning der Wiirmeerzeugung 
d der STOKESschen Verschiebung, die Berechnung der thermischen An- 
egungsenergie und der thermischen Dissoziationsenergie der Elektronen des 
Zentrums usw. 

Unten werden die Ionenkristalle mit kubischem Gitter behandelt, deren 
Elementarzelle aus zwei Ionen besteht. Viele Ergebnisse gelten auch fiir 
beliebige isotrop polarisierbare Jonenkristalle und polare Fliissigkeiten. 

Von den verschiedenen Typen der Normalschwingungen stehen die longi- 
tudinalen Polarisationsschwingungen in stirkster Wechselwirkung mit den 
Elektronen. Die Wechselwirkung der akustischen und der transversalen 
Polarisationsschwingungen mit den Elektronen vernachlassigen wir. Diese 
Schwingungen werden also als konservatives abgeschlossenes Untersystem 
abgetrennt, und man kann sie in der Folge weglassen. Nehmen wir den Radius 
| ws |? der Wolke der Elektronen als hinreichend groB an, so kénnen wir den 
Kristall als dielektrisches triage polarisierbares Kontinuum auffassen. ® (x, ¢) 
sei der Anteil des Vektors der spezifischen Polarisation des Diraenns 
der zu den longitudinalen Polarisationsschwingungen gehért. Dann lat sich 
dig potentielle Energie der Wechselwirkung der Elektronen des Zentrums 
untereinander und mit den Ionen sar nea Oe schreiben : 


9 


Ver.g=SWwI+t TW td+—= Y——e Dy eee ae Dade’. (74) 
i i eae "ft =a Y —T;) 

Hierbei ist W(r) das tibliche periodische Potential des Elektrons in einem 
Kristal] mit unbeweglichen Ionen, W’ (r) die lokale Anderung des periodischen 
Potentials, die von der Stérung (dem Zentrum) hervor gerufen wird, ¢ und 7 
sind die Naimmers der optischen Elektronen des Zentrums und 1;; der Ab- 
‘stand zwischen ihnen, n? ist das Quadrat des Brechungsindex und spielt 
die Rolle der Dielektrizitdtskonstanten fiir die Polarisation der stark ge- 
bundenen Elektronen des Dielektrikums. Das letzte Glied in (74) ist die 
Wechselwirkungsenergie des Elektrons mit dem traigen Anteil der Polarisation 
des dielektrischen Kontinuum. 

Die Polarisation des Dielektrikums liBt sich in einen triigheitslosen und einen 
tragen Anteil aufteilen. Der erstere stellt die Polarisation des Kristalls dar, 
bei der die Ionen in ihren Gitterpunkten bleiben. Diese Polarisation beruht 
auf der Deformation der Elektronenschalen des Dielektrikums; sie tritt 
‘allein auf (d. h. ohne den triigen Anteil) bei einem aiuBeren Feld, dessen 
Frequenz gré8enordnungsmiBbig zwischen den ultraroten Higenfrequenzen 
der Ionenschwingungen und den charakteristischen Frequenzen der ge- 
bundenen Elektronen des Dielektrikums liegt. Die oben erwiihnten schwach 
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gebundenen optischen Elektronen des Zentrums sowie auch die Tae 
elektronen, die sich durch den Kristall bewegen, erzeugen ein Feld m 
eben dieser Frequenz und rufen infolgedessen die trigheitslose Polarisati 
der Ionen hervor. Diese Polarisation folgt tragheitslos der Bewegung des 
Elektrons, das sie erzeugt; infolgedessen wirkt auf das Elektron eine Zusatz- 
kraft von der Periodizitat des Gitters. Wir nehmen an, diese Kraft sei bereits 
in dem periodischen Potential W (x) einbegriffen, das in Formel (74) auftritt. ; 
Unter dem tragen Anteil der Polarisation verstehen wir die Polarisation, die 
auf den Verschiebungen der Ionen aus ihren Gitterpunkten beruht. Sie setzt 
sich zusammen aus den Dipolmomenten, die den Ionenverschiebungen ent- 
sprechen, und aus dem Teil der Polarisation der gebundenen Elektronen, der 
durch diese Ionenverschiebung erzeugt wird. Eben um diese trage Polarisation 
handelt es sich, wenn von den (thermischen) Polarisations-Eigenschwingungen 
der Ionen die Rede ist. Kine triage Polarisation kann man im Kristall erzeugen, 
wenn man zunichst ein diuBeres elektrisches Feld langsam einschaltet, und 
es dann so schnell wieder ausschaltet, da zwar die Ionen keine Zeit haben, 
ihren Ort zu verindern, die stark gebundenen Elektronen aber imstande sind, 
wahrend des Schaltens ihre Polarisation zu indern. Die Polarisation, die gleich 
nach diesem Abschalten des Feldes im Kristall verbleibt, stellt den trigen 
Anteil der Polarisation dar. Variiert die elektrostatische Induktion 9 (tr) 
des eingeschalteten Feldes riumlich hinreichend langsam (indert sie sich auf 
Entfernungen von der Gréfenordnung der Gitterkonstanten nur wenig), sa 
kann man die makroskopische Elektrodynamik anwenden. Beim langsamen 
Einschalten des Feldes entsteht dann eine Gesamtpolarisation 


Bee aD (75° 


(¢ ist die Dielektrizitatskonstante des Kristalls).  _ 
Beim ploétzlichen Abschalten des Feldes betriigt die Anderung der Polarisation 


n? —] 
éf = — ——_.- 76 | 
4 an (? 
Die verbleibende trige Polarisation betrigt also 
C 
% = Bees a Oy = ®, (77° 
4a 
+ 
mi & 1 1 Ad 


Die potentielle Energie der triigen Polarisation des Kristalls ist leicht z: 
berechnen als die Arbeit, die das Feld D beim langsamen Einschalten un 
plotzlichen Abschalten leistet. Diese potentielle Energie betriigt 


Uy = —" [ Bede. (79 


Das Integral erstreckt sich tiber das gesamte Kristallvolumen, 
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ir miissen betonen, daf das letzte Glied in Formel (74) die potentielle 
ergie der Wechselwirkung der Elektronen allein mit dem trigen Anteil 
er Kristallpolarisation 8 (r, ¢) darstellt. Die Wechselwirkung mit der trig- 
eitslosen Polarisation ist, wie schon gesagt, in das erste Glied der Formel (74) 
it einbegriffen. 

ir die weiteren Rechnungen mu8 man die Normalkoordinaten der Schwin- 
ingen des Dielektrikums explizit einfiihren. Hierzu entwickeln wir die triige 
olarisation in eine trigonometrische Reihe von der Form 


x 2 ree rane apie Fen (0) 
ry — Nv —— & = — 
B (t,t) = [el Pa (t) te), Xu (0) \3 ‘a xx fiir x, > 0 coe 


(ZL ist die Abmessung des kubischen Grundgebietes, in dem y, orthonormiert 
ist). Hier und im folgenden ist x ein dreidimensionaler Vektor. Als Normal- 
<oordinaten kann man bekanntlich die Entwicklungskoeffizienten P,, ver- 
wenden. Die Energie der Eigenschwingungen des Dielektrikums wird dann 
folgendermaBen dargestellt : 


2 a 
Se (Ps ae Ps), (81) 
zx @. 


Hierbei sind w, die Eigenfrequenzen der Schwingungen und c,, positive Kon- 
stante. Fiir die langwelligen Schwingungen, bei denen sich die Polarisation 
raumlich langsam andert, mu8 das erste Glied in Formel (81) — die poten- 
tielle Energie der Ionen — mit der Gré8e (79) tibereinstimmen. Hieraus folgt, 
daB bei x + 0 die GréBe c, gegen die Konstante cy geht, die in Formel (78) 
vorkommt. 

Um den in § 2 entwickelten mathematischen Apparat unmittelbar anwenden 
zu kénnen, muB man die Normalkoordinaten so normieren, daB der Energie- 
operator (81) die Form 


1 oP 
— Viho,|g—-+—}. 82) 
2 a! h Ox (« aa) ( ) 


et | AO 83 
Ux \ rae (35) 


. “= a 7 e — | 
P= 210 Gt iol += 2.—, (84) 
i tt 4) 
v9 
—- —- %u (t’) (¥ Y;) 
hOnlx x : 
Vv.) =—e/ eae o—..8 dt’ . (85) 


1) Anm. d. dtsch. Red.: 7 ist Abkiirzung fiir alle t;; s. 8. 372. 
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Die weiteren Rechnungen fiihren direkt zu den Formeln in § 2. Nach Forme 
(13) erhalten wir namlich 


h ee 
Ve =—J pa Does (8 
wobei : 
; x P 
Dus =| De Te] )d (87) | 
7 
und. 


, 
dt,...dt%...dty.  (88)} 


ist. Der Vektor D, ist die elektrostatische Induktion, die durch die |y,|?-Wolke« 
der optischen Elektronen erzeugt wird. Der Skalar D,, ist der Entwicklungs-| 
koeffizient fiir die Entwicklung von 9, in eine trigonometrische Reihe von 
der Form (80). 

Setzen wir (86) in (28) ein, so erhalten wir 


h? | appa : 
J (ps) = i niin Im Ai Vo (0| vs (7) dr ——— Sen Lins. (89) 


Die Wellenfunktion ys, der optischen Elektronen wird durch Extremums 
bildung des Funktionals (89) bestimmt. 

Ist der effektive Radius der |y,|?-Wolke der optischen Elektronen hinreichend 
eroB (beispielsweise so gro wie die Gitterkonstante oder gréBer), so kanm 
man in dem Funktionai (89) folgende Vereinfachungen vornehmen: Unter 
Benutzung der Methode der effektiven Masse des Elektrons im Kristali®” 
(PEKAR (1946a, 1951)| kann man in der HAMILTONfunktion das periodische 
Potential W (r;) fortlassen, wobei man gleichzeitig die Elektronenmasse m7 
durch eine bestimmte effektive Masse mw ersetzt. Ferner ist, weil bei groSer 
Radien der Zustiinde die Induktion D, eine langsam veriinderliche Funktion} 
der Koordinaten ist, ihr FOURIER-Koeffizient D,, nur fiir die langwelliger} 
Oberschwingungen von Bedeutung, also fiir kleine |%|. Im letzten Glied dew 
Ausdrucks (89) iiberwiegen also in der Summe die Glieder mit kleinem | |e 
In diesem Fall kann man den Koeffizienten c, mit seinem Grenzwert ¢y vod 
das Summenzeichen ziehen. Infolge dieser Vereinfachungen kann man da 
Funktional (89) foleendermafen schreiben : 


: A 
ities ws (Pr) t- a eee BW eae 


a fl t t 
| 
+5 2a] inde 2 [or wae (oc 
NW iS; Taj 
Nachdem wir ys (r) durch Extremumsbildung des Funktionals (90) gefunded| 


haben, miissen wir nach den Formeln (88), (87) und (86) entsprechend auc( 1 


1) Anm. d. dtsch. Red.: Zu dieser Methode s. ADAMs, Phys. Rev. 85, 41, 1952 wd 
Winston, Phys. Rev. 89, 317, 1952. 1 | 
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D;. D,,, und V,,, berechnen. Diese Rechnungen miissen sowohl fiir den Aus- 
gangs- als auch fiir den Endzustand des optischen Ubergangs ausgefiihrt 
werden. Dann kann man nach Formel (70) die Gleichgewichtslagen der 
Tonen q,5, und q,, berechnen, ferner nach den Formeln (53), (54), (57) und (58) 
die Hauptparameter der Absorptionsbande: Die Lage des Maximums, den 
Maximalwert des Absorptionskoeffizienten t,, die Halbwertsbreite der 
Bande 6@. Die StoKEssche Verschiebung Aw und die Wirmeerzeugung beim 
optischen Ubergang werden nach Formel (70) berechnet. Bei diesen Rech- 
nungen kann man @, niiherungsweise durch @, ersetzen, denn in den Tonen- 
kristallen des betrachteten Typs gibt es stets einen Zweig longitudinaler 
Polarisationsschwingungen, in dem die Dispersion der Frequenzen w, ge- 
wohnlich gering ist. AuBerdem spielen bei grofBen Radien der Elektronen- 
zustinde in simtlichen oben angegebenen Formeln die Glieder mit kleinen | x | 
die Hauptrolle. 

Die Frequenz @® des reinen Elektroneniiberganges und die thermische An- 
regungsenergie W,, werden nach den Formeln (31) und (73) berechnet. 
Das oben umrissene Programm der Rechnungen ist fiir alle Typen lokaler 
Zentren gemeinsam. Vom Typ des Zentrums, d.h. von seinem Modell, hingt 
nur die Form des Potentials W’ (r;) und die Anzahl der optischen Elektronen 
ab. Bis heute wurden drei Typen lokaler Zentren durchgerechnet?): 


1. Die F-Zentren [PEKAR (1950a), PEKAR und DEIGEN (1948)]. Eine Be- 
griindung des Modells der F-Zentren wird von PEKAR (1947) gegeben. In 
diesem Fall gibt es ein optisches Elektron: 


2 


ze 


Ww’ = — — : 9] 

haar (91) 
2. Die F’-Zentren wurden vom Verfasser gemeinsam mit O. F. TOMASSE- 
WITSCH untersucht [PEKAR und TOMASSEWITSCH (1951, 1952); TOMASSE- 
WITSCH (1951, 1952)]. In diesem Fall gibt es zwei optische Elektronen. W’ stellt 
ebenso wie im vorhergehenden Fall das Potential eines anziehenden COULOMB- 
schen Zentrums dar: 


Yr) = — ae ey (92) 
3. Die F,-Zentren, d.h. gepaarte F-Zentren, die miteinander wechselwirken 
und im Kristall eine Art zweiatomiges Molekiil bilden, wurden von M. F. 
DEIGEN (1951) untersucht. Auch in diesem Fall gibt es zwei optische Elek- 
tronen: 


2e2 l ] é 
W* (tj) = —- = ( Or [u— a ee) 


r; | 


wobei a der Abstand zwischen den positiven Zentralladungen der #-Zentren 
ist. Der Kiirze halber behandeln wir hier nicht alle drei der oben erwiihnten 


1) Anm. d. dtsch. Red.: Modelle fiir F-Zentren behandeln auch KRUMHANSL und 
| SCHWARTZ, Phys. Rev. 89, 1154, 1953; PINCHERLE, Proc. Phys. Soc. (A) 64, 648, 1951; 
Smmpon, Proc. Roy. Soc. (A) 197, 269, 1949. Dort weitere Zitate. 
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Typen von Zentren, sondern erortern nur die Theorie der F-Zentren al 
typisches Beispiel. 


~ 


7. Theorie der F-Zentren und Vergleich mit dem Experiment 


Eine Variante der Theorie der -Zentren, in der die Bewegung des Electro 
quantenmechanisch, die der Jonen dagegen klassisch behandelt wird, wur o 
vom Verfasser gemeinsam mit M. F. DEIGEN entwickelt [PEKAR (1947); 
PEKAR und DEIGEN (1948)}. Die Energieterme und Wellenfunktionen des 
Elektrons im F-Zentrum wurden bestimmt; es gelang jedoch nicht, die Form 
der F-Absorptionsbande, ihre Temperaturabhangigkeit und ihre Halbwerts- 
breite abzuleiten. Diese Variante der Theorie entspricht dem oben umrissenen 
allgemeinen Rechenschema nicht ganz. Wir werden auf sie deshalb nicht 
niher eingehen. | 
Hier stellen wir kurz die Ergebnisse einer spiteren Arbeit des Verfassers | 
zusammen [PEKAR (1950a, 1951)], in der die F-Zentren nach dem im vorigen | 
Paragraphen entwickelten allgemeinen Rechenschema behandelt werden. 
Der Grundzustand des Elektrons wird durch absolute Minimumsbildung des: 
Funktionals (90), (91) mit Hilfe einer direkten Variationsmethode bestimmt. . 
Die Wellenfunktion des Elektrons wird durch den Ausdruck 


g's 
Y;, = Cran ene (94) 


i 
Via. ; 


approximiert. Der Parameter a wird so gewihlt, dafi J [ys,| ei Minimum 
wird. Es ergibt sich 


be? faz 
=Fale +4) q 
pet (2 Cae 
din = i 
; Ae *) i 


Auch andere Approximationen fiir ws, wurden durchprobiert [PEKAR unc 
DEIGEN (1948) |, jedoch erwies sich die Niherungsformel (94) als die giinstigste 
sie fiihrt zum tiefsten Wert fiir das Funktional J {ws,|. Von den angeregter| 
Zustinden des Klektrons wurde derjenige berechnet, in den der optischal 
Ubergang des Elektrons mit gréBter Wahrscheinlichkeit er folgt. Die Wellens} 
funktion dieses Zustandes w urde durch den Ausdruck | 


8 | 
We, = —=—= (af) re—* Ft cos} (97° 


V2 Tt 
approximiert. Diese Funktion ist orthogonal zu der Funktion (94) des Grundé 
zustandes, wie das in der direkten Variationsmethode auch sein muf. De? 
oral H ace ane der Bedingung des absoluten Minimums von J [Weill 
sestimint, Die Ree gen Vv 3 iRM é ur 
bet echnungen von G, J. SILBERMANN (1949a) fiir den Fai 
1+ 0.8914 6¢9 


ee EC 


2p = 
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get fil 2 

Ups a (= + 0,3914 “} (99) 
Die weiteren Rechnungen wurden ohne neue Niherungsansiitze nach den 
allgemeinen Formeln des vorigen Paragraphen durchgefiihrt. 
Die Theorie enthilt vier Parameter: e, n, w) und «. Diese Parameter sind 
nicht nur spezifisch fiir die #-Zentren, die Polaronen (siehe unten) und die 
uibrigen lokalen Zentren, sondern charakterisieren auch den Idealkristall ohne 
lokale Elektronenzustiinde. Die Parameterwerte lassen sich also unabhiingig 
von den oben behandelten Erscheinungen durch Untersuchung verschiedener 
Higenschaften der Idealkristalle bestimmen. Insbesondere gehen die Para- 
meter ¢, n und w, wesentlich in die Formel ein, die die Dispersion des Lichts 
im Idealkristall beschreibt, und lassen sich aus den experimentellen Dis- 
persionskurven bestimmen, falls diese bekannt sind. Dagegen gibt es bis 
heute keine zuverlissige Methode zur unabhiingigen Bestimmung des Para- 
meters uw. Die effektive Masse des Elektrons bleibt also der einzige unbekannte 
Parameter der Theorie. Weiter unten wird w durch Vergleich der berechneten 
Abzisse (Frequenz) des Maximums @,, der F-Absorptionsbande mit dem 
entsprechenden experimentell gemessenen Wert bestimmt. 
Tabelle IT gibt einen Vergleich der berechneten mit den experimentell ge- 
messenen thermischen Dissoziationsenergien der F-Zentren'). Diese Werte 


Tabelle IL 


Cu,0 9 4 0),62 1,81 0,51 0,5—0,6 
NaCl 5,8 2,33 2,65 2,78 22; 1,9 
NaBr 6,39 2,6 2,29 2,96 1,9 1,6 
NaJ 6,6 2,91 Peal 3,25 1,8 1,5 
KCl 4,78 Pea lps: 2,19 1,85 2,0 2,0 
KBr 4,81 2,36 1,96 1,87 1,8 vi 
KJ 5,2 2,65 1,8 2,11 Lee en 
RbCl 5,2 2,18 1,98 1,78 1,7 sé 
RbBr 5,16 2,36 A Gail 1,70 1,5 1,4 
RbJ 5,58 2,6 1,59 1,89 1,4 3 


wurden an Hand der Temperaturabhingigkeit der elektrischen Leitfaihigkeit 
der Kristalle bestimmt, in denen die F-Zentren als Spender fiir die Ladungs- 
triger wirkten. Der Beweis dafiir, daB in diesen Kristallen tatsichlich die 
F-Zentren die Spender fiir die Ladungstriger darstellen, liegt in dem Umstand, 
daB eine Entfernung der F-Zentren aus dem Kristall zu. einer Abnahme der 


1) Anm. d. dtsch. Red.: Andere Arbeiten zu diesem Problem: GooDMAN, LAWSON und 
Scutrr, Phys. Rev. 71, 191, 1947; Kuso, Phys. Rev. 86, 929, 1952 sowie HUANG und 
Ruys l. ¢. 

27 Zeitschrift ,,Fortschritte der Physik‘ 


f 
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elektrischen Leitfaihigkeit um viele GréBenordnungen fiihrt. In den betrach. 
teten Kristallen wird die Temperaturabhingigkeit der elektrischen Leit 


fahigkeit gut durch die Formel . 
W : 
=e 27. a 


wiedergegeben. Der hier auftretende Parameter W la8t sich experimentell 
leicht bestimmen. Ihn meinen wir, wenn wir von dem experimentellen Wert 
der thermischen Dissoziationsenergie der #’-Zentren sprechen, In der Tabelle IJ 

sind die Werte W8™ aufgefiihrt, die SMAKULA (1934) fiir die Alkalihalogenid- 
kristalle erhielt. 

In der ersten Spalte der Tabelle II sind die Werte fiir das Cu,O angegeben, 
Hier bildet das F-Zentrum ein positives ,,Leitfihigkeitsloch“‘, das sich in der 
Nihe eines leeren Kupfergitter- | 
punkts lokalisiert hat. Das Ab- | 
sorptionsspektrum des F-Zen- | 
trums im Cu,O liegt im Ultra. - 
roten. Diese Absorption ist ver- 
bunden mit einer Photoleitfahig- . 
keit. SCHONWALD (1932) unter- . 
suchte die spektrale Verteilung : 
dieser Photoleitfahigkeit; seine : 
Ergebnisse sind in Abb. 8 dar- - 
gestellt. 

Das Maximum der Kurve ent- - 
spricht einer Energie des Licht- - 

quants von 0,62 eV. Die Kurve : 
hat die typische Glockenform, die : 
auch fiir die F-Zentren in Alkali. - 
halogenidkristallen kennzeich- 

nend ist [GYULAI (1925); ARSSE- | 
NEWA (1926)]. 

Die Temperaturabhangigkeit deri 
elektrischen Leitfahigkeit im¢ 


LO NOD TOSI OT™ (CSM Ose OSC a eV 


Abb. 8: Spektrale Verteilung der Photoleitfaihigkeit CuO wurde von vielen Forschern 
von Cu,O, die auf tibersehiissigem Sauerstoff = 1 = 
beruht, untersucht ). Die gewonnenen ti 


Werte fiir die thermische Disso-» 
ziationsenergie streuen etwas. In Tabelle II sind die Angaben von ENGEL-. 
HARD (1933) benutzt, weil sie an denselben Cu,O-Proben gewonnen wurden! 
und sich auf die gleichen Zentren beziehen wie die spektralen Messungeni 
von SCHONWALD. 

Die in der Tabelle If angegebenen berechneten Werte der thermischen Disso-~ 
ziationsenergie wurden unter der Voraussetzung abgeleitet, daB sich infolge: 
der Dissoziation des F-Zentrums ein freies Polaron bildet (siehe unten) ,) 
und daB das F-Zentrum kein Ladungstriger ist [PEKAR (1950a, 1951)]. 


‘) Anm. d. dtsch. Red.: Neueste Messungen: BLANKENBURG und SCHUBART, Ann. di 
Physik 12, 281, 1953. Dort weitere Zitate. 
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Die zur Berechnung notwendigen Werte der Parameter ¢, n und w) wurden 
aus den experimentellen Dispersionskurven fiir die betreffenden Kristalle 
entnommen. Fiir die Alkalihalogenidkristalle wurden diese Parameter von 
K.G. ToLpyGo (1950) ziemlich genau bestimmt. Fiir das Cu,O ist die experi- 
mentelle Dispersionskurve leider nur in einem Spektralbereich von ungeniigen- 
der Gré8e bekannt. In diesem Fall liegen deshalb ¢ und n weniger genau fest; 
dagegen wurde @) aus der Temperaturabhingigkeit der Beweglichkeit der 
Ladungstriger bestimmt (siehe unten). 


Durch einen Vergleich der beiden letzten Spalten der Tabelle II erkennt 
man, da die berechneten Werte fiir die thermische Dissoziationsenergie gut 
mit den gemessenen iibereinstimmen. Schlechter ist die Ubereinstimmung bei 
den drei Natriumsalzen, bei denen der Radius der |p |?-Wolke des Elektrons 
im F-Zentrum zu klein‘ist, so da8 das Kriterium fiir die Anwendbarkeit der 
Theorie verletzt ist. 

Eine sehr charakteristische Eigenschaft unserer Theorie besteht darin, daB 
die Energie der optischen Anregung hw» gréBer ist als die thermische Disso- 
ziationsenergie W, obwohl beim optischen Ubergang das Elektron in einen 
diskreten Term, dagegen bei der thermischen Dissoziation in das kontinuier- 
liche Spektrum iibergeht. Das liegt daran, da& beim optischen Ubergang 
neben der Anregung des Elektrons ein betrichtlicher Teil der Energie des 
Lichtquants fiir Wiairmeerzeugung aufgewandt wird. Aus Tabelle II erkennt 
man, daB das erwihnte Uberwiegen der Energie des optischen Uberganges 
tiber die thermische Dissoziationsenergie tatsachlich experimentell beobachtet 
wurde. 

In §5 dieses Artikels wurde schon gezeigt, daB die berechnete Form und 
Temperaturabhangigkeit der F-Absorptionsbanden gut mit dem Experiment 
iibereinstimmen. Dort wurden aber nicht die absoluten Werte der Halb- 
wertsbreiten der Absorptionsbanden berechnet, da hierzu ein expliziter Aus- 
druck fiir die Elektronenwellenfunktion yw, (r) nétig ist. In diesem Paragraphen 
kénnen wir, nachdem wir diese expliziten Ausdriicke fiir y,, und ys. gewonnen 
haben, die absoluten Werte fiir die Halbwertsbreiten und die mittlere Warme- 
erzeugung beim optischen Ubergang berechnen. 


In Tabelle III werden die berechneten Halbwertsbreiten der /-Absorptions- 
banden fiir Zimmertemperatur [PEKAR (1950a, 1951)| mit dem Experiment 
verglichen. Der Parameter a wird durch Formel (50) definiert. Die Form der 
F-Absorptionsbanden wurde mit Hilfe der Formeln (47) und (51) bestimmt. 
- Aus der Tabelle geht hervor, daB bei den Alkalihalogenidkristallen die theore- 
tischen Werte fiir die Halbwertsbreiten gut mit dem Experiment iiberein- 
stimmen. Beim Cu,0 ist die gemessene Halbwertsbreite um 0,1 — 0,15 eV 
hdher als die berechnete; dies kann man folgendermafen erkliren: die oben 
‘erwihnte Streuung der experimentellen Werte fiir die thermische Disso- 
ziationsenergie der F-Zentren im Cu,O beruht auf der Nichtidealitat der 
Kristalle und einer gewissen Verschiedenheit der F-Zentren. Nimmt man an, 
daB bei den oben erwihnten Cu,O-Proben die Energieterme der F-Zentren 
um etwa 0,1 — 0,15 eV streuten, so ist damit gleichzeitig auch die Streuung 
| der thermischen Dissoziationsenergien und die beobachtete zusiitzliche Ver- 
| breiterung der F-Bande erklirt. 

; 27* 
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Tabelle III 4 
T = 293° K | 


hA ow (eV) 
howber (eV) 0,47 0,40 0,37 0,35 
6 wsem (eV) 


Im Grenzfall tiefer Temperaturen diirfen die Form der Absorptionsbanden } 
und insbesondere ihre Halbwertsbreite nicht von der Temperatur abhingen. | 
Das Experiment bestatigt diese SchluBfolgerung der Theorie. 

In Tabelle IV sind die berechneten 
und die gemessenen Halbwertsbreiten 
der F-Absorptionsbanden im Grenz- 
fall tiefer Temperatur angegeben. 
Der beschrinkte Umfang  unseres 
Artikels zwingt uns hier, auf eine 
weitere Darlegung des Vergleichs von d wPer (eV) 0,24 
Theorie und Experiment zu verzichten ; dm8em (eV) 0,22 
wir weisen den Leser auf die ein- 
schlagige Literatur hin. 

Die Behandlung der Rekombination der Leitungselektronen (der Polaronen) | 
mit den #’-Zentren und die Bildung der #”-Zentren, die Berechnung der ent-: 
sprechenden freien Weglinge fiir das Haften des Elektrons und den “Vergleich! 
dieser theoretischen Uberlegungen mit dem Experiment findet man bei: 
PEKAR (1950b, 1951); PEKAR und PERLIN (1950). Eine Erklairung fiir die¢ 
Verschiebung des Maximums der -Absorptionsbande ins Rote bei Erhitzung¢ 
des Kristalls und ein Vergleich zwischen Theorie und Experiment wird von! 
PEKAR (1950a, 1951) cepoher 

Die Quantenausbeute der Photodissoziation der F- und F’-Zentren, ihre¢ 
Temperaturabhingigkeit, ein Vergleich zwischen Theorie und Experiment und: 
eine Reihe anderer Probleme sind von PEKAR (19% 51) und SILBERMANN (1949b)) 
behandelt worden. 


Tabelle IV 


8. Die Polaronen 


Wir gehen nun zur Behandlung der Quantenzustiinde des Leitungselektrons: 
im idealen Ionenkristall ohne Stérstellen der Periodizitiit, Beimischungen! 
usw. tiber. Urspriinglich erfolgten allerdings die Untersuchungen in um- 
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gekehrter Reihenfolge: zuerst wurde das Leitungselektron im Idealkristall 
behandelt und dabei der mathematische Apparat und das System der Niherun- 
gen entwickelt, die dann zur Behandlung der an Kristallfehlstellen und Fremd- 
atomen lokalisierten Atome verwendet wurden. Diese Reihenfolge hielt der 
Verfasser auch in seiner Monographie [PEKAR (1951)] ein, die der eingehenden 
Darstellung der Originaluntersuchungen gewidmet ist. Hier jedoch haben 
wir den anderen Weg eingeschlagen: das Elektron im Idealkristall wird als 
Spezialfall des Problems des Elektrons in einem Kristall mit FehIstellen auf- 
gefaBt, das in den vorigen Paragraphen abgehandelt wurde. Hierdurch er- 
zielt man eine gréBere ‘Knappheit der Darstellung, mu aber einige Opfer 
in bezug auf Strenge und Uberzeugungskraft bringen. 

Da im betrachteten Fall keine Periodizitaétsstérungen vorhanden sind und 
auch nur ein einziges Elektron behandelt wird, hat man in Formel (74) das 
zweite und dritte Glied fortzulassen. Nimmt man an, da der Radius des 
Elektronenzustandes hinreichend groB ist, und benutzt man die Methode der 
effektiven Masse, so kann man in Formel (74) auch das erste Glied weglassen 
und hat dafiir dem Elektron eine effektive Masse « zuzuschreiben [PEKAR 
(1946a, 1951)]. Die potentielle Energie der Wechselwirkung des Elektrons 
mit den Jonenschwingungen sieht dann folgenderma8en aus: 


V (t,q) = lip Pe tei (101) 


Hierbei ist $$ (r’) der triage Anteil der spezifischen Polarisation des Dielektri- 
kums, der auf den longitudinalen Polarisationsschwingungen beruht. V,, 
wird durch die Formeln (80), (83) und (85) definiert. Die HAMILTONfunktion 
des Systems (6) kann man nun folgendermafen schreiben: 


g2 
— |. 102 
aay 


In der gewéhnlichen Bandertheorie wird das zweite Glied der HAMILTON- 
funktion — die Wechselwirkung des Elektrons mit den Schwingungen — als 
kleine Stérung angesehen und in nullter Naiherung weggelassen. Dadurch 
zerfallt das System in zwei nicht miteinander wechselwirkende Untersysteme: 
das freie Elektron und das System der unabhiingigen harmonischen Oszilla- 
toren; die Wellenfunktion des Systems ist das Produkt der entsprechenden 
beiden Faktoren. Dann wird das Wechselwirkungsglied als kleine Stérung 
eingefiihrt und in erster Niherung die Streuung der freien Elektronen an den 
-Tonenschwingungen berechnet. Durch diese Streuung werden die freie Weg- 
lange des Leitungselektrons und seine Beweglichkeit bestimmt. Die Be- 
trachtung der héheren Niherungen der St6rungstheorie sollte nur kleine 
Korrekturen liefern. Der qualitative Charakter der Lésung sollte der gleiche 
wie in nullte Naherung bleiben. Andernfalls wiire die Anwendung der Stérungs- 
rechnung selbst nicht gerechtfertigt. 

Es bestehen ernste Einwinde gegen die Anwendbarkeit der Stérungstheorie 
in dem betrachteten Fall. Wir werden auf diese Einwinde nicht naher ein- 
gehen, sondern gleich zur Darlegung einer anderen Methode tibergehen, in der 


A 1 
=> A+ 3, z(t) Ge +5 Shox (ae — 
2u Z x ; 
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die Wechselwirkung des Elektrons mit den Schwingungen nicht als klein 
angesehen wird. | 
Zur vorliufigen Orientierung und zur Abschiitzung der GréBenordnung be- 
nutzen wir unmittelbar die in den §§ 2 und 6 entwickelte Theorie. Die Wellen- 
gleichung des Systems ersetzen wir durch das aquivalente Variationsprinzip 
und lésen das Problem zunichst mit Hilfe einer direkten Variationsmethode. 
Approximieren wir die Wellenfunktion des Systems durch einen Produkt- 
ansatz (25), so erhalten wir fiir den Schwingungsanteil der Wellenfunktion 
den Ausdruck (22). Der Elektronenanteil der Wellenfunktion wird bestimmt 
als Extremale des Funktionals (90), das man in dem betrachteten Fall folgen- 

dermafen schreiben kann: | 


i 


i C r , 
os) = Ei [| Vw, |2dr — = | D? [r’, ws] dt . (103) | 


= | 


® wird durch Formel (88) definiert und betraigt im Fall eines einzigen Elektrons 


12 (e! —_ 
D [r’, pe] = { ue C a 2 dt. (104) 


Der Grundzustand y,) des Systems ist die dem absoluten Minimum ent- 
sprechende Extremale des Funktionals J [ws]. Sie wird bei PEKAR (1946b, c) 
bestimmt. Die dort ausgefiihrten Rechnungen liefern 


Wo (7) & 0,1229 ats (1 + ayr + 0,4516 a2r2) e—%", (105) 
mit 
e2 
Oy = 0,6585 SS Co. (106) 


Der entsprechende Minimalwert des Funktionals betrigt 


ret 


J =e — 0.0844 on 


C2. (107) 


Zum Zwecke grdberer Rechnungen kann man zur Minimumsbildung des | 
Funktionals J [ys] nicht die dreigliedrige, sondern eine zweigliedrige Nahe- | 
rungsfunktion verwenden: 


al 
Wo (r) & —— (1+ ar) e-*", (108) | 

V7 mu 

Man erhalt damit 
pe? 
a > R Cy, (109) | 
J, = —0,0536 #& 

So) a as IID he Co. (110) ) 


Es ergibt sich also eine Wellenfunktion des Elektrons vom abklingenden Typ, , 
d. h. ein diskreter negativer Energieterm H, des Elektrons. Berechnet man ! 


| 


| 
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die Gleichgewichtslagen der Ionen q,. nach Formel (17), ebenso die ent- 
sprechende Gleichgewichtspolarisation des Kristalls > (v), so erhalt man: 


Bo (t) = FD [t, wo). (111) 


Durch Vergleich dieser Formel mit (77) iiberzeugt man sich leicht, daB 8, 
gerade diejenige Polarisation ist, die ein aiuBeres elektrostatisches Feld mit 
der Induktion D erzeugen wiirde. 

Der gewonnenen Lésung kann man folgende anschauliche Deutung geben: 
der Kristall mit der lokalen Polarisation S8, (r) bildet fiir das Elektron einen 
Potentialtopf V (r, g)), dessen Form man nach Formel (101) berechnen kann. 
In diesem Topf gibt es diskrete Eigenwerte fiir das Elektron. Ein Elektron, 
das sich in einem dieser Terme lokalisiert hat, erzeugt ein elektrisches Feld, 
das die mittlere Induktion D besitzt und durch Formel (104) bestimmt wird. 
Diese Induktion halt die lokale Polarisation 8) (r) des Kristalls stationir auf- 
recht. Die y-Funktion des Elektrons und die ee are des Kristalls be- 
dingen einander und sind also ,,self consistent‘. Dieser ,,self consistent‘ 
Zustand des Elektrons im trage polarisierten dielektrischen Medium bezeich- 
nete der Verfasser als ,,Polaron‘‘ [PEKAR (1946b, c)]. 

Ein Mangel der oben dargelegten Lésung besteht in der ,,gewaltsamen“ 
Approximation der Wellenfunktion des Systems in Form des Produkt- 
ansatzes (25). Infolgedessen wird der mit Hilfe der direkten Variations- 
methode berechnete Energieterm des Systems etwas zu groB. Diese Un- 
zulinglichkeit wird weiter unten weitgehend behoben werden. 

AuBer dem oben behandelten Grundzustand des Elektrons im Polarisations- 
potentialtopf V (r,q)) gibt es umendlich viele angeregte diskrete Terme 
{[PEKAR (1951)]. Es sind optische Uberginge des Elektrons aus dem Grund- 
zustand in die angeregten Zustiinde moéglich. Setzt man voraus, da} wihrend 
der kurzen Zeit des optischen Ubergangs die schweren Tonen nicht imstande 
sind, sich merklich zu verschieben, so hat man den optischen Ubergang bei 
einem gegebenen festen Potentialtopf zu behandeln, der der Polarisation 
BR, (r) des Ausgangszustandes des Elektrons entspricht. Die Absorption durch 
die pata wird von PEKAR und DEIGEN (1948), PEKAR (1951) dergestalt 
behandelt. Die Formel (18) fiir den Polaronenzustand laBt sich folgender- 
maBen schreiben: 


J, = E, (Ys) + : n= | Bit )dt. (112) 


Hier ist /, (q,) der diskrete Energieterm des Elektrons, der in Gleichung (9) 

auftritt. Das zweite Glied in (112) stellt nach Formel (79) die potentielle 

Energie U, des polarisierten Kristalls (oder die Polarisierungsarbeit) dar. Wie 

im Grundzustand des Polarons, so gilt auch in jedem anderen ,,self consistent“ 
, Zustand dieses Teilchens [PEKAR (1951)]: 


ypu 2 J; (113) 


| 
; und aus (112) folgt 
J LE ipo) gee (114) 


= 
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Die Gesamtenergie des ganzen Systems wird durch Formel (21) gegeben. 
Wir betrachten die oben gewonnene Lésung des Polaronenproblems als grobe 
Orientierung und wenden uns nun der Behandlung der exakteren Lésung zu, 
Aus der Form des Funktionals (103) geht hervor, dab, wenn ws (r) eine Extre- 
male ist, die Funktion ws (t —§), wobei § ein konstanter beliebiger Vektor 
ist, ebenfalls eine Extremale ist und der gleichen Energie des Systems ent- 
spricht wie ys (t). Physikalisch entspricht dies dem Umstand, daB sich das 
Polaron mit gleicher Wahrscheinlichkeit an jedem Punkt des homogenen 
dielektrischen Mediums lokalisieren kann, wobei die Energie des Systems 
jedesmal dieselbe ist. Daraus folgt, da sich das Polaron allmahlich durch 
das Dielektrikum bewegen kann und dak fiir diese Bewegung keine Energie 
aufgewandt zu werden braucht: sie hangt mit der Tragheit zusammen; wir 
werden sie als Translationsbewegung bezeichnen. Bei dieser Translation ver- 
schiebt sich der Polarisationspotentialtopf wellenartig im Dielektrikum, und 
die y-Funktion des Elektrons, das der Bewegung des Potentialtopfes adia- 
batisch folet, fithrt ebenfalls eine Translationsbewegung aus. 
Der Hauptmangel des Produktansatzes (25) fiir die Naherung besteht darin, 
da8B er der Korrelation zwischen den Ionenschwingungen und der Elektronen- | 
bewegung nicht Rechnung tragen kann. Wahrend der Ionenschwingungen | 
variiert die Form des Polarisationspotentialtopfes um eine Gleichgewichts- . 
form V (x —&, q). Die y-Funktion des Elektrons, die dem Potentialtopf / 
adiabatisch folgt, deformiert sich ebenfalls. Das bedeutet, da die Form der * 
y-Funktion des Elektrons von den Koordinaten der Ionen abhingen mub. . 
Diese Abhingigkeit wird in der Naherungsformel (25) nicht beriicksichtigt. 
Wir verzichten deshalb jetzt auf den Produktansatz (25), behalten aber den 1 
Teil der Theorie aus § 2 bei, der sich auf die adiabatische Niherung griindet. . 
Nach dieser Naherung | PAULI (1933)] kann man das Problem in zwei Schritte » 
aufteilen. Zunichst hat man Energie und Wellenfunktion des Elektrons bei j 
beliebiger fester Ionenkonfiguration zu bestimmen. Den zweiten Schritt{ 
bildet die Behandlung der Ionenbewegung, wobei die Gesamtenergie Hs, 
als Summand in die potentielle Energie der Ionen eingeht, wie aus Formel (16) } 
ersichtlich ist. 
In der ersten Etappe dieses Verfahrens kann man die Wellengleichung, die: 
die Bewegung des Hlektrons bestimmt, durch das dquivalente Variations- - 
prinzip ersetzen. Man erhilt dann 

sa] = 


2 dt — [xo D [v, wes | ard. (115). 


Be 


F 


(es te ea ; 
Hg (..-Qu--.) = min, ee | | V ps |(?de + | V(r, q)| pe 


Pa 
= mn, i V Ws 
Baleeuhesy 


Dieses Funktional mu zum Minimum gemacht werden durch Variation: 
von ws bei fester Polarisation, d.h. bei festen q,. Man kann zeigen [PEKARi 
(1951, § 16)], daB, wenn ,,die Bindungsenergie“ des Polarons bedeutend gréRer'; 


ist als die mittlere Energie der thermischen Schwingungen, genauer, wenn 


1,8 | J, | S3hep (7 + et 116) 
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ie Polarisation des Kristalls im Gebiet des Polarons nur schwach um die 
aeichgewichtspolarisation $8) (t —§) variiert, und daB die Wellenfunktion 
Elektrons ebenfalls nur wenig von wy (t —§) abweicht. Man kann also 
ie Polarisation in der Form 


Bt) = Bo (&@—s) + BP (r), (117) 


arstellen, wobei V’ (r) eine kleine Stérung ist, falls man § geeignet wihlt. 
nter Benutzung der iiblichen Stérungsrechnung kann man, wenn man 
» (t —§) als Lésung nullter Niherung auffaBt, die Energie des Elektrons 
nerster Niherung berechnen. Fiir den Grundzustand des Elektrons ergibt sich 


h? Z 
PES Cae a | =3,/ 1 V wo [2 dt — [x (t) Do (ve —§) dt. (118) 


it Do (te — §) =D [r, Wo (t — §)]. 

der Folge baie es, sich auf diese erste Niiherung der Stérungsrechnung 
u beschrinken. 

Als giinstigste Wahl von § ist diejenige zu betrachten, bei der H(...q,.. .) 
m genauesten durch Formel (118) gegeben wird. Nach dem Variations- 
rinzip (115) garantiert diejenige Wahl von § die gréBte Genauigkeit bei der 
Berechnung von £,(...q,...), bei der die Energie (118) ein Minimum ist. 
Die Minimumsbedingung hat die Form 


et...d,.--6 e AD) (x—F§ i 
0 ( aE | (r) 02 dr =0 ( = 1,2, 3). (119) 


Aus diesen drei Gleichungen ergeben sich die Komponenten des Vektors §, 
namlich &,, &, und £, in Abhangigkeit von $8 (r). § ist also ein Funktional 
von J (r), wie auch zu erwarten war. Mit anderen Worten: § ist eine Funktion 
der q,. 

Die oben angegebene Wahl von § kann man folgendermaBen deuten. Wir 
nehmen an, in einem triage polarisierbaren Medium mit der spezifischen 
Polarisation 38 (t,t) befinde sich eine freie masselose Ladung, die mit der 
Dichte e [wo (t — op verteilt ist ; die elektrostatische Induktion dieser Ladung 
betriigt Dy (tr —§) (§ ist die Koordinate des Mittelpunktes der Ladung). 
Behandelt man die Bewegung dieser Ladung auf Grund der klassischen 
Mechanik, so ergibt sich ihre Lage in jedem Zeitpunkt aus der Bedingung des 
Minimums der potentiellen Wechselwirkungsenergie dieser masselosen Ladung 
mit dem polarisierten Medium: 


—/(B (x, t) Do (c —§) dv = min. (120) 


o 


Diese Bedingung ist identisch mit (119). 

Wir gehen nun zur zweiten Etappe der Losung unseres Problems in der 
adiabatischen Naherung, zur Behandlung der Ionenbewegung tiber. Die 
HAMILTONfunktion des Systems, die die Ionenbewegung in Anwesenheit 
des Elektrons in der adiabatischen Niiherung bestimmt, ergibt sich durch 
Hinzufiigung des Gliedes HZ, (...q,...) zur HAMILTONfunktion der Ionen- 
eigenschwingungen; dieses Glied stellt die zusitzliche potentielle Energie 


408 S. I. PEKAR 


der Ionen dar (siehe Formel (15)). Man kann also die HAMILTONfunktion 
der Ionen hinschreiben : 


A if: 2 
Sie hor & @ 7 See wal (121) ) 
wobei | 
qe oes Bas (122) 
ho Co 


Te] (v) dr. (123) 


ist. Zu beachten ist, daB § eine Funktion der q, ist. Diese Funktion wird 
implizit durch die Gleichung (119) gegeben. 
Jetzt wenden wir uns der Lésung der Wellengleichung mit der HAMILTON-- 
funktion (121) zu. An Stelle der alten Variablen q, fiihren wir neue ein: 


} 


poet (124a) | 


Gey = V% ed x, (gees Ee (124 b)| 
(eS ae Ue Pa as ie Tae) eC (124¢) 


Die neuen (gestrichenen) Variablen sind etwas unsymmetrisch eingefiihrt} 
worden: von der beliebig groBen Gesamtzahl » der Koordinaten werden »y — 
durch die Formeln (124a), drei durch die Formeln (124b) bestimmt, wahrend} 
als letzte drei neue Koordinaten die drei Komponenten des Vektors § ge- 
wihlt werden. Die physikalischen Ergebnisse hingen nicht davon ab, welche 
sechs alten Koordinaten in der Transformation (124) gesondert in Form de 
Koordinaten q+,, (J = 1, 2, 3) auftreten. 


In den neuen Koordinaten geschrieben, hat die HAmMILTONfunktion (121) 
folgende Form [PEKAR (1951)]: 


pie pe a OW pete G i Bees Lia, (9 
ee rey O x ii 79. O-ar —- 5 ~ 
he : 2 aR , 2M ;2, 08 

dabei ist 


=, (126: 


A 2 2 
nes. y Odnt Ode OP (127 
aM jo. jm OG} Oe; 00, 08 | 


4S” 
ae 


x 


bezeichnet die Summierung iiber siimtliche Werte von x mit Ausnahmo 


von dreien, nimlich — x, (J = 1, 2, 3). 
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Die Glieder H’ und H” werden als kleine Stérungen aufgefaBt und in nullter 
Naherung weggelassen. Hierdurch nehmen die Eigenfunktionen des Operators 
(125) folgende Form an: 


ign , 
ve = Lh eité jek Pn, (Qn) : (128) 


und die Energieeigenwerte der Ionen und damit auch die des Gesamtsystems 
[PAULI (1933)] sind: 


me Divo ke 
lg Sle — Jo =i ha > Ny -- 3 ate 9M’ (129) 


Die Funktion ®,, wird durch Formel (22) bestimmt. Die Konstante M stellt 
die effektive Masse des Polarons bei seiner Wanderung dar: 


c r [ODo\? pe2 \8 ec 
M=—*,| (= Ic 2 5,8 -10-3(F— sh 1 
iene es dt {5,8 - 10 7B Co oe (130) 


M muB als triage Masse des Polarons aufgefaBt werden, die z. B. in Erscheinung 
tritt, wenn man das Polaron durch ein diuferes elektrisches Feld zu beschleu- 
nigen versucht. Die Formeln (128) und (129) zeigen, dafi der Mittelpunkt § 
des Polarons sich wie ein freies Teilchen mit der Masse M bewegt. Die effektive 
Masse des Polarons kann wesentlich von der Masse des freien Elektrons sowie 
auch von der effektiven Masse uw des Elektrons verschieden sein, die in dem 
betreffenden Kristall die trige Masse des Leitungselektrons in dem Fall 
unbeweglich in den Gitterpunkten sitzender Ionen darstellt. In Tabelle V 
sind die effektiven Massen des Polarons in einigen Kristallen angegeben. 


Tabelle V 


Die Funktion (128) stellt die Welle des Polaronenzustandes dar, die Formel 
(129) die Abhingigkeit der Energie dieser Welle vom Wellenvektor f. Kon- 
struiert man aus diesen Wellen ein Wellenpaket, so wird dessen Gruppen- 
geschwindigkeit durch die tibliche Formel 


] ht 

= =—, 131 
Deets Vt Ayn, i (tol) 
gegeben. Der klassischen Geschwindigkeit der Polarons entspricht also die 
| Gruppengeschwindigkeit » des Wellenpaketes. Mit Hilfe der Forme] (131) kann 


man das letzte Glied in (129) folgendermaBen schreiben: 
a oA 


29 
2M § 2 (132) 
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Diese GriBe stellt die kinetische Energie der Wanderung des Polarons dar. 
Man kann die Formel (131) als DE BRoGLIEsche Beziehung fiir das Polaron| 
auffassen; Mv spielt hierbei die Rolle des Polaronenimpulses. 
Die Striche an den Indizes des Produktes und der Summe (128) und (129) 
bedeuten, daS das betrachtete System drei Schwingungsfreiheitsgrade wenige 
hat als ein Dielektrikum ohne Leitungselektron. Man kann das leicht mit! 
Hilfe folgender Uberlegung verstehen. Die Hinfiihrung eines Elektrons, a a 
Zustand der Ionenbewegung adiabatisch folgt, darf die Gesamtzahl de 
Ionenfreiheitsgrade nicht verindern. Vor der Einfiihrung des Leitungs- 
elektrons waren alle » Freiheitsgrade der Ionen Schwingungsfreiheitsgrade, 
Nach Einfiihrung des Elektrons und Bildung des Polarons treten drei Frei- 
heitsgrade &,, £,, und &, der Ionen auf, deren Veraénderung keinen ee 
sondern Translationscharakter triigt. Diese drei neuen Freiheitsgrade kénnen 
nur durch Verringerung der Anzahl der Schwingungsfreiheitsgrade um drei 
auftreten. 
Die Energie des Polaronengrundzustandes ergibt sich aus der Formel (129). 
wenn man in dieser n = 0, f = 0 setzt. Sie ist in der Tat etwas (um 3/2 ha) 
kleiner als die Energie (21), die mit Hilfe der direkten Variationsmethodd 
mit dem Produktansatz (25) als Naiherungsfunktion gewonnen wurde. Das 
war auch zu erwarten, da die direkte Variationsmethode mit der Naherungs- 
funktion (25), wie bereits erwihnt, zu einem etwas zu groBen Wert fir did 
Energie des Systems fiihren muB. 
Die oben durchgefiihrte Rechnung bezieht sich auf den Fall, daB sich daa 
Elektron in dem Polarisationspotentialtopf im Grundzustand mit der Energie E, 
befindet. Kine véllig entsprechende Rechnung 1liBt sich aber auch dureh~ 
fiihren, wenn sich das Elektron im Polarisationspotentialtopf auf einem an- 
geregten Term , befindet. ) 
Der Kiirze halber gehen wir hier nicht auf die Kriterien fiir die Anwendbar 
keit der Theorie ein. Der Leser findet sie bei PEKAR (1951, §$ 7 und 16). ; 
Die oben gewonnenen Polaronenzustiinde des Leitungselektrons im Ideal. 
kristall unterscheiden sich grundsitzlich von den Ergebnissen der Binders 
theorie, in der das Elektron eine nichtabklingende Wellenfunktion besitzy 
und in der keinerlei Korrelation zwischen der Ionen- und der Elektronen: 
bewegung besteht. Es muB hervorgehoben werden, daB gewisse qualitative 
Ziige der oben dargelegten Polaronentheorie schon friiher von sowjetischer 
Theoretikern in Form gelegentlicher Hinweise angedeutet wurden. So sprach 
L. D. LANDAU (1933) den wichtigen Gedanken der Autolokalisation des 
Klektrons im Idealkristall infolge der Deformation des Gitters durch das 
Feld des Elektrons selbst aus. Diese lokalen Zustinde wurden jedoch al: 
unbeweglich angesehen, und L. D. LANDAU versuchte sie mit den F-Zentrer 
in den Alkalihalogenidkristallen zu identifizieren. Im Jahre 1936 bemerkte 
J. I. FRENKEL (1936), daB das Leitungselektron die benachbarten Atome 
des Kristalls deformieren und daf diese lokale Deformation sich mit dem 
Elektron durch den Kristall bewegen muB. Im Jahre 1936 versuchte D. I 
BLOCHINZEW (1936) auf Grund der Niherung der stark gebundenen Eleki 
tronen herauszubekommen, in welchen Kristallen man die friither von LANDAU 
angedeutete Autolokalisation der Elektronen zu erwarten hat. Damals was 
der Weg zur quantitativen Behandlung der autolokalen Zustiinde noch nich! 
| 
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yefunden, und deshalb gelang es noch nicht, ihre Existenz zu beweisen und 
ihre Higenschaften zu untersuchen. Die erwiihnten Arbeiten hatten jedoch 
nen bestimmenden Einflu8 auf die Arbeiten des Verfassers auf diesem 
Gebiet. 
egen der groBen Wichtigkeit des Polaronenproblems fiir die Grundlagen 
er Theorie der elektrischen, photoelektrischen und optischen Erscheinungen 
Halbleitern und Isolatoren leitete der Verfasser die Ergebnisse dieses 
Paragraphen unabhingig mit Hilfe von vier verschiedenen mathematischen 
ethoden ab, wobei die niherungsweise Abschitzung mit Hilfe der direkten 
Variationsmethode, die sich auf den Produktansatz (25) stiitzt, nicht mit- 
gerechnet ist. Alle Rechnungen fiihrten zu den gleichen Ergebnissen. Die erste 
Methode war halbklassisch: die Elektronenbewegung wurde auf Grund der 
Wellenmechanik, die Ionenbewegung auf Grund der klassischen Mechanik 
behandelt [PEKAR (1946b, c; 1947, 1948b)]. Der Vorzug dieser Methode 
besteht in ihrer Anschaulichkeit. Gemeinsam mit L. D. LANDAU berechnete 
der Verfasser mit Hilfe dieser Methode die effektive Masse M des Polarons 
{LANDAU und PEKAR (1948)]. Die anderen Methoden, bei denen auch die 
Tfonenbewegung quantenmechanisch behandelt wird, findet der Leser bei 
PEKAR (1949, 1951). Wir bemerken, daf unsere Ergebnisse auch von N. N. 
BoGoLjuBow (1950) und S. W. TyaBLikow (195la) gewonnen wurden, die 
ihre eigene mathematische Methode benutzten!). 


9. Die Polaronenleitfahigkeit der Kristalle 


Mit Hilfe der gewdhnlichen Statistik kann man zeigen, daB im Fall des 
thermischen Gleichgewichts die Polaronen eine MAXWELLverteilung ihrer 
Gruppenschwindigkeiten aufweisen [PEKAR (1951, § 19)}. Die im vorigen 
Paragraphen behandelten Wellen der Polaronenzustainde sind nicht die 
exakten, sondern nur angeniherte stationiire Zustiinde des Systems, da sie 


durch Vernachlassigung der Glieder H’ und H” der HaAmittTonfunktion 
herauskamen. Die Beriicksichtigung dieser Glieder als kleine St6rungen in 
erster Naiherung fiihrt zu einer Streuung dieser Wellen an den thermischen 
Tonenschwingungen. Die Streuung der Polaronen an den optischen Ionen- 
schwingungen wird eingehend bei PEKAR (1951) behandelt. Beriicksichtigt 
man diese Streuung, so kann man mit Hilfe der iiblichen kinetischen Glei- 
chung die entsprechende Beweglichkeit des Polarons im aiuBeren elektrischen 
Feld berechnen. Die Rechnung zeigt, daB bei tiefen Temperaturen gewohnlich 


die Streuung, die auf der Stérung H’ beruht, iiberwiegt. Mit wachsender 


Temperatur jedoch steigt auch die relative Rolle der auf dem Glied H" 
beruhenden Streuung. Bei hinreichend hohen Temperaturen haben 
beide Streuungstypen gréBenordnungsmibig vergleichbare Wahrscheinlich- 
keiten. 


a) Anm. d. dtsch. Red.: MARKHAM und SeE1Tz, Phys. Rey. 74, 1014, 1948, haben die 
Bindungsenergie eines selbstlokalisierten Elektrons in NaCl berechnet und fanden mit 
self-consistent Lésung als thermische Dissoziationsenergie 0,13 eV. 


412 S. 1. PEKAR 


Weiter unten werden die Formeln fiir die Beweglichkeit 6’, die auf der Storung 
A’ berubt, und die Beweglichkeit 0’, die auf der Stérung H” beruht, an- 
gegeben. Bei tiefen Temperaturen, d. h. bei | 


kT < ha, (133) 


ist ho 
a 8 vn 
ae at h sab (Pee : (i = A (134) 
MM: Coe V2 Wo (: et 4h °) y 
7 


hw, 
perce wie 288 (kT)'! eS —2+4-- (135) 


8 = (OKONIKS: an (2 aM) W2 We Cy 7 


Bei hohen Temperaturen, d. h. bei Giiltigkeit der zu (133) entgegengesetzten: 
Ungleichung ist | 


Cy jek 
Oh == O88. See 136) 
ia ee ( 
9 2 
1" 288 h (136'y 


a eco (2a MET) 


Die tatsichliche Beweglichkeit des Polarons wird durch die kleinere der 
GroBen b' und b” bestimmt. Um die GréBe der Beweglichkeiten und ihrer 
Zusammenhang deutlicher zu machen, betrachten wir ein Zahlenbeispiet 
mit typischen Mittelwerten fiir die Parameter des Kristalls (siehe Tabelle V1) 
Die in der Tabelle angegebenen berechneten Werte fiir die Beweglichkeit des 
Polarons stimmen gr6BenordnungsgemiB mit den experimentel]l an Hand des 
elektrischen Leitfaihigkeit und der HAttkonstante bestimmten Werter 
iiberein. 
Tabelle VI 
&€=5, n= 1,58, w= mM, @ = 5-104 sec—1, hay = 0,083 eV, c = 0,2 


Tiefe Temperaturen 


; Hohe Temperaturen 
k —— 3 hw kT = 3 hoy 


b’ em?/sec + \ 


b’” em?/see + V 


Die Streuung der Polaronen durch die akustischen Tonenschwingunge¢ 
und die entsprechende Beweglichkeit wurden von J. J. PERLIN (1951) bo 
rechnet. Ks zeigte sich, da& bei Zimmertemperatur die Streuung an de¢ 
Polarisations- und nicht die an den akustischen Schwingungen die Haup/ 
rolle spielt. Bei sehr tiefen Temperaturen jedoch kann die Streuung an de 
vkustischen Schwingungen die Oberhand gewinnen. Bei den meisten Kristalle{ 


Ee 
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ritt dies bei so tiefen Temperaturen ein, daf dann die Streuung der Polaronen 
lurch die thermischen Schwingungen iiberhaupt kaum noch vorhanden ist 
ind daf die Beweglichkeit durch die Streuung an Fremdatomen und Fehl- 
stellen, an der Anharmonizitit der escape “inate es usw. bestimmt wird 
PEKAR (1948b, 1951)}. 

eben den oben behandelten Polaronen-Zustinden, in Tosen sich das Elektron 
m Grundzustand mit der Energie Ey befindet, gibt es auch Polaronen, in 
tenen sich das Elektron auf einem angeregten Term H, befindet. Ferner kann 
nfolge der thermischen Schwankung der Ionenschwingungen die lokale 
dlarisation des Kristalls verschwinden oder sogar ihr Vorzeichen wechseln. 
ierbei verschwindet der Polarisationspotentialtopf des Elektrons, und dieses 
ird zu einem Bandelektron im tiblichen Sinne. Diesen Vorgang bezeichnet 
man als thermische Dissoziation des Polarons. 

Wie die Polaronen im Grundzustand und in angeregten Zustaénden sind auch 
die Bandelektronen Ladungstraiger. Damit erhebt sich die Frage, welcher 
dieser Ladungstrager den Hauptanteil an der elektrischen Leitung hat. Die 
Konzentration der Polaronen mit einem Elektron im Grundzustand ist sehr 
viel gréBer als die Konzentrationen der angeregten Polaronen und der Band- 
elektronen. Das Konzentrationsverhaltnis der unangeregten Elektronen und 
der Bandelektronen 1a8t sich in grober Naiherung durch folgende einfache 
Forme! ausdriicken, die aus den Grundlagen der Gleichgewichtsstatistik 
folgt [PEKAR (1948b, 1951)]: 


ah ES (=) eWp/kT (137) 
NBand Lt 


dabei ist W, die thermische Dissoziationsenergie des Polarons; sie hat den 
VV 
ert 


W, = 0.0544 = pe bey. (138) 


GréBenordnungsmaBig ist Wy gleich einigen Zehntel eV. Die Polaronenzustiinde 
sind also energetisch vorteilhafter als die Bandzustiinde, d.h. die dissoziierten 
Zustinde des Elektrons. Das Verhiltnis Npoi/Npana schwankt bei Zimmer- 
temperatur fiir die meisten Kristalle zwischen 10° und 10°. Bei tiefen Tem- 
peraturen ist dieses Verhaltnis noch weitaus gréBer. Die Bandelektronen 
kénnten also nur dann einen wesentlichen Anteil zur elektrischen Leitung 
stellen, wenn ihre Beweglichkeit mindestens um den angegebenen Faktor 
jgroBer wire als die Polaronenbeweglichkeit. Damit also die Bandelektronen 
in der elektrischen Leitung mit den Polaronen konkurrieren kénnten, miibte 
man ihnen eine Beweglichkeit von unzulissiger GroBe zuschreiben, bei der 
‘schon in einem nicht allzu groBen Feld ihre Geschwindigkeit hdher als die 
'des Lichts wiirde. Diese hohe Beweglichkeit ist um GréBenordnungen hoher 
jals die experimentell bestimmte Beweglichkeit der Ladungstriger. Schreibt 
man den Bandelektronen eine Beweglichkeit von verniinftiger GréBen- 
}ordnung zu, so miissen die Polaronen den Hauptanteil an der elektrischen 
/ Leitung stellen, da ihre Konzentration weitaus tiberwiegt [PEKAR (1948b, 
1 1951)). Aus demselben Grunde spielen von siimtlichen Polaronenzustiinden 


b 
30 
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die unangeregten die Hauptrolle, falls das Kriterium der Anwendbarkei: 
der Theorie (116) erfiillt ist. Wie bereits erwahnt, stimmt die berechnet 
Beweglichkeit der unangeregten Polaronen gréBenordnungsmaBig mit de 
experimentell bestimmten Beweglichkeiten der Ladungstrager tiberein’). 
In Abb. 9 wird die Temperaturabhingigkeit der Beweglichkeit mit dem Exp 
ment verglichen2). Die ausgezogene Kurve stellt die Formel (134) fiir ha 
— 0,031 eV dar. Die Kreuze entsprechen den experimentellen Daten vor 
ENGELHARD (1933) fiir Cu,O. Bei 150° K wurde die theoretische Kurv 
mit Hilfe eines Faktors mit dem entsprechenden experimentellen Punkt 
Deckung gebracht. 
Die SchluBfolgerung, daB die Polaronen die wichtigsten Ladungstrager sind 
findet ihre Bestitigung auch in verschiedenen Anwendungen der Theorid 
cm? So wurden beispielsweise in § 7 die be 
Senet rechneten Werte fiir die thermische Dissc 
ziationsenergie der F-Zentren angegeber 
die mit dem Experiment iibereinstimme; 
(siehe Tabelle II). Diese Energien wurde: 
unter der Voraussetzung berechnet, da. 
sich durch Dissoziation eines #’-Zentrun 
ein Polaron und kein Bandelektron bilde 
Zur Bildung eines Bandelektrons wa 
eine héhere Energie erforderlich. 
Ferner zeigt die Untersuchung, daB jed 
lokalisierte Zustand des Elektrons i 
Kristall, der zu einem Fremdatom od: 
einer Gitterfehlstelle gehért, unendlid 
viele diskrete Elektronenterme besitz 
muB8, wobei diese Terme bei Anniherw 


0 ; an die Energie #, = 0 (an den unter 
He) PPE Ee CEP Se = . (a oe ae 
Tae Rand des Leitungsbandes) unendlich die 

Abb. 9: Temperaturabhingigkeit der Be- : a . 
wegiohwelt; der’ badhingekchper in aneinander riicken. Diese Behauptu 


Cu,0. Die ausgezogene Kurve wird verstaindlich, wenn man _ beden! 

stellt die Formel b = 27 (e363/T_1) : af 
dar, die Kreuze die experimen- 28 jedem lokalisierten Elektronenzustan 
tellen Daten von ENGELHARD. — ein Polarisationspotentialtopf entsprich 
der in groBen Abstinden COULOMBs¢ 
Form hat und infolgedessen ein diskretes Energiespektrum besitzt, dess 
oberer Teil dem Spektrum des Elektrons im Wasserstoffatom ihnelt. Bei d 
Absorption von Licht durch ein solches lokalisiertes Elektron ist am wah 
scheinlichsten der optische Ubergang des Elektrons auf den tiefsten p-Tern 
Ubergiinge unmittelbar ins kontinuierliche Spektrum sind sehr unwalk 
scheinlich. Wiren die Bandelektronen die Ladungstriger, so wiirde ei 
Absorption im Gebiet des Maximums des Absorptionskoeffizienten nic 
unmittelbar zu einer Photoleitfihigkeit fithren, sondern fiir deren Entstehw 


*) Anm. d. dtsch. Red.: Neuerdings ist es gelungen, die Elektronenbeweglichkeit 
NaCl zu bestimmen. REDFIELD, Phys. Rey. 91, 753, 1953, gibt als vorliufige Wé 
250 + 50 cm?/Vsec bei 82° K und etwa 40 cm2/V sec bei 200° K. | 
*) Anm. d. dtsch. Red.: Neueste Bestimmung: Nrexr, Ann. d. Physik 12, 297, 1953. | 
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wire noch eine thermische Schwankung notwendig, die das Elektron vom 
diskreten p-Term ins Leitfihigkeitsband bringt. In diesem Fall miiBte die 
Quantenausbeute an Photoelektronen stets wesentlich von der Temperatur 
abhingen und bei starker Abkiihlung des Kristalls stets gegen Null gehen. 
Das Experiment zeigt jedoch, daB es Fille gibt, wo die Quantenausbeute 
an Photoelektronen bei 7’ —+ 0 nicht gegen Null, sondern gegen Eins geht 
[PEKAR (1951)}. In den Fallen, in denen die Quantenausbeute bei starker 
Abkiihlung des Kristalls exponentiell gegen Null geht, ist die entsprechende 
thermische Aktivierungsenergie kleiner als diejenige, die man unter der 
Voraussetzung zu erwarten hitte, daB das Elektron ins Leitfaihigkeitsband 
gehoben wird. Die erwihnten Tatsachen sprechen dafiir, da8 zur Dissoziation 
eines Fremdzentrums und zur Bildung eines freien Ladungstrigers eine 
kleinere Energie erforderlich ist als zur Uberfiihrung eines Elektrons aus 
einem Fremdzentrum in das Leitfahigkeitsband. Hierdurch wird die Annahme 
bestitigt, da nicht die Bandelektronen, sondern die energieiirmeren Polaronen 
die Ladungstriager bilden. Eingehender wird die Quantenausbeute der Photo- 
dissoziation der F- und #’-Zentren bei PEKAR (1951) behandelt. 

Aus dem Gesagten geht hervor, dafB der Begriff des Polarons nicht nur einer 
der Grundbegriffe der Theorie der elektrischen Leitung ist, sondern da8 er 
auch bei der Behandlung vieler anderer Erscheinungen eine groBe Rolle spielt. 


10. Fall schwacher Wechselwirkung des Leitungselektrons 
mit den lonenschwingungen 


Nach der Veréffentlichung der wichtigsten Arbeiten iiber die Polaronentheorie 
die auf der adiabatischen Niherung beruhten, begann man, die Polaronen 
mit Hilfe anderer mathematischer Methoden zu untersuchen. Wir beschif- 
tigen uns mit einer dieser Methoden eingehender. Sie besteht darin, daB man 
in der HAmiLTONfunktion (102) das zweite Glied — die Wechselwirkung des 
Elektrons mit den Schwingungen — als klein ansieht und die tibliche Stérungs- 
theorie heranzieht. Die nullte und die erste Naherung wurden schon vor 
langerer Zeit behandelt und bilden bekanntlich den Kern der gewoéhnlichen 
, Bindertheorie“ der elektrischen Leitung. Diese Naherungen beriicksich- 
tigen die Polarisation des Gitters durch das Feld des Elektrons und die Riick- 
wirkung dieser Polarisation auf das Elektron noch nicht vollstindig. Zur 
Beriicksichtigung dieser Polarisation mu man die zweite Naherung der 
Stérungstheorie betrachten, was in den letzten Jahren auch durchgefiihrt 
wurde {FROHLICH, PELZER und ZIENAU (1950); TJABLIKOW (1951, 1952)]. 
Es zeigte sich, daB die Korrektur zweiter Ordnung an der Energie des Systems 
nur im Grenzfall tiefer Temperaturen konvergiert (endlich bleibt), ferner 
nur dann, wenn die kinetische Energie des Leitungselektrons kleiner als ha, 
ist. Diese Korrektur an der Energie ist negativ; eine Beriicksichtigung der 
Polarisation fiihrt also zu einer Erniedrigung der Energieterme des Systems. 
Die Rechnung enthiillte auch noch andere Ziige der Polaronentheorie. So 
bewegt sich z. B., wie S. W. TJABLIKOW (1951, 1952) zeigte, das Leitungs- 
elektron in Begleitung der lokalen triigen Polarisation des Kristalls, und die 
effektive Masse M des Ladungstrigers ist infolgedessen gréBer als die effektive 
28 Zeitschrift ,,Fortschritte der Physik‘ 
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Masse des Elektrons in einem Kristall mit unbeweglichen Ionen. Bis auf 
Glieder héherer als zweiter Ordnung gilt 


eam tee 
ety peta) f (139) | 
lu 12 h 7 ha 
Man mu8 jedoch hervorheben, daB die quantitativen Ergebnisse der erwahnten 
Stérungsmethode wesentlich verschieden von den aut Grund der adiabatischen 


Naherung gewonnenen SchluBfolgerungen von § 8 sind. Das ist auch nicht) 
iiberraschend, denn die Stérungsrechnung ist anwendbar, wenn | 


Coe Tid | 
ues 1 | 
h D a | . 
d. h. ‘ (140)| 
lol 2 (ie | | 
ha 


ist, wihrend die adiabatische Naherung gerade bei Giiltigkeit der umgekehrten! 
Ungleichung verwendbar ist [s. Formel (116), in der mj) = 0 zu setzen ist}.| 
Die adiabatische Naiherung und die Stdrungsrechnung haben also keinent 
gemeinsamen Anwendungsbereich und sind bei Giiltigkeit zweier entgegen-i 
gesetzter Grenzfille giltig, namlich der Fille starker bzw. schwacher Wechsel- 
wirkung des Elektrons mit den lonenschwingungen. 

Bei der Untersuchung einer Reihe von Alkalihalogenidkristallen und dee 
Cu,0 [PEKAR (1951, Anhang)| zeigte es sich, daB stets die zu (140) entgegena 
gesetzte Ungleichung erfiilit war. In diesen Kristallen ist also die adiabatische 
Naherung giiltig und die St6rungsrechnung unbrauchbar. Kristalle miti 
einem Cy, das so klein ist, dafgS das Kriterium (140) erfiillt ist, wiirden in ihrem 
Eigenschaften den homdopolaren Kristallen sehr nahekommen, und dié 
Wechselwirkung des Elektrons mit den Polarisationsschwingungen wiirde im 
ihnen offenbar eine untergeordnete Rolle spielen). 


11. Kritik des Beweises der Selbstlokalisation des Elektrons 
von Morr und Gurnrty 


Mott und GURNEY (1940, S. 8Off.) gaben einen Beweis fiir das Vorhanden : 
sein von ,,self consistent’ selbstlokalisierten Elektronenzustinden im Ideal | 


*) Anm.-d. dtsch. Red.: Die Theorie des Polarons wurde neuerdings behandelt von GURARI 
Phil. Mag. 44, 329, 1953 (Weiterfiihrung der Methode von FROHLICH, PELZER una 
ZIENAU |, c.) sowie von LEE und Pines, Phys. Rev. 92, 883, 1953 (dort Zitate friihere: 
Arbeiten; Ubertragung der von Tomonaca fiir Mesonen entwickelten ,,intermediat« 
coupling method*‘). Diese Autoren sind der Ansicht, da ihre Methode besser sei als dili 
von PEKAR und begriinden dies u. a. damit, daB sie eine merklich tiefere Energie be 
der Anwendung des Variationsprinzips erhalten (— 0,16 eV gegentiber — 0,09 eV). Sil 
erhalten eine Polaronenmasse, die nur wenig gréBer ist als die Masse des Elektron: 
(z. B. fiir NaCl 1,9m,); als PeKars Wert wird 17 m, angegeben. Dies stimmt mij 
Tab. V erst tiberein, wenn man den Faktor 3 in (130) beriicksichtigt und mw aus Tab. I 
entnimmt. | 
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kristall, wie sie von LANDAU (1933) vorhergesagt worden waren. Ihr Beweis 
verliuft in unseren Bezeichnungen folgendermaBen. Die trige Polarisation 
des Mediums wird durch das Feld der y-Wolke des lokalisierten Elektrons 
nach den Formeln (104) und (77) bestimmt. Der entsprechende Polarisations- 
potentialtopf des Elektrons hat die Form: 


A re 
V (t) = — Ge | Yo (©) | at « (141) 
kes 
Im Grenzfall groBer r hat dieses Potential CoULOMBsche Form: 
2 
r—> oo Vjyr— =. (142) 


Ferner gibt es nach einem bekannten Satz der Quantenmechanik, falls das 
Potential die asymptotische Form (142) hat, unendlich viele diskrete Energie- 
terme des Elektrons, wobei fiir Zustinde mit hinreichend groBem effektiven 
Radius das Energiespektrum des Elektrons in das Spektrum eines wasser- 
stoffahnlichen Atoms tibergeht. In dem Polarisationspotentialtopf (141), 
den das Feld des Elektrons erzeugt, gibt es also diskrete Terme und folglich 
auch lokale Elektronenzustinde. 
Zu diesem Beweis méchten wir folgendes bemerken: ry sei der effektive Radius 
der wo-Funktion des Elektrons, die in dem Ausdruck (141) vorkommt. Wie 
man leicht sieht, hat dann das Potential (141) nur im Gebiet r > 7) COULOMB- 
sche Form. Der oben erwaihnte quantenmechanische Satz gewiihrleistet die 
Existenz lokaler Zustainde nur fiir hinreichend groke effektive Radien 7,, 
genauer fiir Zustainde, in denen der gréBere Teil der yw-Wolke des Elektrons 
im Gebiet des CouLOMBpotentials gelegen ist; bei diesen Zustinden muB 
T,>7, sein. Der Satz gewiahrleistet durchaus nicht die Existenz eines 
stationiren Zustandes y, mit einem Radius 7) in dem Potentialtopf (141), 
da in diesem Fall die wy -Wolke des Elektrons in ein Gebiet eingeschlossen 
ware, in dem das Potential wesentlich verschieden vom COULOMB- 
potential ist. 
Wie groB auch der Radius 7, der Funktion-y, in (141) sei, der Satz gewahr- 
leistet nur die Existenz lokaler Zustiinde mit weitaus groBeren Radien r, in 
dem Potentialtopf. Die Uberlegungen von Mott und GuRNEY bilden also 
- keinen Beweis fiir die Existenz von ,,self consistent®’ Zustanden. 
_Nimmt man an, da8 infolge einer Fluktuation der Schwingungen oder einer 
zeitweiligen St6rung durch duBere Krafte ein Polarisationspotentialtopf 
von der Form (141) entstanden ist, so kann sich nach dem erwahnten Satz 
in diesem ein Elektron lokalisieren, indem es in den Zustand yp, fallt. Da jedoch 
7, >7, ist, kann man leicht zeigen, daB die y,-Wolke des Elektrons zu sehr 
, verschmiert ist, daB sein Feld nicht ausreicht, um die urspriingliche Polari- 
sation des Kristalls aufrechtzuerhalten. Nach Aufhéren der auBeren Krifte 
beginnt also eine Depolarisation, durch die sich der Radius der Funktion 
y, vergréBert. Das Problem, ob diese Depolarisation zum volligen Verschwin- 
_den des Polarisationspotentialtopfes fiihrt, oder ob sich ein ,,self consistent‘ 
lokalisierter Zustand bildet, laBt sich auf Grund der Erwigungen von MoTT 
und GURNEY nicht lésen. 
28* 
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Ganz Abhnliches li8t sich zu dem Beweis der Behauptung sagen, daf dic 
Vereinigung eines zweiten Elektrons mit einem F-Zentrum und damit 4 
Bildung eines stabilen #’-Zentrums méglich sei [MoTT und GURNEY (1940 


S. 130)]. Ubersetzt von H. VoGEL| 
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1. Einleitung 


Die phinomenologische Theorie der Supraleitung ist von LONDON une 
v. LAUE begriindet worden. Es wird angenommen, dai sich der Strom im 
Supraleiter zusammenseizt aus dem ,,Normalstrom” §” und dem Suprastroni 
°. Es gelten die MAXWELLschen Gleichungen (s. unten, Abschnitt 3) einschlie&! 
lich des OHMschen Gesetzes fiir den Normalstrom. Fiir 3 gelten zusiitzlich dil 
beiden Gleichungen 


y fo" o 

eee) A 
AS =e (1.1) 
crotaS?+ 9 =0. (LE 


Dabei ist die LONDoNsche Konstante A eine fiir den Supraleiter charakteristisehl 
Materialkonstante. 
Von diesen beiden Gleichungen ist die erste als Ersatz fiir das OHMsche Geset: 
verhiltnismaBig evident, die zweite dagegen eine ad hoc eingefiihrte Zusate 
hypothese, die es aber nun gerade erméglicht, Experimente an einfach ung 
mehrfach zusammenhiingenden Supraleitern im Magnetfeld richtig zu be 
schreiben. Die Gleichung (1.2) hat den Charakter einer Nebenbedingung: Di 
Zeitableitung der linken Seite von (1.2) verschwindet als Folge von (1.1) ung 
den MAXWELIschen Gleichungen. Die LoNDoNschen Gleichungen geben (vow 
Skineffekt abgesehen, vgl. PIPPARD 1950) eine gute Beschreibung der meiste 
Experimente, wie das in den Monographien von v. LAUE (1949) und LoNDO} 
(1950) im einzelnen dargestellt ist. 
Trotzdem hat die LONDONsche Theorie in der letzten Zeit verschiedene Uni 
formungen erfahren, die durch modellmibige Vorstellungen nahegeleet wurder 
V. LAUE (1948) hat die Méglichkeit betrachtet, daB die LONDONsche Konstan’ 
ein Tensor ist. IKOPPE (1949, 1950) hat von der HEISENBERGschen Theorie au 
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gehend eine nichtlineare Verallgemeinerung entwickelt; zunichst unter etwas 
einschriinkenden Annahmen, die dann von vy. LAUVE (1949b) beseitigt wurden. 
Spater hat Koppe (1951) die ,,Konstante“ 4 selbst als Feldfunktion betrachtet. 
SehlieBlich hat GINSBURG (1950) eine Theorie vorgeschlagen, die zwar keine 
echte phinomenologische Theorie ist, da sie makroskopisch unbeobachtbare 
GréBen enthalt, die aber wenigstens in gewissen Niiherungen auf eine nicht- 
lineare Theorie mit veriinderlichem 4 hinauslauft. 

Die Schwierigkeit derartiger Verallgemeinerungen liegt darin, daB man die 
beiden Gleichungen (1.1), (1.2) nicht beliebig abindern kann, sondern immer 
zeigen mu, dafs das System nachher noch gewisse Forderungen erfiillt. Diese 
sind: 1. Der Schwellenwert des dicken Supraleiters muB unverindert bleiben. 
2. Der Zustand des einfach zusammenhiingenden Supraleiters im Magnetfeld 
mu durch das Feld allein bestimmt und unabhingig von der Vorgeschichte 
sein. Beide Fragen stehen in engem Zusammenhang mit der Thermodynamik 
der Supraleitung. Man kann sie am besten iibersehen, wenn man von vorn- 
herein von der Thermodynamik irreversibler Prozesse ausgeht. Das soll im 
folgenden geschehen, Dabei ergibt sich eine geschlossene Theorie, die alle bisher 
betrachteten als Sonderfalle enthalt. 


2. Thermodynamik irreversibler Prozesse') 


In einem Kontinuum werden thermodynamische Prozesse beschrieben durch 
einen Wirmestrom 9 und einen Arbeitsstrom YY, die mit der Energiedichte w, 
der Entropiedichte s und der Dichte der Entropie-Erzeugung o zu verkniipfen 
sind. Der erste Hauptsatz, der die Erhaltung der Energie ausspricht, lautet 
offenbar : 


O tite 
div % + div Q + a F—0), (2.1) 
Der zweite Hauptsatz hat die Form: 
ry) ds 
liv {— = =o 0% 29 
div (| se at = (2.2) 
Die wesentliche Aussage ist dabei, daB die Dichte der Entropievermehrung o 
positiv sein muf. Fiihrt man die Dichte der freien Energie f = u — T's ein, 
so bekommt man 
] OT 
= - = div Wf + a Si grad T + s AK Soli Neyt. (2.9) 
Fiir isotherme Prozesse ergibt sich daraus: 
cd ipa (2.4) 


Ot 
Diese Gleichung wird normalerweise zur Ermittlung der freien Energie benutzt. 
Wir wollen im folgenden den Zusammenhang umkehren, indem wir die freie 


Energie f vorgeben und dann feststellen, welche Bewegungsgleichungen auf 
Grund von (2.4) thermodynamisch méglich sind. 


1) Vel. etwa DE Groot (1951). 
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Zur Aufteilung der rechten Seite von (2.4) in eine Divergenz und ein positives | 
Glied ist noch zu bemerken, daB 2% und o Funktionen der Zustandsgr6éBen sein 
miissen. Ohne diese Einschrinkung kénnte man natiirlich immer Vektor- 
felder 9 finden, die (2.4) mit o = 0 erfiillen. So lassen sich aber nur solche 
Terme zu div 9 zusammenfassen, die mindestens linear in den Raumableitungen 
der ZustandsgréBen sind. Dieses Argument wird im folgenden wiederholt 
benutzt werden. | 
| 


3. Freie Energie und Feldgleidhungen des Supraleiters | 


Wir betrachten zunichst isotherme Zustiinde. Ein Zustand ist dann beschrie- | 
ben, wenn ©, §, der Normalstrom 3", der Suprastrom 3* und gewisse Para- 
meter 7;, die spiter eingefiihrt werden, als Funktionen von ¢ gegeben sind. 
Die freie Energie f mu8 demnach eine Funktion dieser Groen sein. 

v. LAUE hat neben dem Suprastrom Y* noch einen Supraimpuls & eingefiihrt. 
In der Mechanik spielt der Impuls eine grundlegendere Rolle als die Geschwin- 
digkeit. Wir benutzen diese Assoziation und betrachten im folgenden © als die 
wesentliche ZustandsgréBe des Supraleiters, die in der obigen Liste an die 
Stelle von $* zu treten hat. 

Um die wesentlichen Gesichtspunkte besser hervortreten zu lassen, machen wir 
im folgenden die iibliche Annahme, da Dielektrizitatskonstante und Perme- 
abilitiit des Supraleiters gleich 1 sind!). Dann folgt aus der MAXWELLschen 
Theorie, daB f die Form 


f=z (V+ 9) + 2G, m) (3.1) 


haben mu, wobei 2 nicht mehr von © und § abhingt. f ist nur bis auf eine © 
Konstante bestimmt. Es sei so definiert, daB f die Differenz zwischen der ab- — 
soluten freien Energie und der freien Energie /,, im feldfreien normalleitenden 
Zustand ist; die wahre freie Energie ist also f + f,. Diese Definition spielt — 
zunachst keine Rolle, wird sich aber spiiter als zweckmaBig erweisen. 

In der urspriinglichen Fassung der LoNDONschen Theorie ist z. B. 


Lined 2 
Ss I §*— A}, (3.2) | 


wobei Af = f, —f, den Phasenumwandlungsanteil der freien Energie be- | 
deutet. Es ist manchmal bequem, Q allgemein in dieser Weise zu trennen; wir ’ 


definieren : 2 (G, ni) = 2° (G, ni) —-Af (ni); Q° =0 fir G =0. (3.3) | 
Wir setzen nun (3.1) in (2.4) ein und formen zunichst den Feldanteil um: : 


oe gp 

dt OF ssuteum Oe Ot” 
*) Im Hinblick auf die Behauptung von GinspurG ( 1952), daB Supraleiter extrem hohe » 
Werte der Dielektrizitiitskonstanten aufwiesen, sei bemerkt, daB es keine Schwierigkeiten | 
bereitet, sich von dieser Annahme frei zu machen. Man mu dann neben € und § noch D und | 
% einfiihren. Neue Effekte (Abhingigkeit der DK vom Suprastrom) treten nur auf, falls; 
auch 2 von den elektrischen FeldgréBen abhiingt. | 
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etzt ziehen wir die MAXWELLschen Gleichungen heran: 


of = —crob € 

0€ ss 

A rte 
AO Gy C. 


uit ergibt sich aus (2.4) . 
Wolo s St ast 
div c[$ x ©] + og & + eg = divYA + To. 
ier sind die beiden ersten Terme links bereits eine Divergenz (des Poynting- 


ektors) und ein positiver Dissipationsterm (die JOULEsche Wirme). Wenn 
ich die beiden verbleibenden Terme ebenso umformen lassen: 


Exe = = div W + To! (3.4) 


ann ist (2.4) mit Y= A’ +c[H x E] und To= To’ + o,, © erfiillt. Wir 
nutzen in Zukunft (3.4), wo € und § schon weitgehend eliminiert sind. 
ir betrachten nun zunichst den Fall, daB 2 nur von & abhangt und alle 
nderen Parameter feste Werte haben. (Der allgemeine Fall wird in Abschnitt 6 
achgeholt.) Aus (3.4) wird dann 


02 06 
Waser eine 1 / 
ES aG at div YW -+ To’. 
abei ist 23 der Vektor mit den Komponenten (a ee 1) 


iese Gleichung muB zweierlei liefern: 1. die ,, Bewegungsgleichung“ fiir ®, und 

. eine Beziehung zwischen & und 9°, die bisher noch offen gelassen worden war. 
a keinerlei Raumableitungen von FeldgréBen auftreten, muB YF’ = 0 sein. 
s liegt nahe, dann auch o’ = 0 zu setzen, und 


6 =E (3.5) 
pi is 8 - 
alae ag (3.6) 


zu wahlen. 
immt man die Rotation von (3.5) und driickt rot © durch § aus, so ergibt 


ich 
0 1 
5 (rot & + 3) — 0 


Wir k6énnen also, ohne in Widerspriiche zu geraten, an Stelle von (1.2) postu- 


lieren: . 
crotG+§=0 |. (3.7) 


An dieser Gleichung, die die allgemeine Formulierung der zweiten LON DoNschen 
Gleichung darstellt, soll im folgenden unter allen Umstanden festgehalten 
werden. Wir hatten schon oben darauf hingewiesen, dai sie den Charakter 
einer Nebenbedingung hat. Aus diesem Grunde konnte sie bei der Ableitung 
der Bewegungsgleichungen nicht mit herauskommen. Es sei an dieser Stelle 
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nochmals betont, da es im Wesen der hier vorgetragenen Methode liegt, dak 
man von der Thermodynamik ausgeht und die Bewegungsgleichungen so zu 
bestimmen sucht, daB sie sich in die Thermodynamik einfiigen. Die letztere 
Forderung stellt eine notwendige Bedingung dar, legt aber die Bewegungs. 
gleichungen nicht eindeutig fest. 
Im Hinblick auf spateres sei angemerkt, daf (3.7) nicht nur mit (3.5), sonderr 
auch noch mit der erweiterten Fassung 

G = € + grad V (3.6 
vertriiglich ist, da die Rotation eines Gradienten verschwindet. | 
Vorlaufig kénnen wir jedoch davon keinen Gebrauch machen, und zwar a 
Grund des folgenden allgemeinen Prinzips: In einem Supraleiter sind statio 
naire Zustiinde méglich, bei denen alle Zeitableitungen verschwinden. AuBen 
dem miissen aber fiir stationare Zustinde alle Dissipationsglieder verschwinden 
und das so entstehende Gleichungssystem mu8 kompatibel sein. Nun haben wil 
als Dissipationsglied bis jetzt nur die JOULEsche Warme, und es folgt, da iv 
stationdaren Zustanden © = 0 sein muB. Dann folgt aber aus (3.8), daB V fiw 
stationire Zustinde konstant sein mu8. Das kann man nur erreichen, went 
man fiir V eine Funktion der Dissipationsglieder wihlt. Im Augenblick komm 
also nur in Frage: V = V(@) und das wiirde in (3.4b) links ein Glies 


— 4g grad V(?) ergeben, welches gleich Z7'o’ > 0, also nichtnegativ sei: 


miBte. Das ist aber nicht méglich, da © und € voneinander unabhangig sina 


4. Spezielle Ansatze 


Wir erhalten verschiedene Theorien der Supraleitung, wenn wir fir 2 (@ 
resp. 2° (G) geeignete Ansitze machen. Sie unterscheiden sich lediglich b 
ziiglich des Zusammenhangs zwischen Y* und & nach (3.6). 


: | 
A. Der bereits erwiihnte Ansatz Q° = FI 2 (4. 


. & : 2 
fiihrt auf ¥* = 7 Setzt man © = AY in (3.5) und (3.7) ein, so kommt ma 
auf die urspriingliche LoNDOoNsche Fassung zuriick. 


B. Eine naheliegende Verallgemeinerung von (4.1) ist 


| 
eo) uj re “4 
SP = Ss nie Gi Ge. . (4. 
ik 4 
Die Koeffizientenmatrix der quadratischen Form kann als symmetris: 
vorausgesetzt werden, da sich ein antisymmetrischer Teil wegheben wiirc 
Aus (4.2) folet 
5 oe ND ¥ 
vi = 2 Nir Ge 
k 

oder, wenn A;x die zu 7;% inverse Matrix ist, 
’ as 
GS; — > Ani Si- 
i | 
Das ist die von v. LAUVE (1948) vorgeschlagene Verallgemeinerung. Da 
Reziproke einer symmetrischen Matrix wieder symmetrisch ist, ergeben si 

A J ay Sur ah ; 4 . 

die y. LAvEschen Sy mmetriebedingungen A;; = A,; hier ganz von selbst. 
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C. Die Gleichungen (3.5) bis (3.7) enthalten bereits die allgemeine Form der 
nichtlinearen Theorie, solange man nur @ als Variable hat. Die Gestalt von Q 
ist nicht ganz gleichgiiltig. Wie unten gezeigt wird, muB die Bedingung 


dQ 02 
(G — @’) te Ga a 0 
erfiillt sein. Anschaulich bedeutet das, da die Zunahme des Suprastromes 
und die Zunahme des Supraimpulses immer einen spitzen Winkel einschlieBen 
miissen ; das heiBt also etwa, daB ein elektrisches Feld wirklich beschleunigend 
wirken muB. 
In den Arbeiten von KoppeE (1949, 1950) und v. LAUE (1949, 1949b)- wird noch 
der Suprastrom 3 als Variable benutzt. Auf Grund von (3.6) kommt das auf 
eine Legendre-Transformation hinaus, und man kann (3.6) ersetzen durch 


(4.3) 


dQ _ 
Aus (4.4) folgt dann unmittelbar die vy. LAUEsche Integrabilititsbedingung 
08; OG, 


OSt ABI 
Die Bedingung (4.3) kann geschrieben werden AQ*. A > 0 und ist deshalb 
gegeniiber der Legendre-Transformation invariant. 


5. Ein Beispiel zur nichtlinearen Theorie’) 


Die nichtlineare Theorie war urspriinglich entwickelt worden, um zu zeigen, 
daB sich die von der HEISENBERGschen Theorie geforderten Stromsattigungs- 
effekte widerspruchsfrei in die phinomenologische Theorie einbauen lassen. Die 
Forderung war dabei, daB der Betrag | $* | des Suprastromes einen temperatur- 
abhangigen Maximalwert Jy nirgends auch nur kurzzeitig tiberschreiten darf. 
Da nach (3.5) der Supraimpuls unbeschrankt zunehmen kann, mu8 man das 
durch eine geeignete Wah] von %* (@%) erreichen. Da es nicht méglich war, 
theoretische Aussagen iiber diese Funktion (resp. tiber 2°) zu bekommen, 
wurde der mathematisch bequeme Ansatz 

ease bi a 

AV1 + C/A? Dig 


gewahlt, der entsteht, wenn man setzt 
nee = ww ivan @ + G?/2? D5)! =a Aly. (552) 


Aus (5.1) folgt unmittelbar, daB | 3% | < Jy. 

Die Nichtlinearitat der Feldgleichungen macht sich im wesentlichen bei der 
Struktur der Eindringschicht bemerkbar. Nach der linearen Theorie hat man 
fiir die Stromverteilung im supraleitenden Halbraum (z > 0): 


as 


Ins 
OF = [ee-2leV2 ; [se RH: (5.3) 


1) Vgl. Koppr (1949, 1950). 
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Dabei ist H° das Magnetfeld an der Oberfliiche. LaBt man H®° wachsen, dann 
kann | $*| groBer als I, werden. Das geschieht vor allem bei tiefen Tempera- 
turen, da dort der Supraleiter einen hohen magnetischen Schwellenwert hat 
und J 4, fiir 7 =Onach Null gehen sollte. Andererseits mu8 wegen der MAXWELL- 


schen Gleichungen der gesamte Abschirmstrom, also das Integral iiber i ol 
gleich cH® sein. Das ist nur méglich, wenn der Suprastrom langsamer abfallt 
und damit das Magnetfeld tiefer in den Supraleiter eindringt. 

Bei der linearen Theorie bestimmt eine Messung der Eindringtiefe die Kon- | 
stante 2, und Messungen werden meist so ausgewertet, daB man den entspre- - 
chenden Wert von / berechnet. Auch im allgemeinen Fall kann man die Aus- | 
wertung in dieser Weise durchfiihren, und man erhalt dann einen scheinbaren | 
Wert Aer. Die oben betrachtete nichtlineare Theorie sollte dann dazu fiihren, 
daB Age von der Stirke des duBeren Magnetfeldes abhangt, in dem die Messung 
durchgefiihrt wird. Eine ins Einzelne gehende Diskussion findet sich bei, 
KOPPE (1950). ! 
Experimentell sind solehe Effekte von PIPPARD (1950b) aufgefunden worden. 
Die Abweichungen selbst sind allerdings im untersuchten Bereich sehr klein» 
und liegen wesentlich unter den von KOPPE angegebenen. Dabei ist aber zu: 
beriicksichtigen, daB den damals berechneten Werten der willkiirliche Ansatz 
(5.2) und eine sehr grobe Abschatzung der Maximalstromstiarke Jj, zugrunde: 
liegt. Es geniigt deshalb eine qualitative Ubereinstimmung. 

Die Moglichkeit, nichtlineare Effekte an Supraleitern durch die Erzeugung von: 
Oberschwingungen elektrischer Wellen (,,Klirrfaktor“‘) nachzuweisen, ist von’ 
HEYWANG (1950) theoretisch untersucht worden. 


6. Allgemeine Feldgleichungen | 


Es diirfte ziemlich wahrscheinlich sein, daB sich in einem Supraleiter eines 
GréBe o; als ,,Dichte der Supraleitungselektronen“* konstruieren lat, obwoh] j 
natiirlich eine Kinteilung in Supra- und Normalleitungselektronen im Hinbliek! 
auf die Ununterscheidbarkeit der Elektronen cum grano salis zu verstehen ist.’ 
Ferner ist anzunebmen, da8 sich ge, auf einen temperaturabhingigen Gleich-\ 
gewichtswert einstellt und da dadurch im wesentlichen die Temperatur-’ 
abhingigkeit von A zustandekommt. Wenn das richtig ist, dann brauchene 
aber A und g, nicht immer den Gleichgewichtswert zu haben, und das wiirde: 
bedeuten, dafS man gg statt als Konstante als FeldgréBe neben & einzufithrens 
hitte. Das ist von KoppE (1951) untersucht worden, wobei er sich zunichsti 
auf stationiire Zustinde beschrinkte. | 
Wir wollen dieses Problem jetzt allgemein behandeln, wobei wir uns wieder aut 
einen rein phinomenologischen Standpunkt stellen. Dazu nehmen wir an, dal! 
der Zustand eines Supraleiters nicht nur von G, sondern noch von einer Reihe| 
skalarer FeldgréBen 4; abhingt. — Uber die Bedeutung der 7; brauchen wit 
nichts zu wissen. In Gleichung (3.1) hatten wir diese Méglichkeit bereits vor: 
gemerkt. Wir kénnen jetzt wieder unmittelbar an (3.4b) ankniipfen und er} 
halten, wenn wir gleich noch aus (3.6) und (3.8) einsetzen: | 
02 . 

if grad V— 3 che Hm = div Y’+ To’. (6.1) 
| 

} 


a6 Zan, 
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Wir wollen zunichst versuchen, ohne eine Abinderung der ersten LONDON- 
schen Gleichung auszukommen, und setzen demzufolge zuniichst V = 0. Es ist 
offenbar unméglich, aus der linken Seite eine Divergenz herauszutrennen 
oder sie identisch zum Verschwinden zu bringen. Folglich mu8 gelten: 9{’’ = 0 
und 

0 


BOD comers 0. 
ON¢ 
Das kann man erfiillen durch 
. 02 0QY 
My Nk 


wenn die & Funktionen V7; der k Veranderlichen so beschaffen sind, daB 
> %M; (x,... te) = 9 (6.2) 


fiir alle Werte der 2; gilt. 
Man wird annehmen kénnen, dafi die Abweichungen der 4; von den Gleich- 


gewichtswerten und damit die a@ immer klein bleiben. In diesem Fall 


Oni 
kann man sich mit einer linearen Naiherung begniigen und 
: 0Q 
Hi = — SMe Mee (6.3) 


setzen. (6.2) fiihrt dann auf die bekannte Bedingung dafiir, daB die quadra- 
tische Form 2 y;.%;, x, nichtnegativ ist. 

Wir wollen uns auf den Ansatz (6.3) und ein einziges 77 beschranken. Wir haben 
dann als zusatzliche Gleichung 


: 0Q 
ay ee 6.4 
" Hay (6.4) 
und neben der JOULEschen Wiirme noch das zusitzliche Dissipationsglied 
02\ 
of Soe 6.5 
pape (Fn): ap" 


' Wir wollen kurz iiberlegen, was die beiden Gleichungen bedeuten, wenn wir fiir 
7 etwa 0s wihlen. Man kann die Entstehung der Supraleitungselektronen in 
einem Supraleiter als eine Art chemische Reaktion auffassen. (6.4) ist dann ein 
naheliegender Ansatz fiir die Reaktionsgeschwindigkeit in der Nahe des Gleich- 
gewichtes. Andererseits ist mit dem Ablauf einer Reaktion eine Entropie- 
erzeugung verbunden; es ist deshalb klar, da wir neben der JOULEschen Warme 

‘noch ein weiteres Dissipationsglied haben miissen. 

» Andererseits gibt gerade die vorstehende Uberlegung AnlaB, mit (6.4) noch 
nicht ganz zufrieden zu sein. Denn da der Suprastrom eine Bewegung der 

' Supraelektronen darstellt, sollte es zu lokalen zeitlichen Anderungen von Qs 
auch durch die Mitfihrung riumlicher Unterschiede kommen. Das heiBt, dak 

* man in (6.4) rechts noch ein Glied von der Form const. div 3% haben michte. 

Eine formale Méglichkeit dazu ergibt sich, wenn wir uns erinnern, daf (6.4) aus 

der speziellen Annahme V = 0 folgte. Wir ziehen jetzt die Méglichkeit V =--0 in 

- Betracht. Nach dem am Ende von Abschnitt 3 Gesagten mub V eine Funktion 


428 H. Koppr 
der Dissipationsglieder sein. BE? war ausgeschlossen worden. Wir kénnen aber | 
jetzt den Ansatz aQ i 
V ten te 7 
on = |i 
untersuchen, indem wir (6.4) zu | 
eet eae ) 


erweitern. Geht man damit in (6.1) ein, so kommt, wenn man noch die De-) 
finition von 9° beriicksichtigt, 


dQ GL « go os ee | 
$8 grad (“5 + yb (a) . an -F=div%A'+ To 


oder nach Umformung des ersten Gliedes: 


div = 3 + Ub (=) + te (F —x div 9*) = div YW’ + Fo". 
4 


Ov On 
Diese Gleichung ist erfiillt, wenn man setzt: 

i= Ge x div $8 (6.6) 

/ 
Sere oe ae (6.7 

OY 

002? 
T Be) pi less : 6.8 
o lu ( F) = (6.8: 


Dabei kénnen uw und x noch beliebige Funktionen der Feldvariablen sein. Die 
Gleichungen gelten unabhiingig von der modellmaibigen Bedeutung von 7. Sie 
miissen deshalb der Form nach invariant sein gegeniiber Substitutioner 
4-7 = 7 (mn). Das ist auch der Fall: Die Funktionen uw und x transformierer 


: ea al , , r / . ass 
sich dabei in w = (dy/dy’)? und x =x (dy /dy). Dagegen gilt “aa 

, 0 ‘1 
=x ~—,,d.h. der Arbeitsstrom und das Zusatzglied zur ersten LONDONScher} 


On 

Gleichung sind gegeniiber einer Anderung von 7 invariant. 
Uber die Bedeutung von 1’ und V wird man sich leicht klar, wenn man einmad 
annimmt, es sel 7 = 9s. Dann muB (6.6) fiir eine gehemmte Phasenumwandlung} 
(uw = 0) in die Kontinuititsgleichung tibergehen, daraus folgt, da man x = | 
m Z ne = a = ue im wesentlichen das 
On dos 00s 

chemische Potential der Supraleitungselektronen, und das Zusatzglied zu 
Lonponschen Gleichung ist eine Diffusionsscheinkraft, wie sie nach EINSTEII 
(1905) allgemein aus einem Gradienten des chemischen Potentials folgt. Da daa 
chemische Potential gespeicherte Arbeit darstellt, mu dann zwangsliufig ei 
der Mitfiihrung des Potentials entsprechender Arbeitsstrom existieren, ded 
durch (6.7) gegeben ist. Formal steht es natiirlich frei, x = 0 zu setzen; mad 
kommt dann auf (6.4) zuriick. | 
Hiir stationiire Zustiinde miissen die Dissipationsglieder verschwinden ; man ha 
aia ‘is (6.4 
y | 


mu setzen hat. Dann ist aber 


Allgemeine phinomenologische Theorie der Supraleitung 429 


wahrend in (6.6) alle Terme einzeln verschwinden. (6.9) ist die Verallgemei- 
nerung der friiher von KoppPe (1951) angegebenen Beziehung 


1 
5 (3—4f =9, 


ie sich aus dem Ansatz (4.1) mit 7 = A ergibt. 

6.9) ist eine Gleichung zwischen 7 und @, die dazu benutzt werden kann, 7 zu 
liminieren. Die Theorie fiir veriinderliche Parameter ist also fiir stationire 
Zustinde ein Sonderfall der allgemeinen nichtlinearen Theorie mit einem 
2 (G), welches durch Enveloppenbildung aus Q (G, 7) (als Funktionenschar 
mit dem Parameter 7 aufgefaft) gewonnen wird. Damit erledigen sich die von 
BECK (1952) erhobenen Einwinde. 


D 


7. Veranderlichhe Temperatur 


Wir haben von Abschnitt 3 ab nur Vorginge bei konstanter Temperatur be- 
rachtet, und uns daher auf die einfache Beziehung (2.4) beschrinken kénnen. 
Venn wir diese Voraussetzung jetzt fallen lassen, so ist 1. nachzupriifen, ob 
die bisher abgeleiteten Bewegungsgleichungen auch noch bei variabler Tem- 
peratur beibehalten werden kénnen, und 2. die Méglichkeit gegeben, galvano- 
hermische Effekte mit in Betracht zu ziehen. Beide Fragen sind einigermaBen 
nabhangig voneinander. Die Thermodynamik bedingt gewisse Beziehungen 
zwischen den galvanothermischen Effekten, aber nicht, deren Existenz. Wir 
wollen deshalb auf ihre Betrachtung verzichten, und lediglich die Konsistenz 
es bisher betrachteten Schemas beweisen. 

Die Gleichung (2.4) entsteht aus der allgemeineren Beziehung (2.3) unter der 
Annahme 7’ = const; d.h. wenn man alle Glieder wegliBt, die Ableitungen 
mach den Raumkoordinaten oder nach der Zeit enthalten. Das sind in un- 
ittelbar ersichtlicher Weise zunichst das zweite und dritte Glied in (2.3). 
% kann (das ist ein allgemeines thermodynamisches Prinzip) nicht von den 
Ableitungen von 7 abhingen. Dagegen kénnen Temperaturausgleichsvorginge 
zu Entropieerzeugung fiihren. Die Glieder 7'o in (2.3) und (2.4) diirfen sich also 
um Terme unterscheiden, die fiir 7’ = const. verschwinden. Um diese Ver- 
haltnisse deutlich zu machen, schreiben wir die beiden Gleichungen in der 
folgenden Form: 


(34) Sak lg aE i f O grad T+ 8 a 4+ To-+ To* =0 (7.1) 
s 


at aT dt 
Af 
(5) + div&%+ To =0. (7.2) 
Ot) 
Dabei bedeutet Ai die Zeitableitung von f bei konstantem 7’, und alle Sym- 
T 


|bole, die in beiden Gleichungen auftreten, haben jetzt exakt dieselbe 
| Bedeutung. 

| Nun sind in den vorangehenden Abschnitten die Bewegungsgleichungen so 
jeingerichtet worden, da (7.2) erfiillt ist. Subtrahieren wir (7.2) von (7.1), so 


aes 
- bleibt f pe 
Of Gay 1 
E i fg kate el ae as ora ™o* =(). 
(s+) ay 7 p © grad P+ To 


+ asain 
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Davon verschwindet die erste Klammer wegen ae = —8 und es verbleibt 
die Bedingung 5 grad 20! 73} 


Das hat aber (solange galvanothermische Effekte auBer Betracht bleiben) mit 
der Supraleitung gar nichts mehr zu tun, sondern ist eine Bedingung fiir de 
Mechanismus der Wirmeleitung. Macht man dafiir z. B. den Ansatz Q) = 
—{4 grad 7’, dann fiihrt (7.3) auf die Bedingung 4 > 0, also die Forderung, da 
die Wiirme von den wirmeren Stellen zu den kilteren flieBt. | 


8. Gleicdhgewichtsbedingungen 


Bisher haben wir stillschweigend vorausgesetzt, daB es im Raum ein festes Ge 
biet gibt, welches supraleitend ist, und innerhalb dessen demzufolge die bishe; 
abgeleiteten Gleichungen gelten. Da ein Metall auf Grund einer Phasenum 
wandlung supraleitend wird, muB natiirlich die Méglichkeit in Betracht ge 
zogen werden, daB diese Phasenumwandlung nicht tiberall eingetreten ist. Eh 
gebe also eine Phasengrenzfliiche T zwischen dem supraleitenden und norma. 
leitenden (d. h. noch nicht supraleitend gewordenen Material, die wir uns nae’ 
dem iiblichen thermodynamischen Schema als fiir Supraleitung semipermeablj 
Membran!) vorstellen kénnen. Die Gestalt von T ist nun ein weiteres Be 
stimmungsstiick des Zustandes des Supraleiters, und es entsteht die Frag¢ 
ob und wie sie sich im Laufe der Zeit andert. Wir suchen hier nur die Bedingunj 
dafiir, da sich T nicht andert, also Gleichgewicht besteht. Dazu rechnen wi 
die Krifte £ aus, die notwendig sind, um T festzuhalten. Gleichgewicht 
herrscht, wenn f iiberall verschwindet, oder wenn T mit der Oberfliche def 
Supraleiters iibereinstimmt, und f nach auBen gerichtet ist. | 
Wir stellen uns dazu vor, da T in einer vorgegebenen Bewegung begriffen isi 
und stellen die Leistungsbilanz fiir ein Gebiet auf, welches ein Stiick von if 
umfaBt2). 
Beziiglich der Grenzfliche ist zu bemerken, daf nur ihre Form physikalisel} 
Bedeutung hat, da eine Verschiebung der Fliche in sich selbst die physikalise} 
Situation nicht aindert. Wir kénnen deshalb von vornherein annehmen, dak dd 
Verschiebungsgeschwindigkeit » und die Kraft pro Oberflicheneinheit if 
immer senkrecht zur Oberfliche sind. £ ist also eigentlich ein Druck. » unc4 
sollen positiv gerechnet werden, wenn sie in das normalleitende Gebiet hineit 
welsen. 
Abb. 1 zeigt das betrachtete Gebiet B = B, +B, mit der Oberfliche C 

B, und B, sind die zur Zeit ¢ supraleitenden resp. nichtsupraleitenden Tei 


1) dh. ,,fiir Metall durchlissig, aber fiir Suprastréme undurchlissig‘*. Natiirlich ist o 
eine praktisch unmégliche Fiktion, aber das macht bekanntlich fiir thermodynamisc 
Betrachtungen nichts aus. 
*) Diese Methode wird (in unpriiziser Form) bereits von LONDON (]950) benutzt. Sie ]| 
den Vorteil, da man auf diese Weise unmittelbar zu der Beziehung (8.4) kommt. Man ka 
sich auch nach v. LAvE (1949) auf sehr langsame Anderungen beschriinken, die als rev | 
sibel betrachtet werden kénnen und dann nach dem iiblichen Schema die Krifte aus | 
Anderung geeigneter charakteristischer Funktionen berechnen. Das kann in gewis\ | 
Fallen vorteilhaft sein [z. B. bekommt man so die spaitere Gleichung (9.1) am einfachst¢a 
wird aber im allgemeinen dann doch recht kompliziert. Vel. auch BecK (19525). . 
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AB das in der Zeit At von T iiberstrichene Gebiet. Die Leistungsbilanz fiir B 
im Zeitabschnitt At lautet : 


—At/Xdo + At(k.v-do=| fit + At) dv—] f (t)dv + At|/ Todv, (8.1) 
fe} T B B B 

Links steht die B wihrend At zugefiihrte Arbeit, die entweder durch die Ober- 
fliche einstrémt, oder durch die Verschiebung von T geleistet worden ist. Sie 
wird entweder in Wiirme verwandelt, oder 
als freie Energie gespeichert. Wir bezeichnen 
fiir den Augenblick einmal mit /() resp. f() 
die Dichten der freien Energie (einschlieB- 
lich Feldanteil) im supraleitenden resp. nicht- 
supraleitenden Zustand. Es gilt also 


fe) = = (G+ §Y) (8.2) 


jm) — _ (G + §?). (8.3) 


Streng genommen hitten wir in beiden 
Gleichungen noch f, zu addieren; da es 
aber bei konstanter Temperatur eine Konstante ist, fallt es hier w ieder heraus. 
Wir haben nun: 

[f(t + At) dv —[ f(t) dv = | ff (t + At) —f (t)} dv 

B B By 
+ [ {7 (t + At) —f™ (dv + | f (t + At) )dv—| i”) (t) dv, 
B,—4B AB AB 
Wenn wir bis zu Gliedern von der Ordnung At entwickeln, und beriicksichtigen, 
daB AB und alle Integrale iiber 1B ebenfalls von der Ordnung At sind, so ergibt 
sich: 


“0 fl) 
fi (¢ + At) an— ff (‘jde=—2t | pi dv +A 7 we dv + At | Qed. 
8 B gai B; " 


n 
Das erste Glied links in (8.1) kann nach dem GAussschen Satz umgeformt 
_ werden. Man hat aber zu beachten, da { durch die Phasengrenzfliche 7 hin- 
durch nicht stetig zu sein braucht. Man mu} demnach den GAussschen Satz fiir 
B, und B, getrennt anwenden, und bekommt dann 


[ Xdo = [ div Adv + [ [MA], do 

fo) B ft 

wobei [%], der Sprung der Normalkomponente von QM ist. Setzt man das alles 
in (8.1) ein. so ergibt sich im Limes At— 0: 


Py r (at ) 
{(WU], -s kv —Qv} do = fy + div % + 73 dv. 
ot 
T B 
29 Zeitschrift ,,Fortschritte der Physik‘ 
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Die rechte Seite verschwindet identisch, und da die Gleichung fiir alle Gebiete | 
B gelten mub, folgt: Le 0p j 


9 setzt sich zusammen aus dem Poyntingvektor und eventuellen Zusatz- | 
gliedern nach (6.7). Der Poyntingvektor ist aber (wegen der Stetigkeit von ©} 
und §) stetig, und von dem Rest kann man leicht einsehen, daf [9], wenigstens | 
wie v? nach Null gehen muB. In der Grenze v = 0 hat man demnach 

Q= LATS: (8.4) ) 
Dabei gilt das <-Zeichen, falls der ganze Supraleiter supraleitend ist. Fiir| 
den urspriinglichen LoNDoNschen Ansatz gibt das die bekannte Bedingung 

1 


p Aneta f | 
ida oT (8.8)) 


9, Gleichgewichtsbedingungen fiir den didken Supraleiter 


Im letzten Abschnitt hatte sich ergeben, daB ein Supraleiter nur dann als 
Ganzes supraleitend sein kann, wenn der Supraimpuls (resp. Suprastrom) an 
der Oberfliche einen durch (8.4) gegebenen Maximalwert nicht iiberschreitet. 
Erfahrungsgemif kann beim dicken Supraleiter diese Bedingung durch die ¢ 
einfachere 


5H < Al (9.1) | 


ersetzt werden. Dabei gilt als ,,dicker‘* Supraleiter einer, dessen Kriimmungg 
iiberall sehr viel kleiner als ¢ VA ist, und der sich deshalb anniihernd wie ein! 
supraleitender Halbraum benimmt. 
Die Beziehung (9.1) laBt sich unmittelbar thermodynamisch ableiten!). Sie 
gilt unter der Voraussetzung, daB sich bei einer Verschiebung der Phasengrenz- | 
fliche der abschirmende Suprastrom als Ganzes mitverschiebt. Dann steht} 
links in.(9.1) die Arbeit, die vom elektrischen Feld geleistet wird (= MAXWELL-} 
sche Spannung X Weg), und rechts der Aufwand an freier Energie, der zur} 
Zerstorung der Supraleitung im iiberstrichenen Gebiet erforderlich ist. 
Wir wollen hier (9.1) unmittelbar an (8.4) anschlieBen. Dafiir ist zu zeigen, dahil} 
fiir die supraleitende Halbebene z > 0 bei z = 0 gilt: 


5 (0? = 2° [6 0)]. (9.1%) 
Fiir den stationiiren Zustand gelten die Gleichungen ' 
crot § + § =0 (9.21 
crot § —J* =0 (9.3) 
Of) 


SALE 9.4) 
0 y ( 4\ | 
Wir multiplizieren (9.2) vektoriell mit $*° und (9.3) vektoriell mit §. Subtrak~ 
tion der beiden Gleichungen gibt die Hilfsformel 


‘ 1 | 
[3° x rot G] = [ x rot ] = grad > $? — (© grad) ©. (9.559 


1) Vel. v. LAvE (1949). 
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Dann berechnen wir den Gradienten von 2 unter der Voraussetzung, daB Q 
von x, y, znur iiber © und 7 abhiingt!). 


grad Q = (%* grad) G + [9* x rot G] + ae grad 7 (9.6) 
oder, mit (9.4) und (9.5): i 
il 
grad ( ary s*) = (98 grad) G6 — ( grad) H. (9.7) 


Nun kénnen beim Halbraum alle GréBen nur pumoner von z sein. Wegen 
div 3* = Ound div § = Osind die z-Komponenten von 9° und § konstant und 
also Null, da sie fiir z = 0 verschwinden. Dann kommen ies in (9.7) rechts nur 
Differentiationen in Richtung der Phasengrenzfliche vor, die Null ergeben. 
Folglich ist: l l 
dante 9? chads a Al ace G*= const. 


Die Konstante laBt sich bestimmen, wenn man mit z nach oo geht. Dann 
gehen (und damit 2°) und § nach Null, und es bleibt const = — Af oder 


° Lops 
@ => § 


fiir z = 0 gibt das die zu beweisende Gleichung (9.1). 

Mit Riicksicht auf Spateres wollen wir noch eine andere Interpretation von 
(9.6) geben. Unter der Voraussetzung, daB alles nur von z abhangt, und den 
daraus gezogenen Folgerungen wird aus (9.6): 


dQ 
dz 
und nach Integration tiber z von 0 bis co: 


2 (0) —2 (co) = 2° (0) = |= [3° x H] de. 


l 
= [3° Xrot G] = — — [8° x 9] 


Ganz rechts steht hier die Kraft, die vom Magnetfeld auf die Suprastréme 
ausgeiibt wird, und die demnach gleich dem Wert von 2° fiir z = 0 ist. Diese 
Kraft erscheint also an der Phasengrenzfliiche als ein Zug. der nach (8.4) vom 

,.Phasenumwandlungsdruck** — Af im Gleichgewicht gehalten resp. an die 
Oberfliche des Supraleiters iibertragen wird. Das gilt, wie nochmals hervor- 
gehoben sei, nur, solange 2 nicht explizit von den Raumloordinaten abhinet. 
Sobald das der Fall ist, treten Volumkrifte auf die Materie auf, die spiiter 
(Abschn. 11) betrachtet werden sollen. 


10. Eindeutigkeitssatz 


. Eine wesentliche Forderung, die an eine phiinomenologische Theorie der Supra- 
leitung zu stellen ist, ist die, daB der stationiire Zustand, der sich in einem 
einfach zusammenhingenden Supraleiterin einem konstanten iuBeren Magnet- 
feld nach einiger Zeit einstellt, unabhiingig von der Vorgeschichte ist. Das ist 
nicht trivial, da ja z. B. Gleichung (1.1) offenbar aussagt, daf} die Elektronen 


1) Diese Voraussetzung ist wesentlich; fiir einen etwa infolge veriinderlicher chemischer 
Zusammensetzung inhomogenen Supraleiter gilt (9.1) nicht! 
29* 
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durch ein elektrisches Feld beschleunigt werden. Nun treten beim Kinschalten 
eines Magnetfeldes elektrische Feldstirken auf, die im einzelnen von der Art 
des Einschaltvorganges abhingen. Es ist dann nicht ohne weiteres einzusehen, | 
daB diese verschiedenen Beschleunigungen immer zu den gleichen End- 
geschwindigkeiten fiihren. In der Tat ist das auch bei ungeschickten Abiinde- 
rungen von (1.1) nicht mehr der Fall. 

Da der Zustand eines Supraleiters im konstanten Magnetfeld von der Vor- 
geschichte unabhangig ist, kann man beweisen, indem man die Kindeutigkeit 
des stationiiren Zustands beweist. Dieser mu& sich schlieBlich einstellen, da alle 
Anderungen mit Energiedissipation verbunden sind, und daher schlieSlich zum 
Stillstand kommen miissen. Der Eindeutigkeitsbeweis ist fiir die lineare Theorie 
von Vv. LAUE (1949, § 12) und fiir die nichtlineare Theorie von BECK (1951) ge- 
geben worden. Wie bereits in Abschnitt 6 bemerkt wurde, geniigt das, da die 
allgemeinste Form fiir den stationiren Zustand auf die nichtlineare Theorie’ 
zuriickfiihrt. | 
Die allgemeinen Beweise sind recht kompliziert; wir beschranken uns deshalb 
hier auf den einfachen Fall des zylindrischen Supraleiters im homogenen ma-. 
gnetischen Feld (etwa einer Spule). Wir haben dann im Innern des Supraleitets 


crotG + © =—0 (10.1) 
und die MAXWELLsche Gleichung 
TOR Sia or (10.2) 


Aus (10.1) folgt als Randbedingung 
: igue 
rot & = — e §° (auf dem Rand) (10.3). 


wobei $° ein vorgegebener konstanter Vektor in Richtung der Zylinderachs 
ist. Aus (10.1) und (10.2) folgt 


s 1 g ( 
rot rot ( +- 7 (Gees On, (10.4) 


Wir nehmen nun an, es giibe zwei Lisungen (, und @, zu der gleichen Rand4 
bedingung. Ks ist also auf dem Rand rot (%, — @,) = 0. 
Wir schreiben nun die Gleichung (10.4) fiir G, resp. G, an, subtrahieren, unc 
multiplizieren skalarmit (®, — ,): 

c? (G, — G,) rot rot (G, — G,) + (G, — G,) (B§ — Ys) = 0. 
Das li Bt sich umformen zu: 
c | rot (6, — G,) |? + (G, — Gy) (3 — 3) =e div {(G, — G,) rot (G, — (,)} H 


Integriert man diese Gleichung iiber den ganzen Querschnitt des Supraleiters: 
so erhalt man: | 


[{e | xot (G — G) 2 + G,— G,) (98 —39} do= 0. (10.55 
Setzt man nun voraus, da’ die Funktion Q so beschaffen ist, daB immer 
(®, — Gs) (9* [G,] —F° [Gy] = 0 (10.6) 


gilt, dann ist der Integrand nichtnegativ, und aus (10.5) folgt &, — G, = Of 
also die Hindeutigkeit der Losung, / | 
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e Bedingung (10.6) ist hinreichend ; es konnte nicht bewiesen werden, daB sie 
otwendig ist. Fiir die lineare Theorie gibt sie die Bedingung 2 > 0. DaB A 
ositiv sein mub, damit die Stromverteilung im Supraleiter wenigstens quali- 
tiv richtig ist, ist ohne weiteres klar, Es laBt sich aber in diesem Fall auch 
icht der Zusammenhang mit der EKindeutigkeit zeigen. Fiir Q° = &2/24 mit 
onstantem / ergeben sich aus (10.4) drei Gleichungen der Form 
Ap+Ey=0; H=—1/ed. 
as ist aber ein bekanntes Eigenwertproblem. Damit die Lésung der Randwert- 
uufgabe eindeutig ist, darf E kein EKigenwert sein. Da alle Eigenwerte positiv 
d, und fiir geniigend groBe Querschnitte beliebig nahe an Null heran- 
ebracht werden kénnen, ist die notwendige und hinreichende Bedingung dafiir 
erade A > 0. 
die tensorielle Theorie liefert (10.6) die bereits von v. LAUE angegebenen 
terminantenbedingungen fiir die A; ,. 
tir den mehrfachzusammenhiingenden Supraleiter hat LONDON einen Satz 
ber das .,Hangenbleiben** des Magnetfeldes in Léchern bewiesen. Fiir den 
eweis') benétigt man nur die Gleichung (3.7), er kann deshalb unverindert 
iir die allgemeine Theorie ttbernommen werden. 


11. Krafte der Suprastréme auf die Materie 


e bisher entwickelte Theorie ist vollstiindig, soweit man sich nicht um die 
églichen Riickwirkungen des Supraleitungsmechanismus auf die Materie 
iimmert. Diese trat ja in den vorangehenden Abschnitten lediglich als 
..Raumgebiet, in dem gewisse Gleichungen g gelten**, auf. Wir hatten schon im 
Abschnitt 9 gesehen, daB die Kriifte, die z. B. ein Magnetfeld auf einen 
Suprastrom ausiibt, in einem homogenen Supraleiter als Zugspannungen an die 
Phasengrenzfliiche iibertr agen werden, die ihrerseits durch den Umwandlungs- 
druck an der Oberfliche des Supr aleiters festgeklebt ist. Wenn man also etwa 
einen zylindrischen Supraleiter in ein longitudinales Feld bringt, dann wirkt auf 
die Oberflache des Zylinders ein allseitiger Druck, und der Zylinder wird etwas 
zusammen gedriickt ‘werden. Natiirlich sind diese Krifte im Verhaltnis zur 
Kompressibilitit der Metalle so klein, da man sich praktisch um diese Defor- 
mationen nicht zu kiimmern braucht. Eine geschlossene Theorie muB aber im- 
stande sein, sie zu berechnen. 

Die Krafte, die an der Phasengrenzfliiche auftreten, sind bereits in Abschnitt 8 
berechnet worden: man erhilt Zugspannungen von der Grobe 2° (%). Am Ende 
von Abschnitt 9 ist aber bereits darauf hingewiesen worden, daB daneben auch 
Volumkriifte auftreten kénnen, die wir jetzt berechnen wollen. 

Wir gehen wieder von einer Leistungsbilanz aus, indem wir die Materie des 
Supraleiters mit in die Grundgleichung (2.4) hineinnehmen?). Die Materie habe 
eine Strémungsgeschwindigkeit » und unterliege den Volumkraften £. Dann 
‘ist der Energiesatz (2.1) zu erweitern zu: 

Ju 


— Ot 
7) Vel. v. LAUVE (1949, § 12). 
2) Das ganze Verfahren ist eine Nachbildung der Berechnung von Kraftwirkungen auf ein 
Dielektrikum; vgl. R. BECKER (1941). 


-+ div (% + 0) + vf=0. 


! 
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Daraus folgt fiir isotherme Vorgiinge ganz wie im Abschnitt 2: P| 
q 
4 div% + To +vf=0. (11.1) 


men, wenn der Supraleitungsmechanismus vom Zustand des Supraleiters ab{ 
hiingt. Betrachten wir als Beispiel den einfachen Ansatz: 
yi - t 


° 1 ORs Cys2 


Nun kann es zu einer Kopplung zwischen Supraleitung und Materie nur <a 


A hangt zunichst einmal von 7’ ab, was hier nicht weiter interessiert. Es ist abe 
ferner bekannt, daB der Sprungpunkt eines Supraleiters von Druck (d. h. als 
der Dichteo ,) und der chemischen Zusammensetzung, gegeben durch gewiss¢ 
Parameter y,; (Konzentrationen) abhingt. Das gleiche ist dann von / zu er) 
warten. Gehen wir gleich wieder zum allgemeinsten Ansatz zuriick, so ist jetzy 
in Betracht zu ziehen, daB 2 auBer von & und den 7; noch von materia 
gebundenen Parametern gy und y; abhiingt. Da eine Dichtedinderung nur ein 
spezielle Deformation darstellt, kénnen wir die Liste gleich noch durch dil 
Komponenten ¢;, des Deformationstensors erginzen. Materiegebunden sol 
dabei heifen, da diese GréBen keine selbstiindigen Bewegungsgleichunge: 
haben, sondern durch das Verhalten der Materie bestimmt werden. Sie ve 
halten sich dabei etwas verschieden und kénnen zweckmaBig in Zusamme: 
setzungs- und Deformationsparameter unterteilt werden. 

Als Beispiel der letzteren behandeln wir zuniichst nur gj. Es gilt offenbar 
Kontinuititsgleichung 


dom 


Ot 


= —divoy » = — oy div » — bv gradoy. (11. 


Dagegen wird die chemische Zusammensetzung bei der Bewegung einfae 
mitgefiihrt (d.h. der substantielle Differentialquotient von y; ist gleich Null } 
und es gilt: | 

>, = —bgrad’y;. (11L.e] 


Das erste Glied von (11.2) fehlt hier, 

Ks sei nochmals darauf hingewiesen, da unter f immer die Differenz zwische 
freier Energie und ihrem Wert im nichtsupraleitenden feldfreien Zustand ver 
standen ist. Diese Definition erweist sich hier wieder als zweckmibig, weil s! 
die Kriifte £, die wir berechnen, die elastischen Spannungen, die auch ohn 
Supraleitung bei einer Deformation auftreten wiirden, nicht enthalten. 
Wir gehen nun auf (11.1) zuriick, und berechnen zuniichst 


Of a 00% dQ. 
a est)| +>=—owm+ J %- 
Ot (55 oy ”, Jom M pa OV; Vi 


Dabei gibt der erste Term rechts die Anderung von f bei konstantem Ow wt 
Vir welche in den vorangehenden Abschnitten bereits in eine Summe’ var 
Divergenzen resp. Dissipationsgliedern umgeformt worden ist. Bezeichnen wi 

| 


| 
| 
| 


| 
| 
| 
| 


Allgemeine phiinomenologische Theorie der Supraleitung 437 


en eventuell noch neu auftretenden Arbeitsstrom mit 9’, so kann man an 
lle von (11.1) einfach schreiben: 

dQ. ,0Q . 

aS ~ ay, y+ div UW + tb =0. 


GM 


zt man hier aus (11.2) und (11.3) ein, so ergibt sich: 


dQ 0 
fv = —— divoyv + > see v grad y; — div Qt’ 
Oom OY; 
{ dQ 
=D _— Ou grad ton t= pe seat bah pee ete 
( Jox i J Joa 
araus folet W =vow | 
(11.5) 
0Q dQ 
£ = — rad ee Se dy,. 
Om gra dou + ao! OV; grad Vy; 


Af. 


sei zuniichst 4 eine Funktion der chemischen Zusammensetzung. Da jede 
ktion der Konzentrationen auch der Gleichung (11.3) geniigt, kOnnen wir 
leich 2 als eines der y; behandeln. Wir haben dann: 


ep”. ; 1 

(= aA S 6) grad 4 = — — 3 grad A. 
Es treten demnach Volumenkrafte an den Stellen auf, an denen sich die Zu- 
sammensetzung der Materie iindert. Daf das so sein mu, macht man sich am 
einfachen Beispiel des supraleitenden Halbraumes klar, bei dem die Zusammen- 
setzung eine Funktion von z sei, wihrend wir uns die Phasengrenzfliche fest- 
gehalten denken. Dann fiihrt eine Verschiebung des Materials bei festgehaltener 
Phasengrenzfliiche und konstantem ‘uf eren Feld zu einer Anderung der 
Verteilung des Suprastromes, in der Abschirmschicht, und damit zu einer 
Anderung der gesamten freien Energie (resp. der fiir den speziellen Fall zu- 
stiindigen charakteristischen Funktion). Kine solche Anderung bei einer 
Verschiebung bedeutet aber eine Kraft in der V erschiebungsrichtung. 
Wir nehmen nun an, daf das Material homogen ist, da aber A von der Dichte ox 
abhingt. Wir haben dann: 
f = oy grad (5) ~ 3°). 


) 2 dou 
; 

In diesem Falle tritt, wie man sich leicht iiberzeugt, keine resultierende Ge- 
‘samtkraft auf, sondern das supraleitende Material wird in den Gebieten de- 
formiert, in denen ein Suprastrom flieBt. 
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Wir wollen noch die Kriifte bei allgemeinen Deformationen angeben. Wenn ein 
materieller Punkt die Verschiebung ¢; erfahrt, dann ist : 


~ “whew a 


Lee (Cen + Cr, e) 


Od; ] 
und wegen 0¢; = v,dt a = 5 (Vie + Vk, i) 


Nach dem gleichen Schema wie oben ergibt sich dann: 
vk 


Wenn man unter f die gesamte freie Energie verstanden hatte, dann enthielte 
die potentielle elastische Energie des Metalles, und der entsprechende Antei» 
von (11.6) ware dann die bekannte Kraftdichte, die aus den elastischen Span 
nungen resultiert. 


Géttingen, Max-Planck-Institut fiir Physik. J 
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1954 besteht das 


CHEMISCHE ZENTRALBLATT 
125 Jahre 


Es ist uns eine Freude, darauf hinweisen zu kénnen, daB der bereits 

im 4. Quartal 1953 um etwa 1600 Seiten vermehrte Textumfang des 

.-Chemischen Zentralblatts** im Jubiliumsjahrgang nicht nur bei- 
behalten, sondern noch weiter erhéht werden wird. 


Der Umfang wird im Jahre 1954 
rund 11200 Seiten Text 


(zuziiglich rund 1500 Seiten Autoren- und Patentregister) 


gegeniiber rund 9600 Seiten Text (zuziiglich rund 1100 Seiten 
Autoren- und Patentregister) im Jahre 1953 und 8080 Seiten Text 
(zuziiglich 971 Seiten Autoren- und Patentregister) im Jahre 1952 


umfassen. 


Der Jubilaumsjahrgang wird also den Umfang der Vorkriegsjahr- 
ginge iibertreffen und damit der umfangreichste Jahrgang des 
,.Chemischen Zentralblatts‘** seit seinem Bestehen sein. 


Der erhéhte Umfang gewahrleistet die Vollstandigkeit der Bericht- 
erstattung, wobei trotz des erweiterten Umfangs auf die Qualitat 
der Referate besonderer Wert gelegt wird. 


Eine geringe Preiserhéhung lieB sich dabei nicht vermeiden. Sie betragt gegeniiber 

den Vorjahren DM 50,—, d.h. rund 12%, wahrend die Umfangserweiterung 40% 

ausmacht. Der Verkaufspreis pro Druckbogen (= 16 Seiten), der bis 1952 DM 0,57 
betrug, verringert sich dadurch auf DM 0,47. 


Der Jahresbezugspreis betragt fiir 1954 also DM 370,— einschlieBlich Autoren- und 
Patentregister; das Ende 1955 erscheinende Sach-und Formelregister zum Jahrgang 
1954 wird, ebenso wie das Sach- und Formelregister zu 1953, DM 96,— kosten. 
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TECHNISCHES ZENTRALBLATT 


Im Auftrage der Deutschen Akademie der Wissenschaften zu Berlin 


herausgegeben von 


PROF. DR. MAXIMILIAN PFLUCKE 


ABTEILUNG MASCHINENWESEN 


Die Abteilung Maschinenwesen vermittelt den Stoff nach folgender Gliederung: 


A. Allgemeines: I. Kraftmaschinen 
Geschichte, Unterricht, Forschung - 
Grundlagenwissenschaften K. Verkehrstechnik, ihre Fahrzeuge, 


Maschinen und Einrichtungen 


L. Feinwerktechnik, Optische Industrie, 


B. Fabrikbetrieb 


C. Werkstoffe, Hilfsstoffe, Betriebsstoffe Medizinische Technik 
' des Maschinenwesens und deren Prii- 
fung M. Maschinen und Gerite fiir die Land-, 
Garten- und Forstwirtschaft, Vieh- 
D. MeBtechnik, Werkstatt-Priiftechnik . ’ 
‘und Regelungstechnik im Maschinen- wirtschaft und Fischereiwesen 
wean N. Anwendungen der Maschinentechnik 


in anderen Industriezweigen: Energie- 


E. Maschinenelemente, ihre Theorie und und Wasserverniorgung > *Cleackeian 


oxtiguns: Technik - Textilfaser-, Textil- und Be-- 
¥. Bearbeitungsverfahren, Werkzeug- kleidungsindustrie - Bauwesen - Berg- 
maschinen, Werkzeuge und Vorrich- und Hittenwesen - sonstige Industrie- 
tungen. zweige 
G. Arbeitsmaschinen O. Unfallschutz - Ausstellungen und Ta- 
gungen - Gewerblicher Rechtsschutz - 
H. Heizungs-, Liiftungs- und Klima- Bibliographien - Buchbesprechungen « 
anlagen, Trockenanlagen, fen Patente 


Als Informationsorgan wird das ,,Technische Zentralblatt — Abteilung Maschinenwesen* 
besonders wertvoll sein fiir: 


Werkzeugmaschinenbauer, Ingenieure des Kraftwerk- ind Kraftmaschinenbaues, Ingenieure 

des Kisenbahn-, Automobil-, Wasser- und Luftfahrzeugbaues, der Feinwerk-, MeB8- und 

Regelungstechnik, Konstrukteure fiir Férderanlagen, Verarbeitungsmaschinen, Heizungs- 

und Liiftungsanlagen, Landmaschinenbauer, ferner Fachleute auf den Gebieten des Fabrik- 
betriebes, der Wairmetechnik, der Strémungslehre u. a. m. 


Gleicherweise werden die Spezialisten in don Entwicklungs-, Konstruktions- und Fertigungs- 

biiros, die schaffenden Krifte in den Betrieben, die Lehrer an den verschiedenen Fach- 

schulen, die Fachleute in Verwaltung und Industrieverbinden, die Leiter und Benutzer von 

Fachbiichereien und beratende Ingenieure aller Sparten wertvolle und vollstiindige Unter- 
richtung tiber ihre Fachgebiete erfahren. 


635 Fachzeitschriften des In- und Auslandes werden zur Zeit (Januar 1954) im Technischen 
Zentralblatt, Abteilung Maschinenwesen, referiert. Die Anzahl wird laufend erhéht. 
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